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SECTION  I 
INTRODUCTION 


1.  OBJECTIVE 

The  primary  objective  of  the  Air  Force  Refuse-Collection  Scheduling  Pro- 
gram is  to  reduce  refuse-collection  costs.  Program  PHASE2  accomplishes  the 
primary  reduction,  that  of  manpower  costs,  by  assigning  streets  to  collection 
sections  in  such  a way  as  to  produce  the  minimum  number  of  trips.  Program 
PHASE3  further  reduces  costs  by  generating  efficient  collection  routes,  which 
may  also  be  minimum-length  routes,  for  each  section. 


2.  SCOPE 

This  section  (Volume  III)  of  the  report  describes  the  workings  of  the 
third  program,  PHASE3.  A program  overview  is  given,  followed  by  a thorough 
description  of  the  logic  involved  in  route  generation.  A skeleton  of  the 
logic  flow  is  provided.  Input  and  output  files  are  described.  Program  re- 
quirements and  restrictions,  error  messages  and  error  handling  techniques, 
definitions  of  important  variables,  and  an  estimate  of  running  time  are  also 
provided . 
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SECTION  II 
PROGRAM  OVERVIEW 


The  purpose  of  program  PHASE3  is  to  produce  routes  that  fully  describe 
the  path  each  collection  vehicle  will  travel.  The  distance  traveled  by  each 
vehicle  to  service  a section  will  be  approximately  the  minimum  possible  dis- 
tance, subject  to  limits  on  U-turns  and  the  original  selection  of  segments  to 
be  included  in  the  section  by  program  PHASE2.  Two  trips  are  produced  for  each 
section,  one  starting  at  the  garage  and  one  starting  at  the  landfill,  enabling 
program  PHASE4  to  further  reduce  the  overall  travel  distance  by  an  appropriate 
selection  of  trips. 

Input  to  program  PHASE3  consists  of  card  data  and  four  disk  filrs.  The 
card  data  must  include  a title  card  and  a card  giving  the  node  numbers  of  the 
landfill  and  garage  and  may  also  include  up  to  500  cards  that  reassign  seg- 
ments to  different  sections.  The  disk  files  contain  segment  and  node  data  for 
the  entire  map,  a list  of  segments  ordered  by  section,  a list  of  vehicle  capa- 
cities, and  pointers  to  the  first  and  last  segments  in  each  section. 

The  output  produced  by  PHASE3  consists  of  printed  output  and  one  disk 
file.  The  disk  file  contains  segment-by-segment  descriptions  of  the  paths  ser- 
vicing each  section.  These  data  will  be  used  by  program  PHASE4  to  generate  the 
route  maps  and  printed  schedules.  The  printed  output  gives  the  initial  and 
final  segment  assignments  and  vehicle  capacities  and  a limited  amount  of  data 
to  be  used  for  debugging  problems  that  may  arise. 

The  program  consists  of  a main  program  named  PHASE3  and  25  subroutines. 

The  flow  of  control  from  one  subprogram  to  another  is  displayed  in  Figure  1. 
Normal  processing  is  shown.  Within  each  subprogram,  only  the  first  call  to 
each  other  subprogram  is  Indicated.  (Subroutine  EXIT  is  a system  subroutine 
and  is  not  included  in  the  description  of  program  PHASE3. ) 

In  its  first  stage,  program  PHASE3  reads  the  problem  title  from  the  first 
data  card  and  the  segments  and  their  section  assignments  from  disk  files.  The 
node  data  are  also  read  from  a disk  file.  The  vehicle  load  is  computed  for 
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Figure  1.  Control  Relationships  Among  Subprograms 
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Figure  1.  Control  Relationships  Among  Subprograms  (Concluded) 


each  section.  Subroutine  ADJUST  Is  called  to  put  Into  the  same  section  all  seg- 
ments that  form  parts  of  the  same  street.  This  reassignment  Is  performed  only 
when  a vehicle's  capacity  will  not  be  exceeded.  Segments  that  do  not  require 
collection  are  assigned  to  section  zero.  PHASE 3 reassigns  additional  segments 
as  user-specified  section  changes  are  read  from  data  cards.  The  changes  and 
final  segment  assignments  are  printed. 

Connected  sections  are  formed  In  the  second  stage  of  PHASE3  processing. 
Subroutine  TREE  is  called  three  times  to  generate  minimum-time  paths  to  all 
nodes  for  travel  away  from  the  garage  and  travel  to  and  from  the  landfill. 

The  path  from  the  landfill  to  the  garage  is  saved  on  file  TAPE9.  Subroutine 
CONNECT  is  called  to  connect  the  pieces  of  each  section.  When  a section  con- 
sists of  more  than  one  piece,  subroutine  CONNECT  searches  for  single-segment 
paths  connecting  the  pieces.  If  pieces  remain  that  cannot  be  connected  by 
single-element  paths,  subroutine  C0N2ST  is  called  to  connect  pieces  with  two- 
step  paths.  Subroutine  CONNST  is  called  to  connect  any  remaining  pieces. 

PHASE3  calls  subroutine  DISCON  to  disconnect  dangling  segments  on  which 
refuse-collection  is  not  required. 

Paths  servicing  the  sections  are  generated  in  the  third  stage  of  PHASE3 
processing.  From  1 to  16  paths  are  generated  for  each  of  two  starting  points: 
the  landfill  and  the  garage.  All  of  the  paths  end  at  the  landfill.  The  short- 
est path  from  each  starting  point  is  saved.  PHASE3  calls  subroutine  EPXP  to  se- 
lect good  entry  points  and  exit  points  for  each  section.  For  each  choice  of 
entry  and  exit  points,  subroutine  CL0SE1  is  called  to  find  return  paths  to  the 
section  from  dangling  or  dead-end  segments.  Subroutine  CL0SE1  allows  U-turns 
in  some  cases  and  adds  additional  segments  to  the  section  in  other  cases. 

Finding  the  travel  path  in  the  collection  region  requires  pairing  the  odd 
nodes  (nodes  at  which  an  odd  number  of  segments  meet).  The  paths  pairing  the 
odd  nodes  must  be  traveled  a second  time.  If  the  pairing  produces  the  minimum 
total  pairing  distance,  and  if  collection  is  required  on  all  segments  in  the 
region,  then  the  travel  path  in  the  collection  region  will  be  of  the  minimum 
possible  length.  Some  streets  in  the  section,  such  as  streets  added  by  subrou- 
tine CONNECT  to  connect  pieces  of  the  section,  may  not  require  collection. 

These  streets  may  cause  the  paths  generated  using  a minimal  odd-node  pairing 
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scheme  to  be  slightly  longer  than  the  minimum  that  could  be  obtained  if  streets 
not  requiring  collection  were  omitted. 

PHASE3  calls  subroutine  GENDM  to  generate  a matrix  of  distances  between 
each  pair  of  odd  nodes.  Subroutine  SELORD  is  called  to  pair  the  odd  nodes  in 
a near-minimum  manner.  Subroutine  SOLV  is  called  to  obtain  an  odd-node  pairing 
with  the  minimum  possible  pairing  length.  PHASE3  computes  the  number  of  times 
each  segment  in  the  section  must  be  traversed. 

Subroutine  TRAVEL  is  called  to  find,  by  trial  and  error,  a path  from  the 
garage  or  landfill  through  the  collection  region  and  then  to  the  landfill. 
Subroutine  OPTPATH  is  called  to  perform  some  improvements  on  the  path.  The 
two  shortest  paths  for  each  section,  one  starting  at  the  garage  and  one  start- 
ing at  the  landfill,  are  written  to  file  TAPE9. 
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SECTION  III 
PROGRAM  LOGIC 

The  logic  for  program  PHASE3  is  described  from  three  viewpoints.  The 
first  description  is  task  oriented.  The  second  is  data-storage  oriented  and 
includes  discussions  of  the  preparation  of  data  for  use  by  program  PHASE4,  the 
use  of  input  data,  and  the  data  structures  used  in  PHASE3.  The  third  view  de- 
scribes each  subroutine  and  the  main  program  in  terms  of  its  purpose  and  the 
manipulations  performed  within  it. 

1.  PROGRAM  TASKS 

The  processing  performed  by  PHASE3  can  be  grouped  into  four  task: : seg- 
ment adjustments,  section  connections,  odd-node  pairing,  and  path  generation 
and  optimization.  The  data  used  during  segment  adjustments  come  from  four 
files.  The  segment  numbers  on  file  TAPE4  are  in  order  by  section.  File  TAPE3 
contains  pointers  marking  the  first  and  last  segments  in  each  section.  All 
remaining  data  for  the  segments  come  from  file  TAPE1 . Node  data  from  file 
TAPE2  are  also  needed  during  the  program-initiated  section  reassignments  of 
the  segments.  The  number  of  segments  meeting  at  each  node  (the  order  of  the 
node)  is  computed  for  each  node  from  the  packed  word  of  bounding-segment 
numbers. 

The  program- i ni tiated  changes  are  performed  by  subroutine  ADJUST.  This 
subroutine  examines  the  segments  bounding  each  order-two  node.  If  both  of  the 
segments  contain  refuse  and  each  is  serviced  in  a different  section,  the  sub- 
routine changes  the  section  assignment  of  one  of  the  segments  (if  the  change 
will  not  cause  the  vehicle  capacity  to  be  exceeded).  The  reassignments  are 
printed  by  the  subroutine  as  they  are  made.  The  amount  of  refuse  in  each  sec- 
tion is  recomputed  and  printed  after  control  returns  to  the  main  program. 

The  user-initiated  section  reassignments  are  controlled  by  data  cards 
read  by  the  main  program.  The  segment  number  and  its  new  section  assignment 
are  read  from  a data  card.  The  segment  is  reassigned,  and  the  corresponding 
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refuse  quantity  Is  transferred  from  the  old  section  to  the  new  section.  The 
program  makes  no  attempt  to  keep  the  refuse  quantity  In  any  section  from  ex- 
ceeding the  vehicle  capacity;  this  Is  the  user's  responsibility.  A message 
is  printed  describing  each  segment  reassignment.  The  final  refuse  quantities 
for  the  sections  are  printed. 

The  first  step  In  the  formation  of  connected  sections  Is  to  verify  that 
each  node  Is  connected  to  the  landfill  and  the  garage.  Subroutine  TREE  is 
called  three  times  to  generate  paths  from  the  landfill  to  each  node  and  from 
the  garage  to  each  node.  The  paths  between  node  and  landfill  are  generated 
twice,  once  for  travel  toward  the  landfill  and  once  for  travel  away  from  the 
landfill.  The  paths  are  built  to  minimize  travel  time.  The  path  Information 
consists  of  an  array  that  gives  the  line  number  of  the  next  node  on  the  path 
to  the  landfill  or  garage  and  an  array  that  gives  the  next  segment  In  the  path 
to  the  landfill  or  garage.  If,  at  any  node  other  than  that  of  the  landfill  or 
garage,  the  line  number  of  the  next  node  in  the  path  Is  zero  or  the  number  of 
the  next  segment  In  the  path  Is  zero,  then  that  node  Is  not  connected  to  the 
landfill  or  garage.  Error  messages  are  printed  when  unconnected  nodes  are 
found. 

Subroutine  CONNECT  Is  called  once  for  each  section  to  connect  separate 
pieces  of  the  section.  The  subroutine  examines  each  segment  In  the  map  de- 
scription until  it  finds  a segment  that  is  part  of  the  desired  section.  The 
starting  node  of  this  segment  is  used  as  the  starting  point  of  a minimum- 
distance  tree,  but  only  segments  within  the  section  are  used  in  building  the 
paths.  The  remaining  segments  in  the  map  description  are  examined.  If  a 
segment  is  In  the  section,  but  one  of  its  nodes  is  not  connected  to  the  path, 
a separate  piece  of  the  section  has  been  found.  Paths  are  built  from  nodes 
in  each  subsequent  piece  to  find  all  of  the  segments  in  each  piece.  When  the 
segments  have  been  marked  according  to  the  piece  in  which  they  occur,  subrou- 
tine CONNECT  adds  to  the  section,  for  travel  purposes  only,  each  single  seg- 
ment that,  connects  two  pieces  of  the  section. 

If  the  section  still  consists  of  more  than  one  piece,  subroutine  C0N2ST 
is  called  to  find  the  best  two-step  paths  between  the  remaining  pieces.  If  the 
section  then  still  has  more  than  one  piece,  subroutine  CONNST  is  called  to  find 
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the  shortest  path  connecting  two  pieces.  These  pieces  are  connected.  If  more 
than  one  piece  yet  remains,  the  next  shortest  path  between  two  pieces  is  found. 
The  process  continues  until  all  pieces  of  the  section  have  been  connected. 

The  next  step  in  the  formation  of  connected  sections  is  the  deletion  of 
dangling  or  dead-end  segments  on  which  no  collection  is  required.  The  main 
program  calls  subroutine  DISCON  to  disconnect  these  segments. 

’ 

The  path  through  the  collection  region  usually  requires  that  some  streets 
be  traversed  more  than  once.  The  path  must  use  two  segments,  or  one  segment 
twice,  each  time  it  crosses  a node.  At  an  even-order  node,  each  street  can  be 
traversed  once.  At  an  odd-order  node,  one  street  must  be  traversed  a second 
time.  Thus,  the  streets  that  must  be  traveled  more  than  once  will  connect 
pairs  of  odd  nodes. 

The  odd-node  pairing  is  performed  In  two  steps:  generation  of  an  initial 
pairing,  and  Iterative  improvement  to  obtain  the  pairing  with  the  minimum  total 
paired  length.  If  collection  Is  required  on  all  streets  in  the  collection  re- 
gion, the  final  odd-node  pairing  will  produce  the  shortest  possible  travel  dis- 
tance in  the  region.  The  presence  of  streets  not  requiring  collection,  however, 
may  cause  the  travel  distance  to  exceed  the  minimum.  Attempts  to  remove  these 
streets  will  be  described  In  the  discussion  of  path  optimization. 

Subroutine  EPXP  is  called  by  PHASE3  to  select  possible  entry  and  exit 
nodes.  The  order  of  each  node  located  where  the  travel  path  enters  and  leaves 
the  collection  region  is  Increased  by  one.  Odd-node  pairing  and  path  genera- 
tion are  performed  for  each  combination  of  entry  and  exit  nodes  because  the 
choice  of  these  nodes  can  change  the  number  of  odd-node  pairs  and  the  total 
travel  distance  in  the  region.  The  path  with  the  shortest  total  distance  will 
be  used.  Subroutine  EPXP  chooses  as  entry  points  the  nodes  closest  to  the 
landfill  or  garage,  provided  that  the  path  from  these  nodes  to  the  landfill  or 
garage  includes  no  other  entry-point  nodes.  The  exit  nodes  are  selected  in  a 
similar  manner,  using  paths  from  the  landfill  to  the  section. 

Since  each  dead-end  or  dangling  segment  has  a node  of  order  one,  a return 
path  to  the  section  from  that  node  is  required.  Rather  than  include  the  gener- 
ation of  these  paths  as  part  of  the  odd-node  pairing,  subroutine  CLOSE  1 is 
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called  by  PHASE3  to  generate'  the  return  paths.  The  subroutine  may  add  addi- 
tional nodes  and  segments  to  the  section  to  provide  a return  path,  or  it  may 
repeat  the  segment,  creating  a U-turn.  The  U-turn  is  used  when  the  shortest 
alternate  path  traverses  more  than  32  segments. 

/ Both  the  initial  odd-node  pairing  and  the  iterative  improvement  require  a 
matrix  of  distances  from  each  odd  node  to  each  other  odd  node.  The  distance 
matrix  is  generated  by  subroutine  GENDM.  This  subroutine  uses  subroutine  TREE 
to  build  paths  starting  at  each  odd  node.  The  distances  to  the  other  odd  nodes 
are  extracted  from  an  array  returned  by  subroutine  TREE. 

Subroutine  SELORD  is  called  by  PHASE3  to  generate  an  initial  pairing  of 
odd  nodes.  Two  pairing  methods  are  coded  in  the  subroutine.  The  fifth  argu- 
ment is  set  to  -1  to  select  one  method,  or  to  +1  to  select  the  other.  The 
first  method  examines  the  distance  from  each  node  to  its  nearest  neighbor. 

The  node  having  the  longest  such  distance  is  paired  with  its  nearest  neighbor. 
If  the  nearest  neighbor  of  any  unpaired  node  has  been  used  as  part  of  an  odd- 
node  pair,  a new  near  neighbor  is  found  for  the  unpaired  node.  The  nodes  with 
the  next  longest  near-neighbor  distance  are  then  paired,  and  the  procedure  is 
repeated  until  all  nodes  are  paired.  The  second  pairing  method  pairs  nodes 
starting  with  the  shortest  near-neighbor  separation.  Again,  the  procedure  is 
repeated  with  the  remaining  unpaired  nodes  until  all  are  paired. 

Since  PHASE3  calls  subroutine  SELORD  with  a fifth  argument  of  -1,  the 
first  method  is  always  used.  The  procedure  is  illustrated  by  the  network  shown 
in  Figure  2(a).  The  six  odd  nodes  are  designated  by  the  letters  A through  E. 
The  lengths  of  the  segments  are  appended  near  the  middle  of  the  segments.  The 
initial  selection  of  near  neighbors  is  shown  in  Figure  2(b).  The  largest  sepa- 
ration has  length  3 and  is  first  encountered  with  nodes  C and  A.  When  nodes 
A and  C are  paired,  node  B loses  its  nearest  neighbor.  A new  near  neighbor, 

F,  is  found  for  node  B.  The  resulting  status  of  the  unpaired  nodes  and  their 
neighbors  is  shown  in  Figure  2(c).  Nodes  B and  F now  have  the  largest  near- 
neighbor separation,  so  they  are  paired.  The  remaining  nodes,  D and  E,  must 
now  be  paired.  The  total  length  of  the  segments  connecting  the  three  pairs 
(A-C,  B-F,  and  D-E)  is  13. 
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Figure  2.  Pairing  of  Odd  Nodes  by  Subroutine  SELORD 


(b)  Node  Neighbors  Before  First  Pairing 


Node  Neighbor 


Separation 


(c)  Node  Neighbors  Before  Second  Pairing 


Node  Neighbor 


Separation 

7 


Notes:  (1)  Letters  indicate  odd  nodes. 

(2)  Numerals  indicate  segment  lengths. 
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To  obtain  pairs  having  a minimum  total  paired  length.  It  Is  necessary  to 
examine  alternating  paths.  An  alternating  path  is  a path  that  alternately 
traverses  segments  connecting  a node  pair  and  then  segments  leading  to  a node 
belonging  to  another  pair.  Using  the  pairing  method  described  above,  A * C 
-*•  D (Figure  2)  is  a two-step  alternating  path.  Similarly,  C-*-A-*B  + F«-E 
is  a four-step  alternating  path.  If  the  alternating  path  starts  and  stops  on 
the  same  node,  it  is  called  an  alternating  cycle. 


The  potential  gain  for  an  alternating  path  is  the  sum  of  the  lengths  of 
the  segments  connecting  paired  nodes  minus  the  sum  of  the  lengths  of  the  seg- 
ments connecting  different  pairs.  All  connecting  lengths  are  taken  along  the 
shortest  path.  For  example,  A > C -*•  0 has  a potential  gain  of  -2,  while  C * 

A -*■  B -*•  F -*  E has  a potential  gain  of  6.  The  alternating  cycle  C •*  A -►  B ♦ F 
-*•  E 0 ■*  C has  the  following  potential  gain:  3-2+7-2+3-5=4.  A 
cycle  with  positive  potential  gain  is  called  a profitable  cycle.  When  a pair- 
ing is  found  that  has  no  profitable  cycles,  that  pairing  has  the  minimum  total 
pairing  length. 

The  improvement  of  the  odd-node  pairing  starts  with  the  formation  of  a 
matrix  of  potential  gains  for  two-step  alternating  paths.  Let  PGp(I,J)  be  the 
largest  potential  gain,  or  profit,  of  any  path  from  1 to  J found  while  comput- 
ing the  N-step  potential  gain  matrix.  Paths  that  end  at  the  partner  of  the 
starting  node  are  assigned  large  negative  potential  gains  (-1000).  Two-step 
paths  that  start  and  end  on  the  same  node  are  assigned  potential  gains  of  0. 
The  shortest  pair-connecting  path  may  traverse  the  segments  connecting  a pair, 
but  is  still  treated  as  an  unpaired  length.  For  example,  A to  E is  A •+  C 
(paired)  -*•  E via  D (pairs  ignored)  for  a potential  gain  (PG)  of  3 - (5  + 3) 

= -5.  The  two-step  alternating  path  from  C to  D is  C -*■  A (paired)  -*■  D via  C 
(unpaired).  The  PG , matrix  for  Figure  2 is  shown  in  Table  1. 


The  PG,  matrix  is  formed  in  subroutine  SOLV.  Subroutine  NEXTM  is  called 
by  subroutine  SOLV  to  double  the  number  of  steps  in  the  PG  matrix.  The  old 
matrix  is  not  destroyed  by  the  new  matrix.  Subroutine  NEXTM  computes  the  di- 
agonal entries  (cycles)  first.  The  entry  from  node  I to  node  1 is  evaluated 
using  the  approach  shown  in  the  following  equation. 
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TABLE  1.  TWO-STEP  PROFIT  MATRIX 


STOP 


N 

A 

B 

C 

D 

E 

F 

A 

0 

-1 

-1000 

-2 

-5 

-6 

B 

-2 

0 

-2 

3 

5 

-1000 

C 

-1000 

1 

0 

-5 

-4 

-6 

D 

-4 

-6 

-5 

0 

-1000 

1 

E 

-5 

-6 

-2 

-1000 

0 

-1 

F 

5 

-1000 

3 

-2 

-2 

0 

PG4(I , I ) = max 


pg2(i.j) 


pg2(j 


.1)  , 


pg2(i,i) 


If  the  cycle  is  profitable  (PG  > 0),  it  is  checked  for  repeated  nodes.  If  any 
nodes  other  than  the  ends  are  repeated,  the  cycle  is  not  allowed.  If  no  valid 
profitable  cycle  exists,  the  remainder  of  the  matrix  is  formed.  The  four-step 
potential  gain  is  computed  as  follows: 

max|pG2(I,K)  + PG2(K,J)j  , PG2(I,J) 

m 
m 

The  PG4  matrix  for  Figure  2 is  shown  in  Table  2. 

Since  no  profitable  cycles  were  found  in  the  PG4  matrix,  subroutine  NEXTM 
is  called  again  to  form  the  PGa  matrix.  In  forming  PG  cycles,  a valid  cycle 
A"C-*-D-*E-*-F-*B-*A  is  found  with  profit  +4.  At  this  point,  the  node 
pairing  in  the  profitable  cycle  is  reversed.  The  new  pairs  are  C-D,  E-F,  and 
A-B.  Subroutine  SOLV  recomputes  the  PG„  matrix.  Subroutine  NEXTM  is  called 
to  compute  the  PG4  matrix.  When  the  PGR  matrix  is  computed,  no  profitable 
cycles  are  found.  Since  the  eight-step  path  Includes  all  the  odd-node  pairs, 
the  present  pairing  has  a minimum  total  pairing  length  of  9. 

] 


When  the  optimum  odd-node  pairing  has  been  obtained,  PHASE3  calls  subrou- 
tine TREE  to  find  the  shortest  path  within  the  section  that  connects  each 


PG4(I ,J)  = max 
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TABLE  2.  FOUR-STEP  PROFIT  MATRIX 


STOP 


A 

B 

c 

D 

E 

F 

A 

0 

-1 

-1000 

2 

4 

-1 

B 

0 

0 

3 

3 

5 

-1000 

C 

-1000 

1 

0 

4 

6 

-4 

D 

6 

-4 

4 

0 

-1000 

1 

E 

4 

-1 

2 

-1000 

0 

-1 

F 

5 

-1000 

3 

3 

0 

0 

odd-node  pair.  The  count  of  times  each  segment  is  to  be  traversed  is  incre- 
mented by  one. 

Main  program  PHASE3  calls  subroutine  TRAVEL  to  generate  a path  from  the 
landfill  or  garage  through  the  collection  region  and  back  to  the  landfill.  The 
trial -and-error  search  used  by  subroutine  TRAVEL  is  adequate  to  find  a path 
because  there  are  usually  many  paths  having  the  minimum  length.  For  example, 
there  are  more  than  140  minimum-length  paths  starting  and  ending  at  node  A in 
Figure  2. 

Segments  that  are  traversed  more  than  once  are  placed  at  the  beginning  of 
the  segment  arrays  to  minimize  the  number  of  U-turns.  Since  the  segments  are 
examined  sequentially  from  first  to  last  for  addition  to  the  path,  the  segments 
requiring  more  than  one  traversal  will  be  traveled  the  first  time  a node  bound- 
ing one  of  them  is  reached.  Unless  the  segment  is  a dead-end,  use  of  this 
method  usually  precludes  U-turns. 

When  no  segments  can  be  added  to  the  end  of  the  path,  subroutine  TRAVEL 
examines  the  segment  data  for  segments  that  must  still  be  traversed.  If  no 
unused  segments  remain,  a suitable  travel  path  has  been  found.  Otherwise, 
the  path  is  retraced  in  reverse  until  an  unused  segment  is  reached.  At  that 
point,  the  trial-and-error  search  resumes. 

After  subroutine  TRAVEL  has  found  a path  for  traversing  the  collection 
region,  PHASE3  calls  subroutine  OPTPATH  to  perform  two  types  of  optimization 
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on  the  path.  The  subroutine  finds  the  shortest  distance  between  the  points 
that  start  and  end  each  travel  stretch.  Because  segments  that  are  not  in  the 
section  may  be  used,  the  travel,  stretch  may  be  shortened.  The  second  type  of 
optimization  removes  twice- traveled  stretches  on  which  collection  is  not  re- 
quired. If  a part  of  the  path  is  traveled  twice  in  the  same  direction,  and  if 
. there  are  no  one-way  streets  on  the  path  between  the  end  of  the  first  traversal 

and  the  beginning  of  the  second,  the  twice-traveled  piece  is  removed  both  times 
and  the  path  between  the  first  and  second  traversals  is  reversed. 

As  travel  paths  are  generated  for  each  combination  of  entry  and  exit 
nodes,  the  total  length  of  each  new  travel  path  is  compared  to  that  of  the 
best  found  thus  far.  When  all  combinations  of  entry  and  exit  nodes  have  been 
tried,  PHASE3  writes  the  best  path  to  file  TAPE9.  At  the  conclusion  of  pro- 
gram PHASE3,  TAPE9  contains  two  choices  of  travel  for  servicing  each  s’ction, 
one  choice  starting  at  the  garage  and  one  starting  at  the  landfill. 

2.  DATA  STORAGE 

Program  PHASE3  obtains  data  from  five  sources:  card  input  and  files 
TAPE1 , TAPE2,  TAPE3,  and  TAPE4.  Files  TAPE1  and  TAPE2  hold  the  segment  and 
node  data  produced  by  program  RCINPT.  Files  TAPE3  and  TAPE4  hold  the  vehicle 
data  and  section  list  produced  by  program  PHASE2.  One  file,  TAPE9,  is  gener- 
ated by  program  PHASE3  and  is  saved  on  disk  for  use  by  program  PHASE4. 

The  problem  title  is  read  from  the  first  data  card  at  the  beginning  of 
main  program  PHASE3.  The  node  numbers  of  the  landfill  and  garage  are  read 
from  the  second  data  card  following  statement  36.  A segment  number  and  new 
section  number  are  read  from  each  remaining  data  card  at  statement  50  of 
• PHASE3.  All  of  the  variables  read  from  data  cards  are  kept  in  storage  local 

to  PHASE3. 

All  of  the  segment  data  on  file  TAPE1  are  read  by  a single  READ  statement 
shortly  after  statement  16  of  PHASE3.  The  count  of  segments,  NSEG,  and  the 
segment  data  array,  STG,  are  stored  in  blank  COMMON.  Variable  AVMD,  the  map 
distance  conversion  factor,  is  not  used  by  the  program.  The  parts  of  the  STG 
array  with  second  subscripts  equal  to  9 and  10  are  not  read  from  TAPE! . 
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All  of  the  node  data  are  read  from  file  TAPE2  by  a single  READ  statement 
shortly  after  statement  18  of  program  PHASE3.  The  total  house  count  and  the 
total  refuse  quantity,  NHTOT  and  TOTREF,  are  kept  in  local  storage.  The  re- 
maining variables  are  stored  in  COMMON  block  NODDATA.  The  data  in  this  COMMON 
block  are  not  changed  during  the  remainder  of  the  program. 

All  of  the  data  on  file  TAPE3  are  read  by  a single  READ  statement.  This 
statement  is  the  second  executable  statement  in  PHASE3.  The  segment  count, 
NSEG,  is  stored  in  blank  COMMON,  but  will  be  overwritten  by  a subsequent  READ 
statement.  The  number  of  sections,  NTRIP,  is  kept  in  storage  local  to  PHASE3. 
Arrays  SINL  and  SFNL  are  stored  in  COMMON  block  TEMSTG.  The  pointers  in  these 
arrays  are  used  only  in  the  loop  through  statement  16  in  the  main  program.  The 
arrays  will  be  overwritten  by  subroutine  CONNECT.  The  array  of  vehicle  capa- 
cities, TC,  is  kept  in  storage  local  to  PHASE3.  These  values  are  not  changed. 

The  segment  numbers  on  file  TAPE4  are  ordered  by  section.  The  segment 
numbers  are  read,  one  at  a time,  by  the  READ  statement  in  the  loop  through 
statement  14.  The  segment  numbers  are  stored  in  the  ISTG  array,  with  the 
second  subscript  increasing  by  1 for  each  segment.  The  section  number  is  also 
stored  in  the  ISTG  array  as  each  segment  is  read.  Following  the  loop  through 
statement  16,  the  segment  numbers  are  sorted  into  increasing  order;  the  sec- 
tion assignments  are  carried  along.  The  segment  numbers  will  be  overwritten 
when  the  remaining  segment  data  are  read  from  file  TAPE1. 

Data  are  written  to  file  TAPE9  from  two  places  in  main  program  PHASE3. 

Two  WRITE  statements  shortly  before  statement  86  write  the  path  from  the  land- 
fill to  the  garage  to  TAPE9.  The  remaining  path  information  is  written  by  six 
WRITE  statements  immediately  preceding  statement  500.  The  statements  write 
header  cards  and  path  descriptions  for  the  three  parts  of  the  travel  path: 
travel  from  the  garage  or  landfill  to  the  collection  region,  travel  and  col- 
lection within  the  collection  region,  and  travel  from  the  collection  region 
to  the  landfill.  These  six  WRITE  statements  are  executed  twice  for  each  sec- 
tion, first  for  the  trip  starting  at  the  garage,  and  second  for  the  trip  start- 
ing at  the  landfill,  as  part  of  the  loop  through  statement  500.  The  loop 
through  statement  800  repeats  the  loop  through  statement  500  for  each  section. 


The  data  for  the  path  from  the  landfill  back  to  the  garage  come  from 
arrays  generated  by  the  calls  to  subroutine  TREE  shortly  after  statement  80 
of  the  main  program.  The  first  call  to  TREE  generates  path  information  for 
paths  from  the  landfill  node  to  each  other  node.  The  subsequent  call  to 
TRACE  puts  the  path  information  from  the  garage  to  the  landfill  in  arrays 
IPS  and  IPN.  Subroutine  FLIP  reverses  the  order  of  this  path.  The  IPN  and 
IPS  arrays  are  written  to  TAPE9.  The  arrays  are  used  for  other  path  infor- 
mation later  in  the  main  program.  Travel -path  data  to  and  from  the  landfill 
are  obtained  in  a similar  manner,  using  subroutines  TRACE  and  FLIP  near  state- 
ment 410  in  the  main  program.  Arrays  ISTS  and  NDTS  are  used  for  travel  to  the 
section,  and  arrays  ISFS  and  NDFS  are  used  for  travel  from  the  section.  The 
data  are  transferred  to  arrays  ISSV  and  NDSV  in  the  loops  through  statements 
430  and  450.  The  data  for  the  path  in  the  collection  region  originate  in  sub- 
routine TRAVEL.  The  path  is  stored  in  the  IPN  and  IPS  arrays.  It  may  be 
changed  by  subroutine  OPTPATH.  It  is  transferred  to  arrays  ISSV  and  i\DS V in 
the  loop  through  statement  440.  All  of  the  arrays  used  to  hold  path  infor- 
mation are  In  COMMON  block  TEMSTG.  These  data  will  be  overwritten  by  the  next 
call  to  subroutine  SOLV. 

The  tree  data  structure  in  subroutine  TREE  and  the  packed  path  matrices 
used  in  odd-node  pairing  are  complicated  enough  to  require  individual  descrip- 
tions. Double-subscripted  array  1ST  is  used  in  subroutine  TREE  to  hold  a bi- 
nary tree  representation  of  paths  to  or  from  a given  starting  node.  Each  TREE 
entry  is  made  up  of  four  words  in  array  1ST  having  the  same  first  subscript: 
the  line  number  of  the  node,  a left  pointer,  a right  pointer,  and  a back 
pointer.  The  pointers  contain  the  line  number  (the  first  subscript)  of  other 
entries  in  the  1ST  array.  The  root  of  the  tree  is  the  starting  node.  Right 
branches  lead  to  nodes  reached  on  the  next  step,  and  left  branches  connect  all 
nodes  reached  from  one  node  by  an  additional  step.  Each  left  pointer  of  the 
leftmost  branches  points  to  a node  one  step  closer  to  the  starting  node.  The 
pointer  contains  the  negative  of  this  node's  line  number  in  array  1ST.  The 
left  and  right  pointers  are  used  to  scan  all  entries  in  the  tree.  The  back 
pointers  are  used  to  facilitate  unlinking  entries  as  the  tree  is  being  built. 

The  stages  of  formation  of  a tree  for  the  network  presented  In  Figure  2 
are  shown  In  Figure  3.  The  four  small  triangles  forming  the  larger  triangle 
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Figure  3.  Stages  of  Tree  Formation 
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shown  in  Figure  3(a)  represent  the  four  elements  in  each  tree  entry.  The  cen- 
ter triangle,  corresponding  to  the  line  number  of  the  node,  contains  the  letter 
designating  the  node.  The  bottom  left  and  right  triangles  represent  the  left 
and  right  pointers.  The  top  triangle  represents  the  back  pointer.  The  network- 
connections  represented  at  various  stages  of  tree  formation  are  shown  to  the 
right  of  the  trees.  Node  A is  the  starting  node. 

The  entry  for  node  A is  set  up  in  subroutine  TREE  preceding  the  level- 
building loop.  Three  nodes,  B,  C,  and  E,  are  added  to  the  tree  by  the  first 
pass  through  the  level -bui lding  loop.  The  order  in  the  tree  depends  on  the 
order  of  the  data,  rather  than  on  lengths  between  nodes.  One  possible  order 
is  shown  in  Figure  3(b).  Nodes  D and  F are  added  to  the  tree  as  second-level 
nodes  to  node  E,  as  shown  in  Figure  3(c).  When  second-level  nodes  are  added 
to  node  C,  as  shown  in  Figure  3(d),  node  D is  disconnected  from  node  E and 
appended  to  node  C because  its  distance  from  node  A is  shorter  via  node  C than 
via  node  E.  The  final  tree  and  all  pointer  connections  are  shown  in  Figure  A. 

The  packed  path  matrix  contains  only  those  portions  of  the  paths  which 
cannot  be  inferred.  Because  the  paths  connect  pairs  of  odd  nodes,  only  one 
node  in  the  pair  need  be  specified;  the  other  is  obtained  from  the  partner 
table.  The  starting  and  ending  nodes  can  be  inferred  from  the  row  and  column 
in  the  matrix;  therefore,  they  need  not  be  specified.  For  example,  in  the 
cycle  A-^C^D+E+F^B-A,  the  information  needed  to  specify  the  path 
is  D,F.  The  sequence  numbers  of  these  nodes  would  be  packed  into  one  word 
in  the  path  matrix.  The  packing  is  accomplished  by  multiplying  the  current 
path  information  by  one  more  than  the  number  of  nodes  and  adding  in  the  new 
node's  sequence  number.  In  the  current  example,  assume  node  D has  sequence 

number  4 and  node  F has  sequence  number  6.  There  are  6 nodes,  so  packing  re- 

quires multiplication  by  7.  The  path  word  would  be  set  equal  to  4 to  indicate 
the  presence  of  node  D.  When  aode  F is  packed  into  the  word,  the  4 is  multi- 
plied by  7,  and  6 is  added.  The  number  of  nodes  packed  into  a word  is  limited 

so  that  48  or  fewer  bits  are  used.  A second  word  is  available  in  a second 
matrix  if  needed. 
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Figure  4.  Final  Tree  for  Sample  Network 


3.  PURPOSE  AND  PERFORMANCE 

In  this  section  the  simpler  subroutines  are  described  first  so  their  work 
ings  will  be  understood  when  they  are  mentioned  again  in  the  descriptions  of 
the  more  complicated  subroutines.  The  main  program  is  described  last.  Logic 
flowcharts  are  given  in  Appendix  A.  Complete  program  listings  are  provided 
in  Appendix  B.  In  Appendix  C,  the  more  important  FORTRAN  variables  are  de- 
scribed for  each  subroutine. 

a.  Subroutine  FLIP 


Subroutine  FLIP  reverses  the  order  of  an  array.  The  first  argument 
is  the  array,  and  the  second  argument  is  the  number  of  items  to  be  reversed. 

Variable  ND2  is  set  to  half  the  number  of  items  In  the  array.  Vari- 
able NP1  Is  set  to  one  more  than  the  number  of  items  in  the  array.  The  loop 
through  statement  10  reverses  the  X-array.  An  item  X(I)  from  the  front  of  the 
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array  Is  saved  In  variable  XT.  An  Item  X(NP1  - I)  from  the  end  of  the  array 
Is  stored  at  X(I).  At  statement  10  the  value  saved  as  XT  is  stored  at  X (NP1 
- I).  The  loop  Is  executed  with  the  Index,  I,  running  from  1 through  ND2. 
Control  then  returns  to  the  calling  program. 

b.  Subroutine  M0VE3 

Subroutine  M0VE3  moves  all  data  starting  at  and  following  a given 
subscript  to  a position  in  the  table  starting  at  a given  higher  subscript. 
Argument  II  is  the  subscript  from  which  the  data  will  be  moved,  and  IF  Is  the 
subscript  to  which  the  data  will  be  moved.  Arguments  A1 , A2,  and  A3  are  the 
arrays  In  which  data  will  be  moved.  Subroutine  M0VE3  starts  by  determining 
whether  the  data  are  to  be  moved  to  a higher  subscript.  If  the  final  sub- 
script, IF,  is  less  than  the  Initial  subscript,  II,  control  returns  to  the 
calling  program.  If  the  initial  and  final  subscripts  are  in  the  proper  rela- 
tion, the  subroutine  executes  the  loop  through  statement  10  to  move  each  item 
starting  at  subscript  II  the  distance  necessary  to  reach  IF.  All  data  between 
subscripts  II  and  IF,  starting  with  the  last  item,  are  moved.  Motion  starts 
at  the  higher  subscript  so  that  no  data  that  must  be  moved  later  are  over- 
written. Zeros  are  stored  at  location  II  of  the  arrays.  Control  returns  to 
the  calling  program. 

c.  Function  IFIND 

Function  IFIND  uses  a binary  search  to  find  a given  number  in  an 
array  and  assigns  the  subscript  of  the  number  as  the  value  of  IFIND.  If  the 
number  Is  not  found,  the  function  sets  the  value  of  IFIND  equal  to  the  nega- 
tive of  the  subscript  at  which  the  number,  to  be  in  numerical  order,  should  be 
Inserted.  (The  array  Is  assumed  to  be  In  Increasing  order.)  There  are  three 
arguments.  Argument  NUM  is  the  number  that  is  sought  in  array  1ARRAY.  The 

second  argument,  IARRAY,  is  the  array  to  be  searched.  Argument  LEN  is  the 

length  of  array  IARRAY. 

Function  IFIND  begins  by  checking  that  LEN  > 0.  If  LEN  < 0,  the 
function  assigns  a value  of  -1  to  IFIND  and  control  returns  to  the  calling  pro- 
gram. This  value  Indicates  that  the  number  sought  is  not  in  the  array  and  would 
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be  stored  as  the  first  entry  In  the  array.  The  binary  search  uses  variables 
II,  IP,  and  IF  as  pointers.  II  is  the  subscript  of  the  front  of  the  region 
being  searched,  IP  is  the  subscript  of  the  item  being  compared  to  the  number 
sought,  and  IF  is  the  subscript  of  the  last  item  in  the  region  being  searched. 
Variable  II  is  set  to  1 at  statement  5.  Variable  IF  is  set  to  the  length  of 
the  array.  The  pointer,  IP,  is  the  subscript  midway  between  II  and  IF  and  is 
computed  at  statement  10. 

Following  statement  10,  NUM  is  compared  to  lARRAY(IP).  If  NUM  < 
lARRAY(IP),  control  transfers  to  statement  20,  indicating  that  the  number  is 
in  the  front  half  of  the  region  being  searched.  At  statement  20  the  final 
pointer  is  moved  to  the  subscript  preceding  the  point  just  searched.  If  NUM 
> IARRAY(IP),  control  transfers  to  statement  30,  indicating  that  the  number 
being  sought  follows  the  subscript  just  inspected.  At  statement  30  tha  initial 
pointer,  II,  is  set  to  the  present  pointer,  IP,  plus  1.  If  the  numbe.  sought 
is  found  at  IARRAY (IP),  control  transfers  to  statement  50,  where  IFIND  is  set 
equal  to  the  current  pointer.  Control  returns  to  the  calling  program.  Where 
NUM  is  unequal  to  IARRAY(IP),  the  initial  or  final  pointer  is  moved  and  con- 
trol resumes  at  statement  40.  At  statement  40  the  final  pointer  is  compared 
to  the  initial  pointer.  If  IF  > II,  control  transfers  to  statement  10,  where 
the  search  resumes  on  the  appropriate  half  of  the  region  examined  previously. 

If  the  final  pointer  becomes  less  than  the  initial  pointer,  the  number  sought 
is  not  in  the  table.  In  this  case,  control  resumes  following  statement  40, 
and  the  value  of  IFIND  is  set  to  the  negative  of  the  current  pointer.  If  the 
number  at  the  current  pointer  is  less  than  the  number  being  sought,  IFIND  is 
set  to  -(IP  + 1)  so  the  number  can  be  inserted  in  the  appropriate  place.  Con- 
trol then  returns  to  the  calling  program. 

d.  Subroutine  SHLSRT3 

Subroutine  SHLSRT3  sorts  one  array  and  carries  two  other  arrays  along 
during  the  sorting.  The  subroutine  uses  Shell's  sorting  algorithm. 

Subroutine  SHLSRT3  has  five  arguments.  The  first,  X,  is  the  array 
to  be  sorted.  The  second  and  third  arguments  are  arrays  that  are  paired  with 
X and  are  rearranged  as  X is  sorted.  The  fourth  argument  is  the  number  of 
words  to  be  sorted.  The  fifth  argument  is  set  to  +1.0  if  the  array  is  to  be 
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sorted  Into  Increasing  order  or  to  -1.0  If  It  is  to  be  sorted  into  decreasing 
order.  The  numbers  are  sorted  by  a procedure  in  which  pairs  of  numbers  are 
compared  and  Interchanged,  if  necessary,  to  put  the  smaller  number  closer  to 
the  beginning  of  the  array.  The  separation  of  the  numbers  compared  is  ap- 
proximately half  the  number  of  entries  in  the  array.  This  spacing  is  halved 
in  subsequent  passes  through  the  array.  When  two  numbers  are  interchanged, 
the  pointers  are  moved  up  so  that  the  smaller  number  is  compared  to  a number 
farther  up  in  the  array. 

The  spacing,  N,  is  set  initially  to  half  the  number  of  words.  The 
number  of  comparisons,  K,  to  be  performed  In  the  loop  through  statement  50  is 
the  total  number  of  words  less  N.  The  loop  through  statement  50  uses  Index  I 
as  one  of  the  pointers.  This  pointer  is  stored  in  variable  J;  the  other 
pointer,  L,  Is  set  to  I+.N.  The  values  of  the  array  to  be  sorted  and  the  arrays 
to  be  carried  along  are  saved  as  XT,  AT,  and  BT.  The  values  of  X at  the  loca- 
tions indicated  by  the  pointers  are  compared.  If  they  are  in  order,  control 
transfers  to  statement  40.  If  not,  the  larger  value  is  stored  closer  to  the 
end  of  the  array.  The  pointers  are  both  moved  up  by  N.  If  the  pointer  having 
the  smaller  value  is  a valid  subscript,  control  transfers  to  statement  20, 
where  another  comparison  is  performed.  At  statement  40  the  saved  values  are 
stored  in  the  appropriate  places  in  the  arrays.  Following  the  loop  through 
statement  50,  if  the  spacing  is  equal  to  1,  the  sort  is  complete  and  control 
returns  to  the  calling  program.  Otherwise,  the  spacing  is  halved  and  control 
transfers  to  statement  10. 

e.  Subroutine  ISHLSRT 

Subroutine  ISHLSRT  sorts  one  array  and  carries  along  another  array 
during  the  sorting.  There  are  five  arguments.  The  first  argument,  IX,  is  an 
array  of  integers  to  be  sorted.  The  second  argument,  IA,  is  an  array  that  is 
reordered  as  IX  Is  sorted.  The  third  argument,  NW,  is  the  number  of  words  to 
be  sorted.  The  fourth  argument,  ISGN,  is  set  to  1 if  the  array  is  to  be  sorted 
into  increasing  order  or  to  -1  if  the  array  is  to  be  sorted  into  decreasing 
order.  The  fifth  array,  INC,  is  the  spacing  of  the  subscripts  on  IX  and  IA. 

The  items  that  will  be  sorted  and  reordered  fall  at  subscripts  1,  INC  + 1, 

2 x INC  ♦ 1,  and  so  forth. 
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The  Internal  workings  of  subroutine  ISHLSRT  are  the  same  as  those 
of  subroutine  SHLSRT3  but  for  three  minor  differences:  the  array  to  be  sorted 
Is  Integer  Instead  of  floating  point,  only  one  array  is  carried  along  with  the 
sorted  array,  and  an  increment  other  than  1 is  allowed  on  the  subscripts.  The 
increment  is  implemented  by  treating  arrays  IX  and  IA  as  double-dimensioned 
arrays  with  a first  dimension  equal  to  INC.  All  references  to  these  arrays 
use  a first  subscript  equal  to  1 and  a second  subscript  which  is  the  same  as 
that  used  in  subroutine  SHLSRT3. 


f.  Subroutine  ADJUST 

Subroutine  ADJUST  puts  segments  meeting  at  order-two  nodes  into  the 
same  section  if  vehicle  capacity  will  not  be  exceeded.  It  has  two  arguments. 
The  first  argument,  TL,  is  an  array  of  vehicle  loads.  The  second  argi  tent, 

TC,  is  an  array  of  vehicle  capacities. 

The  loop  through  statement  70  examines  each  of  the  KNODES  nodes.  If 
the  count  of  segments  bounding  a node,  NUMNBR,  is  not  equal  to  2,  control 
transfers  to  statement  70.  Otherwise,  the  numbers  of  the  segments  bounding 
the  node  are  obtained  from  the  NBRSEG  array  and  are  stored  in  variables  NSGl 
and  NSG2.  The  section  assignments  of  these  segments  are  obtained  from  the 
ISTG  array  and  are  stored  in  variables  NSC1  and  NSC2.  If  the  segments  are  in 
the  same  section,  or  if  either  is  in  section  zero,  control  transfers  to  state- 
ment 70.  Otherwise,  refuse  quantities  R1  and  R2  are  computed  for  the  segments. 
If  both  R1  and  R2  are  zero,  control  transfers  to  statement  70.  Otherwise,  R1 
is  compared  with  R2.  If  R2  < R1 , control  transfers  to  statement  30.  If  not, 
and  if  the  sum  of  R1  and  the  load  of  the  vehicle  servicing  section  NSC2  exceeds 
the  vehicle  capacity,  control  transfers  to  statement  50.  Otherwise,  at  state- 
ment 10  the  load  for  the  vehicle  servicing  section  NSC2  is  Incremented  by  Rl, 
and  the  load  for  the  vehicle  servicing  NSC1  is  decremented  by  Rl . The  section 
assignment  for  segment  NSGl  is  changed  to  NSC2.  A message  is  printed  Indica- 
ting that  the  program  has  adjusted  the  section  assignment.  Control  transfers 
to  statement  70. 

Statement  30  is  reached  when  the  refuse  quantity,  R2,  on  the  second 
segment  is  less  than  or  equal  to  the  refuse  quantity  on  the  first  segment.  If 
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the  sum  of  R2  and  the  load  of  the  vehicle  servicing  section  NSC1  is  greater 
than  the  vehicle  capacity,  control  transfers  to  statement  60.  Otherwise,  at 
statement  40,  R2  Is  added  to  the  load  of  the  vehicle  servicing  secton  NSC1  and 
is  subtracted  from  the  load  of  the  vehicle  servicing  section  NSC2.  The  section 
assignment  of  segment  NSG2  is  changed  to  NSC1,  and  a printed  message  indicates 
the  change.  Control  transfers  to  statement  70.  At  statement  50,  if  the  sum  of 
R2  and  the  load  of  the  vehicle  servicing  section  NSC1  exceeds  the  vehicle  capa- 
city, control  transfers  to  statement  70.  If  not,  control  transfers  to  statement 
40.  At  statement  60,  if  the  sum  of  R1  and  the  load  of  the  vehicle  servicing 
section  NSC2  is  less  than  or  equal  to  the  vehicle  capacity,  control  transfers 
to  statement  10.  If  not,  statement  70  is  reached,  and  the  loop  is  repeated  for 
the  next  node.  After  the  last  node  has  been  processed,  control  returns  to  the 
calling  program. 

g.  Function  CUMDIS 

Function  CUMDIS  evaluates  the  total  length  of  a group  of  segments. 

It  has  two  arguments:  ISEG,  an  array  of  segment  numbers,  and  NSG,  a count  of 
the  segments. 

Variable  SUM  is  initially  set  to  0.  The  loop  through  statement  10 
accumulates  the  lengths  of  the  segments  in  variable  SUM.  If  the  segment  number 
is  greater  than  NSEG,  the  number  of  segments  in  the  map  description,  the  segment 
does  not  have  a street  number,  but  Its  segment  number  is  obtained  from  the  stor- 
age that  normally  would  contain  the  street  number.  Following  the  loop,  CUMDIS 
is  assigned  a value  equal  to  SUM  and  control  returns  to  the  calling  program. 

h.  Subroutine  PRNPCS 

Subroutine  PRNPCS  prints  the  segment  numbers  of  the  segments  in  each 
piece  of  a section.  It  has  two  arguments.  The  first  argument,  NTRIP,  gives 
the  trip  or  section  number.  The  second  argument,  NPIECE,  gives  the  maximum 
number  of  pieces  of  the  section  to  be  printed. 

Variable  IZER0  is  set  to  0.  The  loop  through  statement  40  examines 
each  of  NPIECE  pieces  of  the  section.  Variable  ITRIP  is  set  equal  to  the  trip 


number  plus  an  appropriate  multiple  of  1000.  This  number  is  used  to  identify 
the  piece  of  the  trip.  A carriage  control  variable,  CC,  is  set  equal  to  a 
blank.  Variable  NN,  a count  of  the  segment  numbers  on  the  current  line  of  out- 
put, is  set  to  0.  A heading  giving  the  piece  number  is  printed. 

\ 

The  loop  through  statement  30  scans  each  line  of  the  segment  table. 

If  the  piece  number  for  the  segment  is  not  I TR IP,  control  transfers  to  state- 
ment 30.  Otherwise,  NN  is  incremented  by  1.  The  segment  number  is  printed. 
Carriage  control  variable  CC  is  changed  to  a plus  sign,  causing  subsequent 
segment  numbers  to  be  printed  on  the  same  line  as  the  first.  If  fewer  than 
30  segments  have  been  printed  on  the  line,  control  transfers  to  statement  30. 
Otherwise,  CC  is  set  to  a blank  and  NN  is  reset  to  0.  After  the  loop  through 
statement  40  has  been  completed,  control  returns  to  the  calling  program. 

i.  Subroutine  TRACE 

Subroutine  TRACE  finds  the  segments  and  nodes  in  a path  between  two 
nodes.  It  has  seven  arguments.  The  first  two  arguments,  NSTART  and  NSTOP, 
are  the  node  numbers  of  the  starting  and  stopping  points  in  the  path.  The 
third  argument,  LPREVN,  is  an  array  containing  the  line  number  of  the  next 
node  in  the  path  to  node  NSTOP.  The  fourth  argument,  NPREVS,  is  an  array  con- 
taining the  number  of  the  next  segment  in  the  path  to  node  NSTOP.  The  fifth 
argument,  array  IPS,  will  be  filled  with  the  segment  numbers  in  the  path  by 
subroutine  TRACE.  The  sixth  argument,  array  IPN,  will  be  filled  with  the  node 
numbers  in  the  path.  The  seventh  argument,  NSG,  will  be  set  to  a count  of  the 
segments  in  the  path  from  node  NSTART  to  node  NSTOP. 

Function  IFIND  finds  the  line  number  of  node  NSTART  in  array  NODNUM. 
The  line  number  is  stored  in  variable  L.  The  starting  node  number  is  stored 
in  IPN ( 1 ) . The  line  number  of  the  terminal  node,  NSTOP,  is  found  by  function 
IFIND  and  is  stored  in  LEND.  The  segment  count,  NSG,  is  set  to  0. 

At  statement  10  the  segment  count  is  incremented  by  1.  The  segment 
number  of  the  next  segment  in  the  path  is  obtained  from  the  NPREVS  array  and 
is  stored  in  the  IPS  array.  The  line  number  of  the  next  node  in  the  path  is 
obtained  from  the  LPREVN  array  and  is  stored  in  variable  L.  The  number  of 
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this  node  is  stored  In  the  IPN  array.  If  the  line  number  of  this  node  Is  nega- 
tive or  is  greater  than  the  count  of  nodes,  or  if  the  count  of  segments  in  the 
path  exceeds  200,  control  transfers  to  an  error  message  at  statement  20.  Other- 
wise, if  the  line  number,  l,  is  not  equal  to  the  line  number  of  node  NSTOP,  and 
if  L is  nonzero,  control  returns  to  statement  10.  Otherwise,  control  returns 
to  the  calling  program.  At  statement  20,  an  error  message  is  printed,  giving 
the  information  about  the  path  found  thus  far  between  nodes  NSTART  and  NSTOP. 
System  subroutine  DUMP  is  called  to  give  a core  dump  and  terminate  the  job. 

j.  Subroutine  MOVODD 

Subroutine  MOVODD  moves  odd-order  node  pointers  to  the  front  of  an 
array  and  returns  a count  of  the  odd  nodes.  It  has  four  arguments.  The  first 
argument,  KN,  is  a count  of  all  the  nodes.  The  second  argument,  NODORD,  is  an 
array  of  node  orders.  The  third  argument,  NODPTR,  is  an  array  of  node  pointers. 
The  fourth  argument,  NODD,  is  a count  of  odd  nodes  and  is  returned  by  subrou- 
tine MOVODD. 

Variable  JJ  is  set  to  0.  Variable  KK  is  set  to  one  more  than  the 
number  of  nodes.  These  variables  will  be  pointers  at  opposite  ends  of  the 
NODORD  array.  At  statement  10  the  front  pointer,  JJ,  is  incremented  by  1 . If 
JJ  is  greater  than  or  equal  to  the  final  pointer,  KK,  control  transfers  to 
statement  30.  Otherwise,  NODORD(JJ)  is  examined.  If  it  is  odd,  control  trans- 
fers to  statement  10.  Otherwise,  at  statement  20,  variable  KK  is  decremented 
by  1.  If  JJ  > KK,  control  transfers  to  statement  30.  If  NODORD (KK)  is  even, 
control  returns  to  statement  20.  If  not,  the  orders  and  pointers  to  the  odd 
and  even  nodes  are  interchanged  so  that  the  odd-node  data  are  at  the  front  end 
of  the  arrays.  Control  returns  to  statement  10.  Statement  30  is  reached  when 
pointers  JJ  and  KK  are  equal.  At  statement  30  the  count  of  odd  nodes  is  set 
equal  to  JJ  - 1.  Control  returns  to  the  calling  program. 

k.  Subroutine  TREE 

Subroutine  TREE  builds  four  arrays  for  each  node  that  can  be  reached 


from  a given  starting  node.  The  arrays  give  the  line  number  of  the  previous 
node,  the  number  of  the  previous  segment,  and  the  travel  time  and  distance  be- 
tween the  startina  node  and  the  node  in  Question.  The  subroutine  hac  nine 


arguments.  The  first  argument,  NORG,  is  the  number  of  the  starting  node.  The 
second  argument,  DIST,  is  an  array  of  distances  from  the  starting  node  to  each 
other  node.  The  third  argument,  TIME,  is  an  array  of  travel  times  from  the 
starting  node  to  each  other  node.  The  fourth  argument,  W,  indicates  whether 
time  or  distance  or  a linear  combination  of  the  two  is  minimized  as  paths  from 
the  starting  node  are  built.  The  fifth  argument,  NPREVN,  is  an  array  that  gives 
the  line  number  of  the  previous  node  in  the  path  for  each  node.  The  sixth  argu- 
ment, NPREVS,  is  an  array  that  gives  the  number  of  the  previous  segment  in  the 
path  for  each  node.  The  seventh  argument,  IDIR,  is  set  to  1 if  travel  is  from 
the  starting  node  or  to  -1  if  travel  is  toward  the  starting  node.  The  eighth 
argument,  MAXLVL,  is  the  maximum  number  of  levels  or  steps  allowed  in  the  path 
from  the  starting  node.  The  ninth  argument,  INSECT,  causes  the  path  to  be  re- 
stricted to  segments  in  section  INSECT.  If  INSECT  is  0,  all  segments  may  be 
used. 

The  items  in  arrays  DIST,  TIME,  NPREVN,  and  NPREVS  correspond  to  the 
noaes  in  array  NODNUM.  A double-dimensioned  array,  1ST,  and  two  single- 
dimensioned  arrays,  NIU  and  NIULOC,  are  used  to  build  paths  from  the  starting 
node.  The  1ST  array  is  used  to  construct  a binary  tree  that  corresponds  to 
the  paths  from  the  starting  node.  The  starting  node  is  the  root  of  the  tree. 

All  nodes  that  can  be  reached  in  one  step  (segment)  from  the  starting  node  are 
added  to  the  tree.  The  distances  to  these  nodes  are  recorded  in  the  NIU  array, 
indicating  that  these  nodes  are  in  use.  A large  negative  number  is  stored  as 
the  distance  to  the  starting  node.  Nodes  that  are  not  in  use  have  a 0 in  the 
corresponding  locations  of  the  NIU  array.  All  nodes  that  can  be  reached  in 
one  more  step  from  the  starting  node  are  now  added  to  the  tree.  If  any  node 
just  found  is  already  in  use,  and  if  its  distance  from  the  starting  point  is 
greater  the  second  time  it  is  found  than  the  first  time,  it  is  skipped.  If 
the  earlier  occurrence  of  the  node  has  the  greater  distance,  the  node  is  de- 
leted from  the  tree  and  all  nodes  reached  from  it  are  removed  from  the  in-use 
category.  When  a pass  .is  made  in  which  no  nodes  are  added  and  no  distances 
are  diminished,  the  procedure  is  finished.  (The  procedure  can  be  made  to 
terminate  after  a fixed  number  of  steps  by  using  the  desired  number  as  MAXLVL 
in  the  calling  sequence.) 

In  the  tree  representation,  the  root  is  the  starting  node.  A right 
branch  leads  to  a node  reached  on  the  next  step.  The  left  branches  connect 
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all  nodes  reached  on  the  same  step.  Three  pointers  are  associated  with  each 
element  in  the  tree:  a left  pointer,  a right  pointer,  and  a back  pointer.  The 
back  pointer  points  to  the  root  of  each  subtree  and  is  used  to  facilitate  un- 
linking elements  in  the  tree.  The  right  pointer  points  to  the  element  to  the 
right  in  the  tree;  this  pointer  has  a value  of  0 if  there  is  no  right  branch. 

The  left  pointers  point  to  left  branches.  When  there  is  no  left  branch,  the 
left  pointer  points  to  the  node  one  step  closer  to  the  starting  node.  NIULOC(I) 
gives  the  position  In  the  1ST  array  of  a node  with  number  NOONUM(I). 

Following  the  declarative  statements,  six  functions  are  defined. 
Function  INS  gives  the  section  number  of  a segment.  Functions  NU1  and  NU2  give 
the  starting  and  ending  nodes  of  a segment.  Function  LEN  gives  a segment's 
length.  Function  TIM  gives  the  time  required  to  travel  from  one  end  of  a seg- 
ment to  the  other.  Function  NWAY  gives  the  number  of  ways  of  travel  permitted 
on  a segment. 

Subroutine  TREE  begins  execution  by  computing  CW,  the  coefficient  of 
time  in  the  linear  combination  of  time  and  distance  minimized  as  paths  are  gen- 
erated from  the  starting  node.  The  line  number  of  the  starting  node  is  found 

by  function  IF IND  and  is  stored  In  variable  I.  If  I > 0,  control  transfers  to 

statement  570.  Otherwise,  an  error  message  is  printed  and  subroutine  EXIT  is 
called  to  terminate  the  program. 

Statements  570  through  600  clear  arrays  1ST,  NIL),  and  NIULOC.  At 

statement  600  the  fourth  column  of  array  1ST  is  set  to  the  line  number  of  the 

next  row.  These  numbers  are  used  as  links  for  unused  storage  in  the  1ST  array. 
They  will  be  overwritten  by  back  pointers  later  in  the  subroutine.  Variable 
NEXT,  obtained  from  1ST (1,4),  indicates  the  next  available  line  of  the  1ST 
array.  The  line  number  of  the  starting  node  is  stored  in  I ST (1,1).  The  left, 
right,  and  back  pointers  for  this  line  are  set  to  0.  Variable  CUM,  the 
weighted  function  of  time  and  distance,  is  set  to  0.  The  starting  node  is  de- 
signated in-use  by  the  storage  of  -1000  in  the  location  of  array  NIU  corres- 
ponding to  the  starting  node.  Variable  NMRT  is  set  to  1,  indicating  that  the 
current  activity  is  at  the  first  level  of  the  tree.  The  loop  through  statement 
940  will  build  up  to  MAXLVL  levels  of  the  tree.  Variable  NEWS,  a new  node  in- 
dicator, is  set  to  0.  Variable  K,  a pointer  to  the  current  line  of  the  1ST 
array,  is  set  to  1. 


32 


At  statement  620,  variable  NEXTK  Is  set  equal  to  the  right  pointer  of 
the  current  tree  entry.  If  NEXTK  Is  0,  control  transfers  to  statement  640. 
Otherwise,  the  level  counter,  NMRT,  is  incremented  by  1 . K is  set  equal  to 
NEXTK,  and  control  transfers  back  to  statement  620. 


Statement  640  is  reached  when  an  element  that  has  no  right  branch 
has  been  found.  At  statement  640,  if  the  level  number  is  equal  to  the  loop 
index,  J,  indicating  that  additional  nodes  should  be  added  to  the  tree,  con- 
trol transfers  to  statement  680.  Otherwise,  at  statement  660  NEXTK  is  set 
equal  to  the  left  pointer  of  the  current  tree  entry.  K is  set  equal  to  the 
absolute  value  of  NEXTK.  If  NEXTK  is  negative,  indicating  the  absence  of  a 
left  branch,  control  transfers  to  statement  670.  If  NEXTK  is  0,  indicating 
that  the  NEXTK  pointer  is  back  to  the  root  of  the  entire  tree,  control  trans- 
fers to  statement  930.  If  NEXTK  is  positive,  indicating  the  presence  of  a 
left  branch,  control  transfers  to  statement  620. 

At  statement  670,  the  level  counter,  NMRT,  is  decremented  by  1. 
Control  transfers  to  statement  660.  At  statement  680,  variable  I FIRST  is 
set  to  1 to  indicate  that  the  next  node  added  to  the  tree  will  be  a right 
branch  rather  than  a left  branch.  The  line  number,  N,  of  the  node  being 
built  upon  is  obtained  from  IST(K.l).  Variable  LAST  is  set  equal  to  K,  the 
current  node's  line  in  the  tree. 

Variable  KSEG  is  set  equal  to  the  number  of  neighboring  segments  for 
the  node.  The  loop  through  statement  900  will  examine  these  neighboring  seg- 
ments and  add  to  the  tree  nodes  that  are  one  step  farther  from  the  base  node. 
The  number  of  a neighboring  segment  is  retrieved  and  stored  in  variable  KK. 

If  a section  restriction  on  segments  is  present,  and  if  the  segment  is  not  in 
the  required  section,  control  transfers  to  statement  900.  Otherwise,  the  num- 
ber of  the  ending  node  of  the  segment  is  stored  in  variable  L.  If  the  starting 
node  of  the  segment  is  the  node  already  in  the  tree,  and  if  travel  from  the 
starting  to  the  ending  node  is  permissible,  control  transfers  to  statement  700. 
Otherwise,  if  the  ending  node  is  not  in  the  tree,  or  if  travel  from  the  start- 
ing to  ending  node  is  not  allowed,  control  transfers  to  statement  900.  Other- 
wise, the  number  of  the  starting  node  is  stored  in  variable  L. 
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At  statement  700  the  line  number  of  node  L is  obtained  by  function 
IF IND  and  is  stored  in  variable  L.  The  cumulative  value  of  the  linear  combina- 
tion of  time  and  distance  to  the  node  is  stored  in  variable  CUM.  The  location 
in  the  NIU  array  corresponding  to  the  node  is  examined  for  the  status  of  the 
node.  If  NIU(L)  is  negative,  control  transfers  to  statement  900.  If  It  is  0, 
the  node  is  not  in  use,  and  control  transfers  to  statement  840.  If  it  is  posi- 
tive, the  node  is  already  in  use,  and  control  transfers  to  statement  720. 

At  statement  720,  CUM  is  compared  to  the  value  stored  in  array  NIU. 

If  the  current  value  (CUM)  is  greater  than  or  equal  to  the  previous  value,  con- 
trol transfers  to  statement  900.  If  not,  the  previous  occurrence  of  the  node 
is  disconnected  from  the  tree.  The  location  in  the  tree  of  the  previous  occur- 
rence of  the  node  is  obtained  from  array  NIUL0C  and  is  stored  in  variable  LB. 
The  back  and  left  pointers  are  obtained  from  the  tree  and  are  stored  in  vari- 
ables LA  and  LC.  If  the  node  to  be  disconnected  was  not  the  last  node  added 
to  the  tree,  control  transfers  to  statement  730.  Otherwise,  variable  LAST  is 
set  equal  to  the  back  pointer.  If  the  node  is  the  first  node  of  its  level 
added  to  the  node  of  the  previous  level,  variable  IFIRST  is  set  to  1 . At  state- 
ment 730  this  test  is  repeated  (although  different  FORTRAN  variables  are  used); 
if  the  test  is  true,  control  transfers  to  statement  740.  Otherwise,  the  left 
pointer  of  the  node  indicated  by  the  back  pointer  is  replaced  by  the  left 
pointer  of  the  current  node.  Control  transfers  to  statement  760.  At  state- 
ment 740  the  right  pointer  of  the  node  indicated  by  the  back  pointer  is  re- 
placed by  either  0 or  the  current  node's  left  pointer,  whichever  is  larger. 

At  statement  760,  if  the  left  pointer  of  the  current  node  is  positive,  the 
back  pointer  of  the  node  indicated  by  the  current  node's  left  pointer  is  re- 
placed by  the  current  node's  back  pointer.  The  left  pointer  of  the  current 
node  is  set  to  0.  Variable  LC  is  set  to  the  current  node's  line  in  the  tree. 

The  statements  from  700  through  840  remove  the  portion  of  the  tree 

built  on  the  node  just  disconnected.  At  statement  780,  variable  LB  is  set 
equal  to  variable  LC.  At  statement  800,  LC  is  set  equal  to  the  right  pointer 
of  the  node  at  line  LB.  If  the  right  pointer  is  positive,  control  transfers 
to  statement  780.  If  not,  at  statement  820  the  line  number  of  the  node  in 

array  N0DNUM  is  obtained  from  the  tree  and  stored  in  variable  NA.  The  node's 

entries  in  arrays  NIU  and  NIUL0C  are  set  to  0.  The  back  pointer  is  set  equal 
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to  the  line  number  of  the  next  free  storage  in  the  tree.  Variable  NEXT  is  set 
equal  to  the  line  number  of  the  node  currently  being  disconnected.  The  left 
pointer  is  stored  in  variable  LC,  and  LB  is  set  equal  to  the  absolute  value  of 
LC.  If  LC  is  negative,  control  transfers  to  statement  820.  Control  transfers 
to  statement  840  if  LC  is  zero,  or  to  statement  800  if  it  is  positive. 

A new  node  is  added  to  the  tree  by  the  processing  completed  between 
statements  840  and  900.  The  new  node  indicator,  NEWS,  is  set  to  1.  The  node's 
line  number  in  the  N0DNUM  array,  L,  is  stored  in  the  next  available  line  of  the 
tree.  The  current  line's  left  pointer  is  set  equal  to  the  left  pointer  of  the 
previous  line.  The  line  number  of  the  previous  node  in  the  path  is  obtained 
from  the  tree  and  stored  in  array  NPREVN.  The  segment  number,  KK,  is  stored 

in  the  NPREVS  array.  The  right  pointer  is  set  to  0.  The  location  of  the  next 

free  storage  in  the  tree  is  taken  from  the  current  entry's  back  pointer  and  is 
stored  in  variable  NEWNXT.  The  back  pointer  is  set  equal  to  the  line  number 
of  the  previous  entry  in  the  tree.  The  distance  and  time  to  the  current  node 
are  computed.  The  NIU  array  entry  corresponding  to  the  current  node  is  set 
equal  to  CUM.  The  NIUL0C  array  entry  corresponding  to  the  current  node  is  set 

equal  to  NEXT,  the  line  number  of  the  node  in  the  tree.  If  I FIRST  is  equal  to 

1,  control  transfers  to  statement  860.  Otherwise,  the  previous  entry's  left 
pointer  is  set  equal  to  the  current  node's  line  number.  Control  transfers  to 
statement  880. 

At  statement  860  the  previous  entry's  right  pointer  is  set  equal  to 
the  current  node's  line  number.  The  current  node's  left  pointer  is  set  equal 
to  the  negative  of  the  previous  entry's  line  number.  The  line  number  of  the 
previous  node,  NPREVN(L),  is  obtained  from  the  previous  entry  in  the  tree. 
IFIRST  is  set  equal  to  0. 

At  statement  880,  variable  LAST  is  set  equal  to  the  current  line 
number  in  the  tree.  Variable  NEXT  is  set  equal  to  NEWNXT,  the  line  number  of 
the  next  free  storage  in  the  tree.  Statement  900  marks  the  end  of  the  loop 
that  scans  segments  connected  to  each  node  in  the  previous  level  of  the  tree. 

A search  is  begun  for  another  node  in  the  tree  to  which  additional 
nodes  can  be  added.  If  any  nodes  have  been  added  at  the  current  level,  K is 
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set  equal  to  the  negative  of  the  left  pointer  of  the  last  entry.  At  statement 
920,  K is  replaced  by  the  left  pointer  of  the  tree  entry  at  line  K.  if  K is 
positive,  control  transfers  to  statement  620.  Otherwise,  the  level  counter, 
NMRT,  is  decremented  by  1.  The  sign  of  K is  changed.  If  K is  positive,  con- 
trol transfers  to  statement  920.  Otherwise,  NMRT  is  incremented  by  1 . If  no 
new  nodes  have  been  added  to  the  tree,  control  transfers  out  of  the  loop  on 
statement  940  to  statement  960.  Statement  940  marks  the  end  of  the  loop  that 
extends  the  level  of  the  tree.  Following  statement  960,  control  returns  to 
the  calling  program. 

1.  Subroutine  C0N2ST 

Subroutine  C0N2ST  connects  pieces  of  a section  by  two-segment  paths 
where  possible.  The  shortest  two-segment  path  between  pieces  is  used.  There 
are  two  arguments.  The  first,  NTRIP,  gives  the  section  or  trip  number.  The 
second,  NC0N,  is  the  count  of  remaining  pieces  requiring  more  than  a two-step 
path  to  connect  them.  NC0N  is  computed  and  returned  by  subroutine  C0N2ST. 

ISTG( 1 0, i ) contains  the  section  and  piece  information  for  segment  i.  The  sec- 
tion number  identifies  one  piece  of  the  section,  and  multiples  of  1000  are 
added  to  the  section  number  to  indicate  the  remaining  pieces. 

NS,  a count  of  two-step  connections,  is  set  to  0.  The  loop  on  state- 
ment 10  sets  to  0 the  arrays  used  for  piece  numbers,  segment  numbers,  and  con- 

nection lengths.  The  loop  through  statement  140  examines  each  node  as  a possi- 
ble midpoint  for  a two-step  path  between  pieces.  The  node  number  is  stored  in 
variable  NODE.  The  count  of  segments  connected  to  this  node  is  stored  in  LIM. 
If  only  one  segment  is  connected  to  the  node,  control  transfers  to  statement 
140.  Variable  LIMM1  is  set  equal  to  one  less  than  ti  ° number  of  neighboring 
segments.  The  packed,  neighboring-segment  numbers  are  stored  in  NBRS1 . 

The  loop  through  statement  130  examines  all  but  the  last  neighboring 
segment  as  possible  first  segments  of  the  two-segment  path.  The  last  segment 

will  be  considered  only  as  a possible  second  step  in  a later  loop.  A segment 

number  is  retrieved  and  stored  in  variable  ISEG.  If  the  segment  is  in  section 
NTRIP,  control  transfers  to  statement  140.  Otherwise,  the  line  number,  L,  of 
the  node  at  the  other  end  of  the  segment  is  found.  The  count  of  segments 
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bounding  the  second  node  is  stored  in  variable  LIM2.  The  packed,  bounding- 
segment  numbers  are  stored  in  NBRS2. 

The  loop  through  statement  20  examines  each  of  the  segments  connected 
to  the  second  node.  If  a segment  is  found  in  section  NTRIP,  control  transfers 
to  statement  30.  Otherwise,  control  transfers  to  statement  130. 

The  loop  through  statement  120  examines  the  remaining  segments  bound- 
ing the  first  node.  The  segment  number  is  retrieved  and  stored  in  variable 
JSEG.  If  the  segment  Is  In  the  current  section,  control  transfers  to  statement 
140.  Otherwise,  the  line  number  of  the  other  node  is  found  and  Is  stored  in 
variable  L.  The  count  of  segments  bounding  this  second  node  Is  stored  in  LIM2. 
The  packed  segment  numbers  are  stored  in  NBRS2. 

The  loop  through  statement  40  examines  each  segment  bounding  the 
second  node.  If  any  of  the  segments  are  in  the  current  section,  control  trans- 
fers to  statement  50.  Otherwise,  control  transfers  to  statement  120. 

At  statement  50,  the  piece  numbers  connected  by  the  two-step  path 
are  compared.  If  the  piece  numbers  are  equal,  control  transfers  to  statement 
120.  If  JSECT  is  larger  than  KSECT,  control  transfers  to  statement  60.  Other- 
wise, control  transfers  to  statement  70.  At  statement  60  the  values  of  JSECT 
and  KSECT  are  Interchanged  so  that  the  smaller  piece  number  will  be  In  JSECT. 

At  statement  70,  if  either  of  the  segments  is  a two-way  segment,  control  trans- 
fers to  statement  80,  where  the  length  of  the  two-segment  path  is  stored  in 
variable  DIST.  Otherwise,  If  the  two  one-way  segments  point  toward  or  away 
from  one  other,  control  transfers  to  statement  120. 

If  no  two-step  paths  have  yet  been  found,  control  transfers  to  state- 
ment 100.  Otherwise,  the  loop  on  statement  90  examines  arrays  NP1  and  NP2  to 
see  whether  a path  connecting  pieces  JSECT  and  KSECT  has  already  been  found. 

If  a shorter  connection  has  been  saved,  control  transfers  to  statement  120.  If 
the  new  connection  is  shorter,  control  transfers  to  statement  110.  If  no  path 
connecting  the  pieces  has  been  saved  as  yet,  the  loop  Is  exited. 
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At  statement  100,  the  number  of  entries  saved  Is  Incremented  by  1. 

The  five  FORTRAN  statements  starting  at  statement  110  save  the  piece  numbers, 
the  segment  numbers,  and  the  length  of  the  two-step  path.  Statement  120  marks 
the  end  of  the  loop  that  seeks  the  second  segment  In  the  path.  Statement  130 
marks  the  end  of  the  loop  that  seeks  the  first  segment  of  the  path.  Statement 
140  marks  the  end  of  the  loop  on  the  node  connecting  the  two  segments. 

The  remaining  statements  in  the  subroutine  replace  the  higher  piece 
number  by  the  lower  for  all  segments  in  pieces  that  have  been  connected  by  two- 
step  paths.  A count  of  the  pieces,  NCON,  is  initially  set  to  0.  The  loop 
through  statement  170  examines  each  of  the  connections.  The  smaller  piece  num- 
ber from  the  NP1  array  is  stored  in  variable  NNEW,  and  the  larger  piece  number 
from  array  NP2  is  stored  in  variable  NOLD.  The  smaller  piece  number  is  stored 
in  the  ISTG  array  for  each  of  the  segments.  If  the  two  piece  numbers  are  equal, 
control  transfers  to  statement  170.  (The  numbers  cannot  be  equal  on  the  first 
pass  through  the  loop,  but  may  be  equal  after  the  first  pass.)  Otherwise, 

IFOUND  is  set  to  0. 

The  loop  through  statement  150  examines  all  of  the  segments  in  the 
map  description  and  replaces  the  larger  piece  number  by  the  smaller  wherever 
it  occurs.  If  any  numbers  are  replaced,  IFOUND  is  set  to  1.  After  statement 
150,  the  number  of  pieces,  NCON,  is  incremented  by  IFOUND.  If  the  last  connec- 
tion is  being  processed,  control  transfers  to  statement  170.  Otherwise,  IP1 
is  set  equal  to  the  line  number  of  the  next  connection. 

The  loop  through  statement  160  examines  the  remaining  connections 
and  replaces  occurrences  of  piece  number  NOLD  by  piece  number  NNEW.  The  num- 
bers are  rearranged,  if  necessary,  to  make  the  number  in  the  NP1  array  smaller 
than  the  number  in  the  NP2  array.  Following  the  loop  through  statement  170, 
control  returns  to  the  calling  program. 

m.  Subroutine  FNDPTH 

Subroutine  FNDPTH  finds  up  to  five  shortest  paths  connecting  pieces 
of  a given  trip  to  a given  node  in  another  piece  of  the  trip.  FNDPTH  has  ten 
arguments.  The  first  argument,  NTRIP,  is  the  trip  or  section  number.  The 
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second  argument,  NORG,  is  the  number  of  the  node  at  which  the  paths  must  end. 

The  third  argument,  IP,  is  the  number  of  the  piece  containing  node  NORG.  The 
fourth  argument,  NDSV,  is  the  number  of  paths  saved.  The  fifth  argument,  LNF, 
is  the  line  number  of  the  node  that  starts  the  path.  The  sixth  argument,  NPF, 
is  the  number  of  the  piece  in  which  the  starting  node  lies.  The  seventh  argu- 
ment, LNT,  and  the  eighth  argument,  NPT,  are  arrays  containing  the  line  number 
of  the  terminal  node  and  the  number  of  the  piece  in  which  it  lies.  The  ninth 
argument,  TIM,  Is  an  array  of  travel  times  for  the  paths.  The  tenth  argument, 
IPATH,  is  a double-dimensioned  array  that  contains  the  numbers  of  the  segments 
in  up  to  five  paths  to  node  NORG. 

The  maximum  number  of  segments  allowed  in  a path,  MAXSTG,  is  eval- 
uated by  using  the  LOCF  function  on  two  words  in  the  IPATH  array  whose  second 
subscripts  differ  by  one.  The  loop  through  statement  10  seeks  the  line  number 
of  node  NORG  and  sets  to  0 arrays  DIST,  TIME,  LPREVN,  and  NPREVS.  The  loop 
through  statement  40  sets  the  piece  and  line  number  arrays,  NPT  and  LNT,  to 
the  values  appropriate  for  node  NORG.  The  embedded  loop  on  statement  30  clears 
the  IPATH  array. 

Subroutine  TREE  is  called  to  build  a tree,  using  travel  time  between 
node  NORG  and  each  other  node  as  the  variable  to  be  minimized.  The  direction 
of  travel  Is  away  from  node  NORG. 

NDSV,  the  number  of  nodes  found  in  other  pieces,  is  set  to  0.  The 
loop  through  statement  90  will  examine  the  piece  assignment  of  each  segment. 
Variable  JP  is  set  to  the  piece  number  of  the  segment.  If  the  segment  is  not 
in  the  current  section,  NTRIP,  or  if  it  is  in  piece  IP,  control  transfers  to 
statement  90.  Variable  L is  set  equal  to  the  line  number  of  the  end-point 
node  closer  to  node  NORG.  Variable  T is  set  equal  to  the  travel  time  to  this 
node. 

Variable  JJ,  a pointer  to  the  line  at  which  nodes  will  be  saved,  is 
initially  set  to  1 . If  no  nodes  have  been  saved  as  yet,  control  transfers  to 
statement  60.  Otherwise,  the  loop  through  statement  50  Is  executed.  This 
loop  examines  the  piece  number  of  each  saved  node.  If  a node  is  found  in  piece 
JP,  the  travel  time  to  that  node  is  compared  with  the  travel  time  to  the  current 
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node.  If  the  node  already  saved  is  reached  in  a shorter  time  than  is  the  cur- 
rent node,  control  transfers  to  statement  90.  Otherwise,  control  transfers  to 
statement  70.  When  no  nodes  in  the  table  are  in  piece  JP,  control  continues 
following  the  loop  through  statement  50.  Variable  JJ  Is  set  to  one  more  than 
the  number  of  nodes  already  saved.  If  the  maximum  permissible  number  of  nodes 
has  not  been  saved,  control  transfers  to  statement  60.  If  the  maximum  has  been 
reached,  and  if  the  travel  time  to  the  current  node  is  greater  than  or  equal  to 
the  largest  travel  time  already  saved,  control  transfers  to  statement  90. 
Otherwise,  JJ  is  set  equal  to  the  line  number  of  the  largest  time  in  the  table. 
Control  transfers  to  statement  70.  At  statement  60,  the  number  of  saved  nodes 
is  Incremented  by  1.  At  statement  70,  the  piece  number  of  the  current  segment 
is  saved.  The  line  number  of  the  current  node  and  the  travel  time  to  that  node 
are  saved. 


The  loop  through  statement  80  examines  the  travel  times  to  each  of 
the  saved  nodes.  The  largest  travel  time  is  saved  in  variable  TMAX,  and  the 
line  number  of  this  entry  is  saved  in  variable  LTMAX.  Statement  90  marks  the 
end  of  the  loop  that  searches  for  segments  in  pieces  of  the  current  trip. 

The  loop  through  statement  120  retrieves  the  paths  from  each  saved 
node  to  node  N0RG.  Variable  NP,  the  number  of  segments  in  the  path,  is  ini- 
tially set  to  0.  Variable  L is  set  equal  to  the  line  number  in  the  NODNUM 
array  of  the  saved  node.  At  statement  100,  NP  is  incremented  by  1.  The  seg- 
ment number  of  the  next  segment  in  the  path  is  saved,  both  in  the  IPATH  array 
and  as  variable  J.  The  line  number  of  the  next  node  in  the  path  is  stored  in 
variable  L.  If  the  segment  is  in  piece  IP  [which  has  been  saved  in  NPT(I)], 
or  if  the  next  node  is  node  N0RG,  control  transfers  to  statement  120.  If  not, 
and  if  the  number  of  segments  in  the  path  is  fewer  than  the  maximum  allowed, 
control  returns'  to  statement  100.  Otherwise,  an  error  message  is  printed  in- 
dicating that  the  maximum  number  of  steps  has  been  exceeded.  The  first  15 
segments  in  the  path  are  also  printed.  System  subroutine  EXIT  is  called  to 
terminate  the  program.  Statement  120  marks  the  end  of  the  loop  that  obtains 
paths  from  the  saved  nodes.  Control  returns  to  the  calling  program. 


n.  Subroutine  CONNST 
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Subroutine  CONNST  connects  pieces  of  section  NTRIP  that  could  not  be 
connected  by  one-  or  two-step  paths.  It  has  two  arguments:  the  trip  number, 
NTRIP,  and  the  number  of  pieces,  NPIECE,  remaining  at  the  completion  of  subrou- 
tine CONNST. 

Variable  ITER,  the  count  of  iterations  through  subroutine  CONNST,  is 
set  to  0.  The  loop  through  statement  10  seeks  a segment  that  is  in  section 
NTRIP.  When  a segment  is  found,  the  starting  node  number  is  saved  in  variable 
NORG  and  control  transfers  to  statement  20.  At  statement  20,  ITER  is  incre- 
mented by  1 . 

Subroutine  FNDPTH  is  called  to  find  up  to  five  nodes  in  other  pieces 
of  section  NTRIP.  The  number  of  nodes  found  is  stored  in  variable  LIM.  Vari- 
able NDC  is  set  to  1 . 

The  loop  through  statement  40  uses  each  of  the  nodes  found  by  the 
first  call  to  FNDPTH  as  starting  points  for  additional  calls  to  FNDPTH.  The 
paths  found  in  these  subsequent  calls  are  saved  in  array  IPATH.  The  total 
number  of  paths  is  saved  in  variable  NDC. 

The  loop  through  statement  60  searches  the  array  of  path  times.  The 
number  of  the  path  with  the  smallest  time  is  saved  in  variable  LTM  If  no 
paths  have  been  found  between  pieces  of  section  NTRIP,  an  error  message  is 
printed  and  control  returns  to  the  calling  program.  Otherwise,  at  statement 
80  the  time  for  the  shortest  path  is  set  to  a very  large  number.  The  piece 
numbers  connected  by  the  path  are  stored  in  variables  NPS  and  NPL.  The  number 
of  pieces,  NPIECE,  is  decremented  by  1.  The  loop  through  statement  90  replaces 
each  occurrence  of  the  larger  piece  number  in  arrays  NPT  and  NPF  by  the  smaller 
piece  number.  The  loop  on  statement  100  replaces  the  larger  piece  number  by 
the  smaller  piece  number  in  the  segment  data  array,  ISTG.  The  loop  through 
statement  110  changes  the  section  assignments  of  segments  on  the  path  between 
the  two  pieces  to  that  of  the  smaller  numbered  piece. 
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The  loop  through  statement  130  examines  the  piece  numbers  in  the  NPF 
and  NPT  arrays.  If  any  numbers  remain  that  are  different  from  the  smaller 
piece  number,  control  transfers  to  statement  50.  Otherwise,  if  fewer  than  10 
iterations  have  been  performed  and  if  more  than  one  piece  remains  in  the  sec- 
tion, control  transfers  to  statement  20.  Otherwise,  if  10  or  more  iterations 
have  been  performed,  an  error  message  is  printed.  Control  returns  to  the  call- 
ing program. 

o.  Subroutine  CONNECT 

Subroutine  CONNECT  connects  pieces  of  a section.  It  has  two  argu- 
ments. The  first,  NTRIP,  is  the  section  number.  The  second,  NPIECE,  is  the 
number  of  separate  pieces  remaining  after  the  subroutine  has  connected  all 
the  pieces  it  can. 

Variable  II  and  the  piece  count,  NPIECE,  are  each  set  equal  to  1. 

The  section  number,  NTRIP,  is  stored  in  variable  ITRIP.  The  loop  on  statement 
10  clears  arrays  used  by  subroutine  TREE  for  path  information.  At  statement 
20  the  loop  through  statement  30  is  started.  Variable  II  is  set  equal  to  the 
current  value  of  the  loop  index,  I.  If  the  Ith  segment  is  in  section  ITRIP, 
control  transfers  to  statement  40.  If  the  loop  is  completed  and  no  segments 
in  section  ITRIP  are  found,  program  execution  terminates  with  a STOP  30  state- 
ment. 


At  statement  40  the  number  of  the  starting  node  of  the  segment  is 
stored  in  variable  NORG.  Subroutine  TREE  is  called  twice,  first  to  build  paths 
using  only  segments  in  section  ITRIP  starting  at  node  NORG,  and  again  to  find 
paths  to  node  NORG  using  only  segments  in  section  ITRIP.  The  paths  are  built 
to  minimize  the  distance  to  or  from  node  NORG. 

Variable  DJ  is  initially  set  to  0.  Variable  II  is  incremented  by  1. 
If  II  is  greater  than  the  number  of  segments,  control  transfers  to  statement 
80.  Otherwise,  a loop  is  begun  through  statement  70  to  mark  segments  in  other 
pieces  of  the  current  section.  If  the  segment  designated  by  the  loop  index  is 
not  in  the  current  section,  control  transfers  to  statement  70.  Otherwise,  the 
line  number  of  the  starting  node,  LINE,  is  found.  If  a path  exists  between 
node  NORG  and  the  current  node,  control  transfers  to  statement  70.  Otherwise, 
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the  section  number  of  the  segment  is  incremented  by  1000.  Variable  OJ  is  set 
to  1,  indicating  that  an  unconnected  node  has  been  found.  If  DJ  is  0 when  the 
loop  through  statement  70  is  completed,  control  transfers  to  statement  80. 
Otherwise,  NPIECE  is  incremented  by  1,  ITRIP  is  increased  by  1000,  and  control 
returns  to  statement  20. 

Following  statement  80,  if  NPIECE  equals  1,  control  returns  to  the 
calling  program.  If  not,  the  loop  through  statement  130  is  executed.  The  loop 
examines  each  segment  to  determine  whether  it  connects  two  pieces  of  section 
NTRIP . The  section  number  is  retrieved  and  stored  in  variable  JSECT.  If  the 
segment  is  in  section  NTRIP,  control  transfers  to  statement  130.  Otherwise, 
the  starting  and  ending  node  numbers  are  saved  in  array  NNN. 

The  loop  through  statement  110  examines  each  of  the  nodes.  The  count 
of  segments  connected  to  the  node  is  obtained  from  the  NUMNBR  array  and  is 
stored  in  variable  LIM.  If  LIM  is  less  than  or  equal  to  1,  control  transfers 
to  statement  130.  Otherwise,  the  loop  through  statement  100  is  executed.  This 
loop  examines  each  neighboring  segment.  If  the  segment  is  the  same  as  that  in- 
dicated by  the  index  of  the  loop  through  statement  130,  control  transfers  to 
statement  100.  Otherwise,  the  section  number  of  the  segment  is  stored  in  array 
ISECT.  If  the  segment  is  in  section  NTRIP,  control  transfers  to  statement  110. 
Otherwise,  control  passes  to  the  end  of  the  loop.  Following  the  loop  through 
statement  100,  control  transfers  to  statement  130.  Following  statement  110,  if 
the  two  end-point  nodes  are  in  the  same  section,  control  transfers  to  statement 
130.  Otherwise,  the  segment  is  assigned  a section  number  formed  by  concaten- 
ating the  piece  numbers  of  the  end-point  nodes.  The  smaller  number  precedes 
the  larger.  Statement  130  marks  the  end  of  the  loop  that  searches  for  single- 
segment connections  between  pieces. 

After  the  single-segment  connections  have  been  marked,  the  loop 
through  statement  160  assigns  the  lower  piece  number  to  all  segments  in  the 
connected  pieces.  The  connecting  segment  is  also  assigned  to  the  piece. 

The  loop  through  statement  170  counts  the  number  of  separate  pieces. 


The  result  is  saved  as  NPIECE.  If  NPIECE  is  less  than  or  equal  to  1,  control 
returns  to  the  calling  program.  Otherwise,  subroutine  C0N2ST  is  called  to 
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make  two-step  connections  between  separate  pieces  of  section  NTRIP.  Only  the 
shortest  connection  between  pieces  is  made.  If  the  section  has  been  merged 
into  one  piece,  control  returns  to  the  calling  program.  Otherwise,  execution 
continues  within  subroutine  CONNECT. 

The  loop  on  statement  210  clears  arrays  NNSN  and  NNTN.  The  loop 
through  statement  240  examines  each  segment.  The  section  number  is  stored 
in  variable  K.  If  the  section  is  not  equal  to  NTRIP,  control  transfers  to 
statement  240.  Otherwise,  the  loop  through  statement  230  is  executed.  This 
loop  saves  the  piece  number  and  accumulates  the  number  of  houses  in  the  piece. 
Following  statement  240,  variables  NMAX  and  NC0N  are  set  to  0.  The  loop 
through  statement  270  examines  the  number  of  houses  in  each  piece.  If  the  num- 
ber of  houses  is  greater  than  0,  control  transfers  to  statement  260.  If  not, 
the  piece  number  is  stored  in  variable  K.  Variable  NC0N  is  incremented  by  1. 
The  loop  on  statement  250  removes  from  section  NTRIP  all  segments  in  a piece 
that  contains  no  houses.  Control  transfers  to  statement  270.  Statement  260 
saves  as  variable  NMAX  the  largest  piece  number  encountered  in  the  loop  through 
statement  270. 

The  loop  through  statement  280  locates  the  line  in  the  NNFN  array 
having  the  lowest  piece  number.  This  number  should  be  the  section  number, 
NTRIP.  If  so,  control  transfers  to  statement  290.  If  not,  program  execution 
terminates  at  a ST0P270  statement. 

At  statement  290  the  house  count,  NNTN,  is  checked.  If  NNTN  is  posi- 
tive, control  transfers  to  statement  310.  Otherwise,  the  loop  through  state- 
ment 300  replaces  the  piece  numbers  of  segments  in  the  highest  numbered  piece 
with  the  section  number. 

At  statement  310,  the  number  of  deleted  pieces,  NC0N,  is  subtracted 
from  the  total  piece  count.  If  only  one  piece  remains,  control  returns  to  the 
calling  program.  Otherwise,  subroutine  C0NNST  is  called  to  connect  all  of  the 
remaining  pieces.  Control  then  returns  to  the  calling  program. 


Subroutine  DISCON  removes  from  a section  all  segments  that  form 
dead-end  streets  on  which  no  collection  is  required.  The  subroutine  has  four 
arguments.  The  first,  NTRIP,  gives  the  section  number.  The  second,  KN,  gives 
the  number  of  nodes  in  the  section.  The  third,  NODPTR,  is  an  array  of  pointers 
to  the  node  numbers.  The  fourth,  NODORD,  is  an  array  of  node  orders. 

Variable  IDEL,  a count  of  deleted  nodes,  is  set  to  0.  At  statement 
10  the  repeat  indicator,  IRP,  is  set  to  0.  The  loop  through  statement  60  ex- 
amines each  of  the  nodes  in  the  section.  If  the  node  order  is  not  equal  to  1, 
control  transfers  to  statement  60.  Otherwise,  variable  LINE  is  set  equal  to 
the  line  number  of  the  node.  The  number  of  segments  connected  to  the  node  is 
obtained  from  array  NUMNBR  and  is  stored  in  variable  LIM.  A word  coni,  ining 
the  packed  segment  numbers  is  obtained  from  array  NBRSEG  and  is  stored  in  vari- 
able NBRS. 

The  loop  through  statement  20  examines  the  section  assignment  of  each 
neighboring  segment.  When  a segment  is  found  in  the  current  section,  NTRIP, 
control  transfers  to  statement  30.  If  no  segments  are  found  in  the  current 
section,  an  error  exists  in  the  program.  An  error  message  is  printed,  and  the 
program  terminates  with  a call  to  DUMP.  Otherwise,  at  statement  30,  if  the 
segment  requires  collection,  control  transfers  to  statement  60.  If  collection 
is  not  required,  the  segment  is  assigned  to  section  zero.  The  order  of  the 
node  connected  to  the  segment  is  set  to  0.  The  count  of  deleted  nodes,  IDEL, 
is  incremented  by  1.  Variable  NODE  is  set  equal  to  the  number  of  the  node  at 
the  other  end  of  the  segment.  The  loop  through  statement  40  searches  for  this 
node.  If  the  node  is  found,  control  transfers  to  statement  50.  If  not,  an 
error  exists  in  the  program.  An  error  message  is  printed,  and  program  execu- 
tion terminates  with  a call  to  DUMP. 

At  statement  50,  the  order  of  the  node  is  decreased  by  1.  If  this 
node  has  already  been  examined  in  the  loop  through  statement  60,  and  if  it  now 
has  an  order  equal  to  1,  the  repeat  indicator,  IRP,  is  set  to  1.  Statement 
60  marks  the  end  of  the  loop  that  examines  the  nodes  in  the  section.  If  the 
repeat  indicator  is  equal  to  1,  control  transfers  back  to  statement  10. 
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When  no  other  order-one  nodes  remain,  the  loop  through  statement  70 
closes  up  the  empty  entries  in  the  NODPTR  and  N0D0RD  arrays.  The  node  count, 
KN,  is  decreased  by  IDEL.  Control  returns  to  the  calling  program. 

q.  Subroutine  EPXP 

Subroutine  EPXP  selects  nodes  that  are  suitable  points  for  entry 
into  a section  from  the  landfill  or  garage  and  nodes  that  are  suitable  points 
from  which  to  leave  the  section  for  the  landfill.  The  subroutine  has  four 
arguments.  The  first,  KN,  gives  the  number  of  nodes  in  the  section.  The  sec- 
ond, NODPTR,  is  an  array  of  pointers  to  the  node  numbers.  The  third,  N0D0RD, 
is  an  array  of  node  orders.  The  fourth,  I,  is  the  section  number. 

Variable  KLIM,  the  maximum  number  of  entry  or  exit  nodes  to  be  found, 
is  computed  on  the  basis  of  the  number  of  nodes  in  the  section.  KLIM  will  have 
a value  from  5 through  10.  Variable  KKEEP  is  set  to  6.  Variable  KLPK  is  set 
to  the  number  of  nodes  or  to  the  sum  of  KLIM  and  KKEEP,  whichever  is  smaller. 

Data  are  read  from  file  TAPE7  using  a BUFFER  IN  statement.  These 
data  give  the  paths  from  each  node  to  the  landfill  and  to  the  garage.  TAPE7 
is  rewound.  The  loop  on  statement  240  clears  the  arrays  used  for  the  entry- 
and  exit-node  numbers.  The  loop  through  statement  290  finds  three  types  of 
suitable  nodes:  entry  nodes  for  paths  coming  from  the  landfill,  entry  nodes 
for  paths  coming  from  the  garage,  and  exit  nodes  for  paths  going  to  the  land- 
fill. The  loop  on  statement  250  moves  the  travel  times  from  the  garage  or 
landfill  into  the  DISTS  array  for  each  node  in  the  section.  Subroutine  SHLSRT3 
is  called  to  sort  the  times  into  increasing  order.  The  NODPTR  and  N0D0RD  ar- 
rays are  carried  along  during  the  sort. 


Variable  JN,  the  count  of  saved  nodes,  is  set  initially  to  0.  The 
loop  through  statement  280  examines  the  KLPK  nodes  closest  in  time  to  the 
garage  or  landfill.  The  line  number  of  a node  is  stored  in  variable  LINE.  The 
node  number  is  stored  in  variable  NS.  If  no  nodes  have  been  saved,  control 
transfers  to  statement  270.  Otherwise,  the  double  loop  through  statement  260 
examines  the  path  from  the  current  node  to  the  garage  or  landfill.  If  the 
path  includes  a node  already  saved  as  an  entry  or  exit  point,  control  transfers 
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to  statement  280.  Otherwise,  the  count  of  saved  nodes  is  incremented  at  state- 
ment 270.  The  node  number  is  stored  in  the  NTF  array.  Statement  280  marks  the 
end  of  the  loop  that  saves  the  entry  or  exit  nodes.  Statement  290  is  the  end 
of  the  loop  on  the  three  types  of  nodes.  At  statement  300  a message  is 
printed  giving  the  entry  and  exit  nodes  for  the  section.  Control  returns  to 
the  calling  program. 

r.  Subroutine  CLOSE  1 

Subroutine  CL0SE1  adds  to  a section  a path  back  to  the  section  from 
an  order-one  node.  The  subroutine  has  five  arguments.  The  first,  NTRIP,  is 
the  section  number.  The  second,  KN,  is  the  count  of  nodes  in  the  section.  The 
third,  KNX,  is  the  node  count  after  new  paths  have  been  added  to  the  section. 
The  fourth  argument,  NODPTR,  is  an  array  of  pointers  to  the  node  numbers.  The 
fifth  argument,  NODORD,  is  an  array  of  node  orders. 

Variable  KNX  is  set  equal  to  KN.  If  the  trip  number  of  the  current 
call  is  unequal  to  the  trip  number  of  the  previous  call,  or  if  no  path  is  cur- 
rently saved  in  array  IPSSV,  control  transfers  to  statement  8.  Otherwise,  the 
path  in  IPSSV  and  any  additional  segments  artificially  created  for  the  path  are 
set  to  0.  At  statement  8,  the  number  of  steps  in  the  path  is  set  to  0.  The 
current  section  number  is  saved  in  variable  NTRSV. 


The  loop  through  statement  200  examines  each  node  in  the  map  descrip- 
tion, seeking  order-one  nodes  in  the  current  section.  If  the  node  order  is  not 
equal  to  1,  control  transfers  to  statement  200.  Otherwise,  the  line  number  of 
the  node  is  obtained  from  the  NODPTR  array.  The  packed  segment  numbers  of  the 
neighboring  segments  are  stored  in  variable  NBRS.  The  number  of  neighboring 
segments  is  stored  in  variable  LIM.  If  LIM  is  greater  than  1,  control  trans- 
fers to  statement  40.  Otherwise,  the  street  segment  is  a dead  end. 

At  statement  10,  the  number  of  times  the  segment  is  to  be  traversed 
is  incremented  by  1.  The  number  of  the  node  at  the  end  of  the  segment  opposite 
the  order-one  node  is  obtained  and  stored  in  variable  NODE.  The  loop  through 
statement  20  finds  the  line  number  of  this  node,  JJ,  in  the  array  of  nodes  for 
this  section.  At  statement  30  the  order  of  the  node  is  incremented  by  1.  The 
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order  of  the  order-one  node  is  set  equal  to  2.  Control  transfers  to  state- 
ment 200. 


At  statement  40,  variable  MAXLVL,  the  maximum  number  of  steps  allowed 
in  a return  path,  is  set  to  4.  The  loop  through  statement  50  examines  the 
neighbors  of  the  order-one  node  to  obtain  the  segment  that  lies  in  the  current 
section.  The  end-point  numbers  and  the  one-way  indicator  for  the  segment  are 
saved. 

The  loop  through  statement  130  seeks  paths,  first  from  the  section 
to  the  order-one  node  and  then  from  the  node  to  the  section.  If  travel  on  the 
path  opposes  the  direction  of  travel  on  the  segment  already  in  the  section, 
control  transfers  to  statement  130.  Otherwise,  the  data  for  the  segment  in  the 
ISTG  array  are  adjusted  so  that  the  direction  of  travel  on  the  segment  opposes 
the  direction  of  travel  selected  for  the  path.  The  loop  through  statement  75 
clears  the  arrays  used  by  subroutine  TREE  to  build  paths  to  or  from  a node. 
Subroutine  TREE  is  then  called  to  build  paths,  first  to  and  then  from  the 
order-one  node.  The  paths  are  limited  to  MAXLVL  steps. 


The  loop  through  statement  120  examines  each  node  in  the  section  to 
find  paths  to  the  order-one  node.  If  a node  from  the  section  is  not  on  the 
path,  control  transfers  to  statement  120.  Otherwise,  the  next  segment  in  the 
path  is  examined.  If  the  segment  is  not  in  the  section,  control  transfers  to 
statement  100.  Otherwise,  the  line  number  of  the  next  node  in  the  path,  L, 
and  the  number  of  the  next  segment  in  the  path,  M,  are  obtained.  If  the  seg- 
ment number  is  positive,  control  returns  to  statement  80.  Otherwise,  an  error 
exists  in  the  subroutine  and  an  error  message  is  printed.  Program  execution 
terminates  with  a call  to  subroutine  DUMP. 

At  statement  100,  the  time  of  travel  from  the  node  at  line  L to  the 
order-one  node  is  compared  to  the  current  value  of  TMIN.  If  TMIN  is  less  than 
or  equal  to  the  travel  time,  control  transfers  to  statement  120.  Otherwise, 
TMIN  Is  set  equal  to  the  travel  time.  The  number  of  segments  in  the  path  is 
set  equal  to  0.  The  loop  through  statement  110  saves  the  path  from  the  node 
back  to  the  order-one  node.  The  segment  numbers  are  stored  in  array  IPATHS, 
and  the  node  numbers  are  stored  in  array  IPATHN.  If  a segment  that  is  already 
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in  the  section  is  encountered,  control  returns  to  statement  90.  Otherwise, 
the  segments  in  the  path  are  counted  and  the  count  is  stored  in  variable  NPSV. 
If  the  order-one  node  is  encountered,  control  transfers  to  statement  120. 
Otherwise,  the  loop  continues  until  MAXLVL  steps  have  been  saved. 

Statement  120  marks  the  end  of  the  loop  that  examines  the  nodes  in 
the  section  as  starting  points  for  the  path  to  the  order-one  node.  Statement 
130  marks  the  end  of  the  loop  that  examines  the  two  types  of  paths:  paths 
traveled  away  from  the  node  and  paths  traveled  toward  the  node.  Following 
statement  130,  if  a valid  path  has  been  found,  control  transfers  to  statement 
140.  Otherwise,  MAXLVL  is  doubled.  If  MAXLVL  is  less  than  or  equal  to  32, 
control  returns  to  statement  70.  Otherwise,  program  execution  terminates  on  a 
ST0P1401  statement. 

% 

At  statement  140,  the  segment  connecting  the  order-one  node  to  the 
section  is  restored.  NPSV  is  set  equal  to  1,  and  the  segment  number  is  saved 
as  the  return  path  to  the  section.  Variable  NBRS  is  set  equal  to  the  segment 
number.  Control  returns  to  statement  10. 

At  statement  150,  the  order  of  the  order-one  node  is  changed  to  2. 

The  loop  through  statement  160  finds  the  node  in  the  section  at  which  the  path 
begins.  The  order  of  that  node  is  incremented  by  1.  The  loop  through  state- 
ment 180  transfers  the  segments  in  the  path  from  array  IPATHS  to  array  IPSSV. 

In  addition,  the  segments  are  assigned  to  the  current  section,  and  the  number 
of  times  for  traversal  is  set  to  1 for  each  segment.  The  count  of  nodes  in  the 
section,  KNX,  is  incremented  by  1 for  each  node  in  the  path.  The  order  of  each 
of  these  nodes  is  set  to  2. 

Statement  200  marks  the  end  of  the  loop  that  searches  for  order-one 
nodes  in  the  section.  Control  then  returns  to  the  calling  program. 

s.  Subroutine  GENDM 

Subroutine  GENDM  generates  a matrix  of  distances  from  each  odd  node 
in  a section  to  each  other  odd  node.  The  subroutine  has  four  arguments.  The 
first,  ONDMTX,  is  a double-dimensioned  array  used  for  the  distances  from  each 
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node  to  each  other  node.  The  second,  NODD,  is  the  number  of  nodes.  The  third, 
NODPTR,  is  an  array  of  pointers  to  the  nodes.  The  fourth,  NTRIP , is  the  sec- 
tion number. 


The  loop  through  statement  30  selects  each  node  to  be  the  starting 
point  of  a distance  tree.  The  node  number  is  stored  in  variable  NORG.  The 
loop  through  statement  10  clears  the  arrays  used  by  subroutine  TREE.  Subrou- 
tine TREE  is  called  to  build  a distance  tree  from  node  NORG,  using  only  seg- 
ments in  the  current  section.  The  loop  on  statement  20  transfers  the  distances 
from  array  DIST  to  matrix  ONDMTX.  When  the  distance  matrix  is  complete,  con- 
trol returns  to  the  calling  program. 


t.  Subroutine  SELORD 

* 

* 

Subroutine  SELORD  pairs  odd  nodes  in  a manner  that  comes  close  to 
minimizing  the  total  distance  between  pairs.  SELORD  has  five  arguments.  The 
first,  NNN,  is  the  dimension  of  the  distance  matrix.  The  second,  NODD,  is  the 
number  of  odd  nodes.  The  third,  DIST,  is  a matrix  of  distances  from  each  node 
to  each  other  node.  The  fourth,  IPART,  is  an  array  of  sequence  numbers  for  the 
partners  of  the  nodes.  This  array  is  generated  by  the  subroutine.  The  fifth 
argument,  IDIR,  selects  the  method  used  to  pair  the  odd  nodes.  If  IDIR  equals 
-1 , nodes  are  paired  starting  with  the  pair  having  the  largest  near-neighbor 
distance.  If  IDIR  equals  +1,  nodes  are  paired  starting  with  those  having  the 
smallest  near-neighbor  distance. 


The  loop  through  statement  20  clears  arrays  MINJ  and  IPART.  Array 
MIND  is  set  to  1,000,000  in  this  loop. 

The  loop  through  statement  60  selects  each  odd  node  as  a part  of  a 
minimum-distance  pair.  The  loop  through  statement  40  examines  the  remaining  odd 
nodes  to  complete  the  minimum-distance  pairing.  If  the  loop  indices  are  equal, 
control  transfers  to  statement  40.  Otherwise,  the  distance  between  the  nodes 
with  sequence  numbers  I and  J is  retrieved  and  stored  in  variable  D.  If  D is 
greater  than  or  equal  to  the  current  minimum  distance  in  the  MIND  array,  con- 
trol transfers  to  statement  40.  Otherwise,  the  sequence  number,  J,  is  saved 
in  the  MINJ  array  and  the  distance  is  saved  in  the  MIND  array.  When  the  loop 
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through  statement  60  Is  complete,  the  sequence  number  of  each  node's  nearest 
neighbor  will  be  in  array  MINJ.  Note  that  some  nodes  may  be  near  neighbors  t.o 
more  than  one  node. 

The  pairing  of  odd  nodes  begins  with  the  largest  near-neighbor 
distance  if  IDIR  equals  -1,  or  with  the  smallest  near-neighbor  distance  if  1D1R 
equals  +1.  The  loop  through  statement  200  produces  an  odd-node  pair  with  each 
execution  of  the  loop.  Variable  IEXTI  is  set  to  0.  Variable  IEXTD  is  set  to 
either  0 or  2,000,000,  depending  on  the  value  of  IDIR.  The  loop  through  state- 
ment 80  examines  each  node.  If  a node  already  has  a partner,  or  if  the  dis- 
tance to  its  nearest  neighbor  makes  these  two  nodes  a less  desirable  choice  for 
pairing  than  some  other  pair,  control  transfers  to  statement  80.  Otherwise, 
the  sequence  numbers  of  the  node  and  its  nearest  neighbor  are  saved  in  variables 

IEXTI  and  IEXTD.  If  no  suitable  pair  has  been  found,  control  transfers  to 

statement  200.  Otherwise,  the  sequence  numbers  of  the  nodes  are  store  ! in  the 
IPART  array,  making  each  node  the  partner  of  the  other  node. 

The  formation  of  an  odd-node  pair  may  cause  one  of  the  two  nodes  in 
a minimum-distance  pair  to  be  used.  In  this  case,  a new  near  neighbor  must  be 

found  for  the  remaining  node.  The  loop  through  statement  160  searches  for 

nodes  to  which  new  near  neighbors  must  be  assigned.  If  the  node  indicated  by 
the  loop  index  has  a partner  specified  in  the  IPART  array,  control  transfers 
to  statement  160.  Otherwise,  if  the  neighbor  of  the  node  does  not  have  a part- 
ner specified  in  the  IPART  array,  control  transfers  to  statement  160.  Other- 
wise, the  near-neighbor  distance  for  the  node  is  set  temporarily  to  1,000,000. 
The  loop  through  statement  140  selects  a near  neighbor  for  the  node  from  those 
nodes  not  assigned  partners  in  the  IPART  array. 

Statement  160  marks  the  end  of  the  loop  that  selects  new  minimum- 
distance  pairs.  Statement  200  marks  the  end  of  the  loop  that  generates  odd- 
node  pairs.  Control  returns  to  the  calling  program. 

u.  Subroutine  PATH 

Subroutine  PATH  takes  the  packed  representation  of  two  paths  meeting 
at  a node  and  produces  both  a packed  and  an  unpacked  version  of  the  path  from 
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the  first  to  the  last  node.  The  subroutine  has  six  arguments.  The  first, 
ISTPO,  Is  a three-dimensional  array  of  packed  paths.  The  second,  1STPN,  Is 
a three-dimensional  array  of  packed  paths  produced  by  subroutine  PATH.  The 
third  argument,  II,  Is  the  sequence  number  of  the  node  that  starts  the  first 
path  to  be  examined  by  the  subroutine.  The  fourtn  argument,  KK,  is  the  se- 
quence number  of  the  node  that  ends  the  first  path  and  starts  the  second  path. 
The  fifth  argument,  JJ,  Is  the  sequence  number  of  the  node  that  ends  the  sec- 
ond path.  The  sixth  argument,  NG,  Is  an  error  indicator.  NG  Is  set  to  -1 
when  the  new  packed  path  exceeds  the  available  storage,  or  to  0 when  no  prob- 
lems are  found  in  the  combined  path.  NG  Is  set  to  1 if  a node  in  the  first 
path  occurs  a second  time  in  the  second  path. 

Variable  IDIM,  the  number  of  words  available  for  each  packed  path, 

Is  set  to  2.  Variable  JDIM,  the  maximum  number  of  nodes  used  in  the  paths,  is 
evaluated  using  the  locations  of  the  JSTP  and  l ST P arrays.  The  loop  on  state- 
ment 10  clears  array  1STP,  which  Is  used  to  hold  the  node  numbers  in  the  un- 
packed path.  Variable  NT,  the  number  of  nodes  in  the  path,  is  set  to  0.  Vari- 
able I is  set  equal  to  the  sequence  number  of  the  starting  node  of  the  first 
path,  and  variable  J Is  set  equal  to  the  sequence  number  of  the  node  ending 
the  first  path. 

The  loop  through  statement  100  unpacks  each  of  the  two  packed  path 
words.  Variable  MID  is  set  equal  to  one  more  than  the  number  of  nodes  stored 
In  the  ISTP  array.  The  nodes  corresponding  to  variable  1 and  the  paired  part- 
ner to  that  node  are  stored  In  consecutive  locations  in  array  ISTP.  The  count 
of  nodes,  NE,  is  incremnted  by  2.  Variable  NW,  a pointer  to  the  appropriate 
packed  path  word,  is  set  equal  to  1 . At  statement  20  the  packed  path  word  is 
moved  from  matrix  1S1P0  to  variable  IWD.  At  statement  40,  if  IWP  is  positive, 
control  transfers  to  statement  60.  Otherwise,  NW  is  Incremented  by  1 . If  NW 
is  less  than  or  equal  to  IDIM,  control  transfers  to  statement  20.  Otherwise, 
control  transfers  to  statement  80. 

The  two  statements  starting  at  statement  60  retrieve  the  next  node, 
M,  from  the  packed  path.  Node  M and  its  partner  are  stored  in  the  next  two 
available  locations  of  array  ISTP.  The  count  of  nodes  Is  incremented  by  2. 

The  remainder  of  the  packed  path  is  transferred  from  variable  N1WD  to  vari- 
able IWD.  Control  transfers  back  to  statement  40. 
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At  statement  80,  pointer  J is  transferred  to  variable  I.  Variable 
J is  set  equal  to  the  sequence  number  of  the  final  node  of  the  second  path. 
Statement  100  marks  the  end  of  the  loop  that  examines  the  two  packed  path  words. 

Following  statement  100,  the  count  of  nodes  is  incremented  by  1,  and 
JJ  is  stored  as  the  last  entry  in  the  ISTP  array.  Variable  MIDM  is  set  equal 

to  one  less  than  the  number  of  nodes  in  the  first  path,  and  MIDP  is  set  equal 

to  one  more  than  the  number  of  nodes  in  the  first  path. 

Error  indicator  NG  is  set  to  1.  The  double  loop  through  statement 
120  searches  the  second  path  for  a node,  other  than  the  starting  or  stopping 
node,  that  has  already  occurred  in  the  first  path.  If  a node  occurs  in  both 
paths,  control  returns  to  the  calling  program.  Otherwise,  following  statement 
120,  the  error  indicator  is  set  to  0. 

The  loop  through  statement  160  combines  the  two  packed  paths  into  one 
packed  path.  The  loop  is  executed  once  for  each  word  in  the  packed  path. 
Pointers  to  the  starting  and  stopping  nodes  in  a packed  path  word  are  computed 
and  stored  in  variables  KI  and  KF.  The  packed  path  word  in  the  ISTPN  array 
is  set  to  0.  Variable  IWD  is  set  to  0.  If  the  final  subscript  is  less  than 
the  initial  subscript,  control  transfers  to  statement  160.  Otherwise,  the 
loop  through  statement  140  builds  the  packed  path  word,  accumulating  the  result 
in  variable  IWD.  At  statement  160,  variable  IWD  is  transferred  to  the  appro- 
priate location  in  the  ISTPN  array.  Following  statement  160,  if  the  number  of 
nodes  in  the  packed  path  causes  the  available  storage  for  packed  paths  to  be 
exceeded,  variable  NG  is  set  to  -1.  Control  returns  to  the  calling  program. 

v.  Subroutine  NEXTM 

Subroutine  NEXTM  doubles  the  maximum  number  of  steps  in  a path  while 
seeking  more  profitable  paths  from  one  node  to  another.  The  subroutine  has  six 
arguments.  The  first,  MTXO,  is  a matrix  of  profits  in  the  paths  between  nodes. 
The  second,  MTXN,  is  the  new  matrix  of  profits  produced  by  the  subroutine  when 
the  maximum  number  of  steps  in  each  path  is  allowed  to  double.  The  third  argu- 
ment, N0DD,  is  the  number  of  nodes.  The  fourth  and  fifth  arguments,  ISTP0  and 
ISTPN,  are  matrices  of  packed  paths  for  the  original  paths  and  the  new  paths 
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produced  by  the  subroutine.  The  sixth  argument,  NRES,  indicates  whether  or 
not  a profitable  cycle  has  been  found. 

Variable  NRES,  the  profitable-cycle  indicator,  is  set  to  0.  The  num- 
ber of  nodes  is  stored  in  variable  NEDGE.  Variable  MAXD  is  set  to  0. 

The  loop  through  statement  160  selects  each  node  in  turn  as  the  start- 
ing and  stopping  node  of  a cycle.  The  next  four  statements  transfer  the  profit 
and  path  for  the  cycle  from  the  old  matrices  to  the  new  matrices. 

The  loop  through  statement  140  examines  each  node  as  the  midpoint  of 
a cycle.  If  the  midpoint  node  is  equal  to  the  starting  node,  control  transfers 
to  statement  140.  Otherwise,  the  profit  of  the  path  from  the  starting  node  to 
the  center  and  then  back  to  the  starting  node  is  computed  and  stored  in  vari- 
able MPR.  If  MPR  is  less  than  or  equal  to  the  current  profit  of  the  cycle, 
control  transfers  to  statement  140.  Otherwise,  subroutine  PATH  is  called  to 
produce  the  packed  path  for  the  new  cycle.  If  error  indicator  NG  is  nonzero, 
control  transfers  to  statement  140.  Otherwise,  the  new  profit  is  stored  in  the 
profit  matrix.  If  the  new  profit  is  less  than  or  equal  to  MAXD,  control  trans- 
fers to  statement  140,  Otherwise,  MAXD  is  set  equal  to  the  new  profit.  NES  is 
set  equal  to  the  number  of  steps  in  the  path.  The  loop  on  statement  120  trans- 
fers the  nodes  in  the  current  profitable  cycle  from  the  ISTP  array  to  the  JS TP 
array.  Statement  140  marks  the  end  of  the  loop  that  selects  the  midpoint  node 
of  each  cycle.  Statement  160  marks  the  end  of  the  loop  on  the  nodes. 

If  MAXD  is  essentially  zero  (it  is  compared  to  a small  number  to  take 
roundoff  errors  into  account),  control  transfers  to  statement  180.  Otherwise, 
the  profitable-cycle  indicator  is  set  equal  to  1 and  control  returns  to  the 
calling  program. 

Following  statement  180,  a double  loop  through  statement  240  is  be- 
gun. The  loop  seeks  more  profitable  paths  from  one  node  to  another.  If  the 
two  nodes  are  the  same,  control  transfers  to  statement  240.  Otherwise,  MAXP  is 
set  equal  to  the  current  profit  in  the  path  from  node  I to  node  J.  The  two 
words  in  the  packed  path  are  transferred  from  matrix  ISTP0  to  matrix  ISTPN.  The 
profit  is  transferred  to  matrix  MTXN.  The  loop  through  statement  120  selects 
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each  node  to  be  the  midpoint  of  a new  path  from  node  I to  node  J.  If  the  node 
selected  is  either  the  starting  or  ending  point,  control  transfers  to  state- 
ment 220.  Otherwise,  variable  MPR  is  set  equal  to  the  profit  of  the  new  path. 

If  MPR  is  less  than  or  equal  to  MAXP,  control  transfers  to  statement  220. 
Otherwise,  subroutine  PATH  is  called  to  produce  the  new  packed  path  between 
the  nodes.  If  the  error  indicator  is  nonzero,  control  transfers  to  statement 
220.  Otherwise,  the  new  profit  is  stored  in  both  matrix  MTXN  and  in  variable 
MAXP.  Statement  220  marks  the  end  of  the  loop  on  the  midpoint  nodes.  State- 
ment  240  marks  the  end  of  the  double  loop  on  the  nodes.  Control  returns  to 
the  calling  program. 

w.  Subroutine  SOLV 

Subroutine  SOLV  improves  the  odd-node  pairing  until  the  mini  ium  total 
pairing  distance  is  obtained.  The  subroutine  has  three  arguments.  The  first, 
DIST,  is  an  array  of  distances  between  nodes.  The  second,  NODD,  is  the  number 
of  odd  nodes.  The  third,  IERR,  is  an  error  indicator.  The  subroutine  examines 
all  closed  paths  of  a given  length  constructed  from  odd-node  pairs.  If  the 
cycle  produces  a shorter  total  pairing  distance  when  the  end  node  of  each  pair 
is  paired  with  the  first  node  of  the  next  pair,  the  new  pairing  is  used.  The 
algorithm  then  repeats,  starting  with  two-step  paths.  If  no  profitable  cycle 
is  found,  the  algorithm  doubles  the  number  of  steps  allowed  in  the  path.  If 
no  profitable  cycle  has  been  found  when  all  odd-nodes  pairs  have  been  used  in 
the  path,  the  optimum  odd-node  pairing  has  been  obtained. 

Variable  NCH,  a count  of  iterations  through  subroutine  SOLV,  is  set 
to  0.  Variable  NNPW,  the  number  of  nodes  that  can  be  packed  into  a word,  is 
computed.  Variable  NN  is  set  equal  to  2.  The  loop  through  statement  40  doubles 
NN  until  it  exceeds  the  number  of  odd  nodes.  When  NN  exceeds  the  number  of  odd 
nodes,  control  transfers  to  statement  60.  Variable  LIMK  has  a value  of  one 
less  than  the  loop  index.  The  double  loop  through  statement  220  generates  path 
and  profit  matrices  for  the  most  profitable  paths,  involving  two  steps,  between 
any  two  nodes.  The  profit  from  any  node  to  its  partner  is  set  equal  to  -1000. 

The  loop  through  statement  300  doubles  the  number  of  steps  in  the 
path  until  the  maximum  number  of  steps  is  obtained.  If  the  loop  index  is  even. 
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control  transfers  to  statement  240.  Otherwise,  subroutine  NEXTM  is  called  to 
double  the  number  of  steps  used  in  the  paths.  Matrices  MTX1  and  ISTP1  contain 
the  old  profits  and  paths;  matrices  MTX2  and  ISTP2  receive  the  new  profits  and 
paths.  If  the  profitable-cycle  indicator  is  nonzero,  indicating  a profitable 
cycle,  control  transfers  to  statement  400.  Otherwise,  control  transfers  to 
statement  300. 


At  statement  240,  subroutine  NEXTM  is  called.  This  time,  however, 
MTX2  and  ISTP2  are  the  old  matrices,  and  MTXl  and  ISTP1  are  the  new  matrices. 

If  a profitable  cycle  is  found,  control  transfers  to  statement  400.  Otherwise, 
the  loop  continues,  following  the  loop,  control  returns  to  the  calling  program. 

Following  statement  400,  the  iteration  count  is  incremented  by  1. 

If  fewer  than  20  iterations  have  been  performed,  control  transfers  to  state- 
ment 440.  Otherwise,  an  error  message  is  printed.  The  error  indicator,  IERR, 
is  set  to  1,  and  control  returns  to  the  calling  program. 

At  statement  440,  a loop  through  statement  460  is  begun.  This  loop 
selects  an  alternative  pairing  ot  odd  nodes  from  the  nodes  in  the  profitable 
cycle.  Control  returns  to  statement  120. 

x.  Subroutine  TRAVEL 

Subroutine  TRAVEL  finds  a path  that  travels  all  of  the  segments  in 
a section,  starting  and  stopping  at  specified  nodes.  It  has  12  arguments.  The 
first  and  second,  1ST  and  ISP,  are  the  starting  and  stopping  node  numbers.  The 
third,  NSEG,  is  the  number  of  segments.  The  fourth  and  fifth  arguments,  NN1 
and  NN2,  are  arrays  of  the  starting  and  ending  nodes  of  the  segments.  The 
sixth  argument,  NS6,  is  an  array  of  segment  numbers.  The  seventh,  NWAY,  is  an 
array  of  one-way  indicators.  The  eighth,  NTIMES,  is  an  array  that  specifies 
the  numbo”-  of  times  each  segment  must  be  traversed.  The  ninth  argument,  NN, 
is  the  number  of  segments  in  the  final  path.  The  tenth  and  eleventh  argu- 
ments, IPh  and  IPS,  are  arrays  for  the  node  and  segment  numbers  in  the  path. 
These  arrays  arc  generated  by  the  subroutine.  The  twelfth  argument,  IERR,  is 
an  error  indicator. 
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Error  indicator  I ERR  and  pointer  IPT1  are  set  initially  to  1.  Vari- 
able M is  set  equal  to  KSDIM,  the  dimension  of  the  arrays  used  in  bidding  the 
path.  Pointer  IPT2  is  set  equal  to  the  number  of  segments. 

At  statement  10,  if  IPT2  is  less  than  IPT1,  control  transfers  to 
statement  40.  Otherwise,  if  the  segment  indicated  by  IPT1  is  to  be  traversed 
only  once,  control  transfers  to  statement  15.  Otherwise,  pointer  IPT1  is  in- 
cremented by  1.  Control  transfers  back  to  statement  10.  At  statement  15,  if 
the  segment  indicated  by  IPT2  is  to  be  traversed  twice,  control  transfers  to 
statement  30.  Otherwise,  pointer  IPT2  is  decremented  by  1.  Control  transfers 
to  statement  10. 

The  15  statements  starting  at  statement  30  interchange  the  data  for 
the  segments  indicated  by  pointers  IPT1  and  IPT2.  The  interchange  put  all 
segments  that  are  to  be  traversed  twice  at  the  front  of  the  arrays.  F inter 
IPT1  is  incremented  by  1,  and  pointer  IPT2  is  decremented  by  1 following  the 
interchange.  Control  transfers  back  to  statement  10. 

At  statement  40,  a loop  on  statement  55  is  begun.  This  loop  trans- 
fers the  number  of  times  the  segments  are  to  be  traversed  from  the  NTIMES 
array  to  the  KTIMES  array.  The  loop  on  statement  205  clears  the  KS  array. 

This  array  counts  the  number  of  times  the  subroutine  backs  up  to  each  step  in 
the  path  during  the  trial-and-error  search  for  the  path. 

The  starting  node  number  is  stored  as  the  first  entry  in  the  IPN 
array.  The  number  of  steps  currently  in  the  path,  NN,  is  set  to  1.  The  first 
entries  in  the  KLINE  and  NEXTII  arrays  are  set  to  0.  The  KLINE  array  indica- 
tes the  line  number  in  the  segment  data  from  which  a particular  step  in  the 
path  comes.  The  NEXTII  array  gives  the  line  number  in  the  segment  data  at 
which  a search  will  begin  for  a different  segment  as  a particular  step  in  the 
path.  At  statement  65,  II  is  set  initially  to  1. 

At  statement  90,  a loop  begins  through  statement  60.  This  loop  will 
examine  segments,  starting  at  line  II  of  the  segment  data,  to  find  a segment 
that  can  continue  the  current  path.  If  the  segment  indicated  by  the  loop  index 
cannot  be  used  again,  or  if  its  use  would  create  a U-turn,  control  transfers 
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to  statement  60.  Otherwise,  if  the  starting  node  of  the  segment  is  the  same 
as  the  last  node  currently  in  the  path,  control  transfers  to  statement  70.  If 
not,  and  if  the  segment  is  a one-way  segment,  or  if  the  ending  node  of  the  seg- 
ment is  not  the  last  node  currently  in  the  path,  control  transfers  to  statement 
60.  Otherwise,  the  current  segment  is  added  to  the  path.  Its  number  is  stored 
in  array  IPS,  and  its  starting  node  is  stored  in  array  IPN.  The  count  of  steps 
in  the  path  is  incremented  by  1.  The  line  number  of  the  segment  is  stored  in 
the  location  in  array  KLINE  corresponding  to  the  current  step  in  the  path.  The 
location  of  the  NEXTII  array  corresponding  to  the  current  step  in  the  path  is 
set  equal  to  the  current  loop  index  plus  1.  The  availability  count  for  the  seg- 
ment is  decremented  by  1.  Control  returns  to  statement  65. 

At  statement  70  the  segment  is  added  to  the  path.  Its  number  is 
stored  in  the  IPS  array.  The  node  count  is  incremented  by  1.  The  ending  node 
number  is  stored  in  the  IPN  array.  The  line  number  of  the  current  segment  arid 
the  loop  index  plus  1 are  stored  in  the  appropriate  locations  of  arrays  KLINE 
and  NEXTII.  The  availability  count  for  the  segment  is  decremented  by  1.  Con- 
trol returns  to  statement  65.  Statement  60  marks  the  end  of  the  loop  that 
searches  the  segment  data  for  a suitable  segment  to  be  added  to  the  path. 

When  control  passes  beyond  the  loop  through  statement  60,  no  suitable 
segment  has  been  found  for  addition  to  the  path.  The  line  number  of  the  last 
segment  currently  in  the  path  is  stored  in  variable  I.  If  this  segment  is  not 
available  to  be  used  again,  control  transfers  to  statement  75.  Otherwise,  the 
segment  is  added  as  the  next  step  in  the  path,  creating  a U-turn.  The  node 
count  is  incremented  by  1,  and  the  appropriate  node  is  appended  to  the  IPN 
array.  The  line  number  of  the  segment  is  stored  in  the  appropriate  place  in 
the  KLINE  array.  The  count  of  segments  plus  1 is  stored  in  the  NEXTII  array, 
indicating  that  no  segments  remain  as  possibilities  for  this  step.  The  avail- 
ability count  for  the  segment  is  decremented  by  1,  and  control  transfers  to 
statement  65. 

At  statement  75,  if  the  last  node  in  the  path  is  not  the  appropriate 
stopping  node,  control  transfers  to  statement  85.  Otherwise,  the  loop  through 
statement  80  examines  the  availability  of  each  segment.  If  each  segment  has 
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been  used  the  appropriate  number  of  times,  control  transfers  to  statement  80. 
Otherwise,  the  line  number  of  the  segment  is  stored  in  the  IPN  array,  starting 
at  the  end  of  the  array.  A pointer,  M,  is  decremented  by  1.  Following  state- 
ment 80,  if  M does  not  point  to  the  end  of  the  array,  control  transfers  to 
statement  185.  Otherwise,  the  error  indicator,  IERR,  is  cleared.  Following 
statement  201,  if  an  error  condition  is  present,  an  error  message  is  printed. 
Control  returns  to  the  calling  program. 

At  statement  85,  variable  I is  set  equal  to  the  line  number  of  the 
last  segment  in  the  path.  If  this  line  number  is  less  than  or  equal  to  0, 
control  transfers  to  statement  201.  Otherwise,  the  availability  count  of  the 
segment  is  incremented  by  1.  Variable  II  is  set  equal  to  the  value  of  NEXTII 
corresponding  to  the  current  step  in  the  path.  The  location  in  the  KS  array 
corresponding  to  the  current  step  's  incremented  by  1.  If  the  subroutine  has 
backed  up  to  this  step  in  the  path  more  than  MAXKS  times,  control  transfers 
to  statement  201.  Otherwise,  the  -er  of  steps  in  the  path  is  decreased  by 
1.  If  any  segments  remain  as  possiole  replacements  at  this  step  in  the  path, 
control  transfers  to  statement  90.  Otherwise,  control  transfers  to  statement 
85. 


At  statement  185,  pointer  M is  incremented  by  1 . It  now  points  to 
the  first  line  number  of  an  unused  segment  in  the  IPN  array.  At  statement  186, 
the  line  number  of  the  last  step  in  the  path  is  stored  in  variable  I.  If  the 
line  number  is  less  than  or  equal  to  0,  control  transfers  to  statement  201. 
Otherwise,  the  availability  count  of  the  segment  is  incremented  by  1.  Vari- 
able II  is  set  from  the  NEXTII  array.  The  appropriate  location  in  the  KS  ar- 
ray is  Incremented  by  1.  If  this  value  exceeds  MAXKS,  control  transfers  to 
statement  201 . 

At  statement  210,  NN  is  decremented  by  1.  The  loop  through  state- 
ment 83  examines  the  unassigned  segments  to  determine  whether  one  of  them  can 
be  used  as  the  next  step  in  the  path.  If  a suitable  segment  is  found,  control 
transfers  to  statement  84.  Otherwise,  control  transfers  to  statement  186. 

At  statement  84,  M is  restored  to  its  maximum  value,  KSDIM.  If 
variable  II  Indicates  that  segments  remain  to  be  examined,  control  transfers 
to  statement  90.  Otherwise,  control  transfers  to  statement  85. 
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y.  Subroutine  OPTPATH 


Subroutine  OPTPATH  performs  two  types  of  optimization  on  the  traver- 
sal path  produced  by  subroutine  TRAVEL:  first,  it  finds  the  shortest  path  from 
the  start  to  the  end  of  each  travel -only  sequence  of  segments;  and  second,  it 
removes  sequences  of  travel -only  segments  that  are  traveled  twice  in  the  same 
direction.  The  subroutine  has  five  arguments.  The  first  and  second,  NNP  and 
NSP,  are  arrays  of  the  nodes  and  segments  in  the  travel  path.  The  third,  TRTY , 
is  an  array  that  indicates  for  each  segment  whether  both  travel  and  collection 
or  travel,  only,  will  be  required.  This  array  is  generated  by  the  subroutine. 
The  fourth  argument,  N,  is  the  number  of  nodes  in  the  path.  The  fifth  argu- 
ment, ISC,  is  the  section  number. 

Variable  NM1  is  set  equal  to  N - 1.  The  loop  through  statement  20 
sets  the  travel  type  for  each  of  the  NM1  segments.  The  travel  type  is  ini- 
tially set  to  T to  indicate  travel  only.  If  collection  is  not  required  on 
this  segment,  control  transfers  to  statement  20.  Otherwise,  the  travel  type 
is  changed  to  a C to  indicate  that  collection  is  required.  If  the  segment  is 
traveled  more  than  once,  the  number  of  times  of  traversal  is  set  to  0.  This 
prevents  the  segment  from  being  marked  for  collection  twice.  Statement  20 
marks  the  end  of  the  loop  that  sets  the  travel  type. 

The  loop  through  statement  30  seeks  the  first  segment  on  which  col- 
lection is  required.  When  this  segment  is  found,  control  transfers  out  of  the 
loop  to  statement  40.  At  statement  40,  the  segments,  nodes,  and  travel  types 
in  the  pat’.,  are  moved  forward  so  that  the  first  collection  segment  will  be  the 
first  segment  in  the  path.  The  number  of  nodes  in  the  path  description  is  re- 
computed. 


Variable  IF  is  set  initially  to  1.  At  statement  60,  a loop  begins 
through  statement  70.  This  loop  examines  each  segment,  seeking  a segment  on 
which  only  travel  is  required.  If  such  a segment  is  found,  control  transfers 
to  statement  80.  Otherwise,  control  transfers  to  statement  200. 

At  statement  80,  a loop  through  statement  90  is  begun.  This  loop 
seeks  the  next  segment  requiring  collection.  If  such  a segment  is  found. 
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control  transfers  to  statement  100.  Otherwise,  control  transfers  to  state- 
ment 200. 


At  statement  100,  the  number  of  consecutive  traveled  segments,  IS, 
is  computed.  If  IS  is  equal  to  1,  control  transfers  to  statement  60.  Other- 
wise, variable  NS  is  set  to  0.  If  the  traveled  segments  start  and  end  at  the 
same  node,  control  transfers  to  statement  115.  Otherwise,  the  loop  through 
statement  110  clears  the  arrays  used  by  subroutine  TREE.  Subroutine  TREE  is 
called  to  find  the  shortest  (in  time)  paths  starting  at  the  node  that  ends 
the  travel  stretch.  Subroutine  TRACE  is  called  to  obtain  the  node  numbers, 
segment  numbers,  and  number  of  segments  in  the  shortest  path  from  the  start- 
ing to  the  stopping  node. 

If  the  path  neither  starts  nor  ends  with  a U-turn,  control  transfers 
to  statement  115.  Otherwise,  function  CUMDIS  is  used  to  find  the  travel  dis- 
tances of  the  original  and  the  new  travel  paths.  If  the  new  path  saves  less 
than  half  a mile  and  the  old  path  is  less  than  three  times  the  length  of  the 
new  path,  control  returns  to  statement  60.  Otherwise,  at  statement  115  vari- 
able IDEL  is  set  equal  to  the  number  of  segments  by  which  the  new  path  is 
longer  than  the  old  path.  The  number  of  nodes  in  the  trip  using  the  new  path 
is  computed  and  stored  in  variable  NEWN.  If  IDEL  is  negative,  control  trans- 
fers to  statement  140;  if  it  is  zero,  control  transfers  to  statement  160;  and 
if  it  is  positive,  control  transfers  to  statement  120. 

Statement  120  and  the  loop  through  statement  130  move  path  informa- 
tion in  arrays  NNP,  NSP,  and  TRTY  to  make  room  for  the  additional  segments  in 
the  new  path.  Control  then  transfers  to  statement  160. 

Starting  at  statement  140,  the  loop  through  statement  150  moves  the 
data  in  the  travel  path  forward  by  the  number  of  segments  by  which  the  new 
path  is  shorter  than  the  old  path. 

At  statement  160,  if  there  are  no  steps  in  the  new  path,  control 
transfers  to  statement  180.  Otherwise,  the  loop  through  statement  170  moves 
the  new  travel  path  on  top  of  the  old  travel  path  in  arrays  NNP,  NSP,  and 
TRTY.  At  statement  180,  variable  IF,  the  number  of  the  last  step  in  the  path 
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examined  thus  far,  is  recomputed  to  take  into  account  the  change  in  the  path 
length.  Variable  N is  set  equal  to  variable  NEWN.  Variable  NM1  is  recomputed. 
Control  returns  to  statement  60. 

At  statement  200,  variable  DO  is  set  to  0.  Variable  NM1  is  set  equal 
to  N - 1,  and  NM3  is  set  equal  to  N - 3.  The  loop  through  statement  260  se- 
lects, in  turn,  each  segment  in  the  travel  path  except  the  last  two.  If  the 
segment's  travel  type  indicates  collection,  control  transfers  to  statement  260. 
Otherwise,  variable  II  is  set  equal  to  the  loop  index,  variable  JI  is  set  equal 
to  the  loop  index  plus  1,  and  variable  NS  is  set  equal  to  the  segment  number. 

The  loop  through  statement  250  examines  the  remaining  segments.  If 
the  segment  number  is  not  equal  to  NS,  control  transfers  to  statement  250.  If 
the  segment  number  equals  NS,  and  if  the  segment  requires  collection  or  is 
traveled  in  the  opposite  direction,  control  transfers  to  statement  250.  Other- 
wise, variable  L is  set  equal  to  1.  The  loop  through  statanent  210  increments 
L for  each  identical  step  in  the  two  paths  starting  at  the  two  occurrences  of 
segment  NS. 

At  statement  220,  variable  IF  is  set  equal  to  the  number  of  the  first 
step  following  the  first  of  the  two  identical  travel  stretches.  Variable  JF 
is  set  equal  to  the  number  of  the  last  step  before  the  beginning  of  the  second 
travel  stretch.  The  loop  through  statement  230  examines  the  steps  between  IF 
and  JF,  seeking  a one-way  segment.  If  a one-way  segment  is  found,  control 
transfers  to  statement  250.  Otherwise,  variable  DN  is  set  equal  to  0.  The 
loop  on  statement  240  computes  the  length  of  the  duplicated  travel  stretch. 

If  the  length  of  the  stretch,  DN,  is  less  than  or  equal  to  DO,  control  trans- 
fers to  statanent  250.  Otherwise,  DO  is  set  equal  to  DN.  Pointers  are  saved, 
marking  the  starts  of  the  two  identical  travel  stretches  and  the  starts  of  the 
segments  between  them. 

Statement  250  marks  the  end  of  the  loop  seeking  the  second  occurrence 
of  a segment.  Statement  260  marks  the  end  of  the  loop  selecting  the  first  seg- 
ment. If  no  duplicated  travel  stretch  has  been  found  in  these  loops,  control 
transfers  to  statement  310. 
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If  a duplicated  travel  stretch  has  been  found,  pointers  are  set  to 
the  first  and  last  steps  in  each  travel  stretch.  The  11  statements  starting 
at  statement  270  reverse  the  path  between  the  two  travel  stretches.  The  first 
travel  stretch  is  overwritten  during  the  reversal.  The  loop  through  statement 
280  moves  the  remainder  of  the  travel  path  forward  in  arrays  NSP,  TRTY,  and 
NNP,  so  that  it  immediately  follows  the  portion  of  the  path  that  has  been  re- 
versed. The  count  of  nodes  in  the  entire  path  is  recomputed.  Control  returns 
to  statement  200. 

At  statement  300,  the  count  of  nodes  is  decremented  by  1.  Variable 
NM1  is  recomputed.  At  statement  310,  if  the  last  segment  in  the  travel  path 
does  not  require  collection,  control  returns  to  statement  300.  Otherwise, 
control  returns  to  the  calling  program. 

z.  Program  PHASE3 

The  problem  title  is  read  from  the  first  data  card.  The  number  of 
segments,  number  of  sections,  pointers  to  the  first  and  last  segment  in  each 
section,  and  vehicle  capacities  are  read  from  file  TAPE3.  The  problem  title 
is  printed,  followed  by  the  heading  for  the  initial  segment  assignments.  The 
loop  through  statement  16  controls  the  reading  of  segments  in  each  section. 
Variables  JI  and  JF  are  set  equal  to  the  sequence  numbers  of  the  initial  and 
final  segments  in  the  section.  The  loop  through  statement  14  reads  the  seg- 
ments in  the  current  section.  The  segment  number  is  storeu  in  the  portion  of 
the  ISTG  array  whose  first  subscript  is  equal  to  1 . At  statement  14,  the  sec- 
tion number  is  stored  in  that  portion  of  the  ISTG  array  whose  first  subscript 
is  equal  to  9 (symbolically  NSECT).  At  statement  16  the  section  number  and 
the  numbers  of  the  segments  in  that  section  are  printed.  Subroutine  ISHLSRT 
is  called  to  sort  the  segment  numbers  into  increasing  order.  The  section  num- 
bers are  carried  along  during  the  sort. 

The  segment  data  are  read  from  file  TAPE1  into  the  STG  array.  Note 
that  the  section  assignments  are  not  overwritten.  The  node  data  are  read  from 
file  TAPE2.  The  loop  through  statement  19  examines  each  segment.  If  the  num- 
ber of  houses  on  the  segment  is  0,  the  section  assignment  is  set  to  0.  At 
statement  19,  negative  house  counts  are  made  positive.  The  loop  through  state- 
ment 20  accumulates  the  total  refuse  quantity  for  each  section. 


The  loop  through  statement  36  computes  the  number  of  neighboring 
segments  for  each  node.  The  computation  examines  the  size  of  the  word  con- 
taining the  packed  neighboring-segment  numbers. 

The  node  numbers  of  the  landfill  and  garage  are  read  from  the  second 
data  card.  A heading  is  printed  for  section  assignment  changes.  Subroutine 
ADJUST  is  called  to  put  segments  that  meet  at  ordar-two  nodes  into  the  same 
section,  wherever  possible.  The  section  assignment  changes  are  printed  by  sub- 
routine ADJUST.  The  load  and  capacity  of  each  vehicle  are  printed. 

At  statement  50,  the  segment  number  and  the  new  section  assignmer.: 
for  the  segment  are  read  from  the  next  data  card.  If  an  end-of-file  card  is 
encountered  during  the  read,  control  transfers  to  statement  30.  Otherwise, 
control  resumes  at  statement  70.  At  statement  70,  the  old  section  assignment 
for  the  segment  is  transferred  to  variable  NSCO.  The  new  section  assignment 
is  stored  in  the  ISTG  array.  If  the  new  section  assignment  is  less  than  or 
equal  to  0,  control  transfers  to  statement  72.  Otherwise,  the  refuse  quantity 
on  the  segment  is  computed.  The  vehicle  loads  for  the  old  and  new  sections  are 
adjusted  accordingly.  At  statement  72,  the  section  reassignment  is  printed. 
Control  returns  to  statement  50. 

Following  statement  80,  the  final  vehicle  loads  and  capacities  are 
printed.  Segment  data,  including  the  latest  section  assignments  for  the  seg- 
ments, are  also  printed. 

Subroutine  TREE  is  called  three  times.  On  the  first  call,  minimum- 
ti.,e  paths  from  the  landfill  node  to  the  other  nodes  are  generated,  with  the 
travel  direction  away  from  the  landfill.  On  the  second  call,  minimum-time 
paths  from  the  garage  node  to  each  other  node  are  generated,  with  the  direction 
of  travel  away  from  the  garage.  On  the  third  call,  minimum-time  paths  from  the 
landfill  node  to  each  other  node  are  generated,  with  the  direction  of  travel 
toward  the  landfill. 

Subroutine  TRACE  is  called  to  find  in  the  data  the  path  from  the 
garage  to  the  landfill,  where  travel  direction  will  be  away  from  the  landfill. 
Subroutine  FLIP  is  called  twice;  first,  to  reverse  the  order  of  the  segments 
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in  the  path  and  second,  to  reverse  the  order  of  the  nodes  in  the  path.  The 
distance  from  the  landfill  to  the  garage  is  evaluated  by  function  CUMDIS  and 
is  stored  in  variable  DIST.  The  path  from  the  landfill  to  the  garage  is 
written  to  file  TAPE9. 

The  data  describing  the  paths  to  and  from  the  landfill  and  garage 
are  written  to  file  TAPE7  by  a BUFFER  OUT  statement.  File  TAPE7  is  rewound. 

The  node  numbers  of  the  landfill  and  garage  are  printed,  followed 
by  the  data  describing  paths  to  and  from  these  nodes. 

Carriage-control  variable  CC  is  set  equal  to  the  character  1.  The 
loop  through  statement  180  examines  each  of  the  nodes.  If  the  node  is  connec- 
ted to  the  landfill,  control  transfers  to  statement  140.  Otherwise,  ^n  error 
message  is  printed.  The  carriage  control  is  changed  to  a blank.  At  statement 
140,  if  the  node  is  connected  to  the  garage,  control  transfers  to  statement 
180.  Otherwise,  an  error  message  is  printed.  Following  the  loop  through 
statement  180,  the  carriage  control  is  examined  to  determine  whether  an  error 
has  occurred.  If  no  disconnected  nodes  have  been  found,  control  transfers  to 
statement  200.  Otherwise,  an  error  message  is  printed  and  program  execution 
terminates  with  a call  to  subroutine  EXIT. 

The  loop  through  statement  800  generates  the  routes  for  each  section. 
The  loop  on  statement  202  transfers  the  permanent  section  assignments  to  stor- 
age in  the  IST6  array  for  temporary  section  assignments.  Subroutine  CONNECT 
is  called  to  connect  separate  pieces  of  the  current  section.  In  this  subrou- 
tine, the  temporary  section  assignments  of  segments  needed  to  connect  pieces 
of  the  current  section  are  set  equal  to  the  current  section  number.  The  loop 
on  statement  210  clears  the  node  order  and  node  pointer  arrays. 

Variable  KN,  the  count  of  nodes  in  the  section,  is  initially  set  to 
0.  The  loop  through  statement  230  examines  each  segment.  If  the  temporary 
section  number  is  not  equal  to  the  current  section  number,  control  transfers 
to  statement  230.  Otherwise,  the  loop  through  statement  220  examines  the 
segment's  end-point  nodes.  If  the  line  number  of  the  node  is  already  stored 
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in  the  NODPTR  array,  control  transfers  to  statement  220.  Otherwise,  subroutine 
M0VE3  is  called  to  make  room  for  the  node  in  the  NODPTR  and  NODORD  arrays.  The 
line  number  of  the  node  is  stored  in  the  NODPTR  array.  At  statement  220,  the 
order  of  the  node  is  incremented  by  1.  Statement  230  marks  the  end  of  the  loop 
that  examines  each  segment. 

Subroutine  DISCON  is  called  to  remove  from  the  section  any  dangling 
(dead-end)  segments  on  which  no  refuse  collection  is  required.  During  the  exe- 
cution of  this  subroutine,  nodes  may  be  deleted  from  the  node  order  and  pointer 
arrays.  When  control  returns  from  the  subroutine,  variable  KN  will  contain  a 
count  of  the  nodes  in  the  current  section. 

Subroutine  EPXP  is  called  to  select  good  entry  and  exit  nodes  for 
the  section.  The  entry  and  exit  nodes  are  stored  by  the  subroutine  in  array 
NTF  in  COMMON  block  TEMSTG. 

The  loop  through  statement  240  saves  the  current  values  of  the  node 
orders  and  pointers  in  arrays  N0D0RDS  and  NODPTRS.  The  current  number  of  nodes 
in  the  section  is  saved  in  variable  KNS. 

The  loop  through  statement  500  selects,  in  turn,  the  two  possible 
starting  points  for  the  trip:  irst,  the  garage  and  second,  the  landfill. 

Variable  SAVDIS  is  set  equal  to  a very  large  number.  The  loop  through  state- 
ment 470  examines  up  to  JMAX  entry-point  nodes.  The  entry-point  node  number 
is  stored  in  variable  ISTART.  If  ISTART  is  0,  control  transfers  to  statement 
480.  Otherwise,  the  loop  through  statement  460  selects  up  to  KMAX  exit  nodes 
from  the  section.  The  exit-node  number  is  stored  in  variable  IST0P.  If  I STOP 
is  0,  control  transfers  to  statement  470.  Otherwise,  the  count  of  nodes  is  re- 
stored to  variable  KN  from  variable  KNS. 

The  loop  through  statement  250  examines  each  node  in  the  section. 

The  node's  order  and  pointer  are  restored  to  arrays  NODORD  and  NODPTR.  If  the 
node  is  the  entry  node  for  the  section,  its  order  is  incremented  by  1 . At 
statement  250,  if  the  node  is  the  exit  node  from  the  section,  its  order  is  in- 
cremented by  1 . 
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The  loop  through  statement  260  examines  each  segment  in  the  map  de- 
scription. The  availability  count  (number  of  traversals)  is  initially  set  to  0. 

If  the  temporary  section  assignment  indicates  that  the  segment  is  needed  in  the 
current  section,  the  availability  count  is  set  equal  to  1.  Subroutine  CL0SE1 
is  called  to  generate  paths  from  order-one  nodes  back  to  the  section.  This 
subroutine  may  cause  additional  nodes  and  segments  to  be  added  temporarily  to 
the  section.  The  availability  counts  are  incremented  accordingly.  The  node 
count,  KN,  is  set  equal  to  the  count,  KNX,  returned  by  subroutine  CLOSE  1 . 

Subroutine  M0V0DD  is  called  to  move  the  odd-order  nodes  to  the  front 
of  the  node  number  and  pointer  arrays.  If  the  count  of  odd  nodes  is  less  than 
or  equal  to  the  maximum  (32),  control  transfers  to  statement  270.  Otherwise, 
an  error  message,  including  node  numbers  and  orders,  is  printed  to  indicate 
that  the  current  section  contains  too  many  odd  nodes.  Subroutine  PRNPCS  is 
called  to  print  the  numbers  of  the  segments  in  the  section.  Control  ansfers 
to  statement  460. 

Following  statement  270,  subroutine  GENDM  is  called  to  generate  a 
matrix  of  distances  between  pairs  of  odd  nodes.  Subroutine  SEL0RD  is  called 
to  pair  the  odd  nodes  in  a manner  that  approximately  minimizes  the  total 
pairing  distance.  Error  indicator  IERR  is  set  to  0.  If  there  are  more  than 
two  odd  nodes,  subroutine  S0LV  is  called  to  generate  the  odd-node  pairing 
having  the  minimum  pairing  distance.  If  an  error  has  been  encountered,  con- 
trol transfers  to  statement  460.  Otherwise,  execution  continues  with  the 
loop  through  statement  290,  which  examines  the  odd  nodes.  If  the  sequence 
number  of  an  odd  node  exceeds  that  of  its  partner,  control  transfers  to  state-  j 

ment  290.  Otherwise,  pointers  to  the  nodes  are  stored  in  variables  Nl  and  N2. 

The  loop  through  statement  280  clears  arrays  to  be  used  by  subroutine  TREE. 

j 

Subroutine  TREE  is  called  to  generate  minimum-distance  paths  from  one  of  the 
odd  nodes  to  the  remaining  nodes  in  the  section.  The  number  of  traversals 
for  each  segment  on  the  shortest  path  between  the  odd  nodes  is  incremented  by 
1.  Statement  290  marks  the  end  of  the  loop  that  examines  the  odd  nodes. 

The  loop  through  statement  300  examines  each  segment  in  the  map  de- 
scription. If  the  number  of  traversals  is  less  than  or  equal  to  2,  control 
transfers  to  statement  300.  Otherwise,  if  the  street  is  not  a one-way  street, 
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the  number  of  traversals  is  reduced  by  a multiple  of  2 until  it  is  either 
1 or  2. 


The  loop  on  statement  310  clears  the  IPN  and  IPS  arrays.  Variable 
JSEG  is  set  equal  to  0.  The  loop  through  statement  380  examines  each  segment. 

If  the  segment  is  not  needed  in  the  current  section,  control  transfers  to  state- 
ment 380.  Otherwise,  the  segment  number,  end-point  node  numbers,  availability 
count,  and  number  of  ways  of  travel  are  transferred  to  arrays  JNSG,  JNN1 , JNN2, 
JTIM,  and  JNWY . Subroutine  TRAVEL  is  called  to  find  a traversal  path  from  the 
entry-point  node  to  the  exit-point  node,  using  the  segments  in  these  arrays.  If 
an  error  is  encountered,  an  error  message  is  printed.  If  no  error  has  been  en- 
countered, subroutine  OPTPATH  is  called  to  improve  the  traversal  path.  The 
node  numbers  of  the  starting  and  ending  nodes  are  stored  in  variables  ISTART 
and  I STOP. 

The  data  describing  paths  to  and  from  the  landfill  and  garage  are 
read  from  TAPE7  by  a BUFFER  IN  statement.  If  IERR  indicates  the  presence  of 
an  error,  control  transfers  to  statement  460.  Otherwise,  if  the  current  trip 
starts  at  the  landfill,  control  transfers  to  statement  410.  If  not,  subrou- 
tine TRACE  is  called  to  obtain  a path  from  the  starting  node  to  the  garage, 
where  travel  will  be  away  from  the  garage.  Control  transfers  to  statement 
420.  At  statement  410,  subroutine  TRACE  is  called  to  obtain  the  path  from  the 
starting  node  to  the  landfill,  where  travel  will  be  away  from  the  landfill.  At 
statement  420,  subroutine  FLIP  is  called  to  reverse  the  order  of  the  segments 
in  the  path.  It  is  called  again  to  reverse  the  order  of  the  nodes.  Subrou- 
tine TRACE  is  called  again  to  obtain  the  path  from  the  exit  node  to  the  land- 
fill. 

Function  CUMDIS  evaluates  the  distance  for  the  three  parts  of  the 
trip.  The  total  distance  is  stored  in  variable  DIST.  If  DIST  is  greater  than 
or  equal  to  SAVDIS,  control  transfers  to  statement  460.  Otherwise,  SAVDIS  is 
set  equal  to  DIST.  Each  of  the  loops  through  statements  430,  440,  and  450 
saves  one  of  the  three  parts  of  the  trip.  The  number  of  segments  in  each 
part  and  the  length  of  each  part  are  also  saved. 
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Statement  460  marks  the  end  of  the  loop  that  selects  the  exit  node. 
Statement  470  marks  the  end  of  the  loop  that  selects  the  entry  node. 

At  statement  480,  If  a valid  trip  was  found,  control  transfers  to 
statement  490.  Otherwise,  an  error  message  is  printed  and  control  transfers 
to  statement  500.  Starting  at  statement  490,  the  three  parts  of  the  saved 
trip  are  written  to  file  TAPE9.  This  trip  has  the  shortest  total  distance 
for  all  of  the  entry-  and  exit-node  combinations  that  have  been  examined. 

Statement  500  marks  the  end  of  the  loop  that  selects  either  the 
garage  or  the  landfill  as  the  start  of  the  trip.  Statement  800  marks  the 
end  of  the  loop  that  generates  routes  for  the  sections.  An  end-of-file  is 
written  on  TAPE 9.  Subroutine  EXIT  is  called  to  terminate  program  execution. 
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SECTION  IV 
INPUT  AND  OUTPUT 


1.  INPUT 

Input  to  program  PHASE3  consists  of  card  input  and  four  disk  files.  Two 
of  the  disk  files  are  generated  by  program  RCINPT,  and  two  are  generated  by 
program  PHASE2.  The  five  files  are  read  at  or  near  the  beginning  of  main  pro- 
gram PHASE3. 

a.  Card  Input 

The  form  and  contents  of  the  three  types  of  data  cards  are  shown  in 
Table  3.  Types  1 and  2 are  required.  The  first  contains  a title  that  is 
printed  on  the  first  line  of  the  output.  The  second  contains  the  node  numbers 
of  the  landfill  and  garage.  The  third  card  is  optional  and  may  be  repeated  up 
to  500  times.  It  contains  a segment  number  and  a section  number.  The  section 
number  indicates  the  section  to  which  the  indicated  segment  is  to  be  reassigned. 

b.  Disk  Files 

Disk  file  TAPE1  contains  segment  data  and  the  map  distance  conversion 
factor  for  the  overall  map.  All  of  the  data  are  read  by  one  binary  READ  state- 
ment. The  following  is  the  list  used  in  the  READ  statement: 

NSEG, ( (STG(J, I ) , J=1 ,8) .DUMMY .DUMMY ,STG(NSF, I ) ,1=1 ,NSEG) ,AVMD 

The  first  word  is  the  count  of  the  segments.  The  segment  data  follow,  11  words 
per  segment.  After  the  segment  data  comes  the  overall  distance  conversion  fac- 
tor. The  STG  array  is  equivalenced  to  the  ISTG  array,  so  some  of  the  contents 
in  the  STG  list  are  integers. 


TABLE  3.  PHASE3  DATA  CARDS 


Columns 

Format 

1-80 

8A10 

1-5 

15 

6-10 

15 

Contents 


Node  number  of  landfill 
Node  number  of  garage 


The  following  card  is  optional  and  may  be  repeated  up  to 
500  times.  I I 


Segment  number 
New  section  number 


TAPE1  contains  the  following  data: 

NSE6  = number  of  segments 

ISTG(1,J)  = street  number 
ISTG(2,J)  = starting  node  number 
ISTG(3,J)  = ending  node  number 
STG(4,J)  = street  length,  in  miles 

ISTG(5,J)  = number  of  houses  These  are  repeated 

STG(6,J)  = speed  limit,  in  mph  for  each  segment. 

ISTG (7 , J ) = number  of  ways  of  travel 

STG(8,J)  = refuse  quantity  adjustment  factor 
STG(9,J)  = x-coordinate  of  segment  midpoint 
STG( 1 0, J ) = y-coordinate  of  segment  midpoint 

ISTG(ll.J)  = shape  code  __ 

AVMD  = map  distance  conversion  factor 

The  index,  J,  corresponds  to  the  Jth  segment  and  is  also  the  segment  number. 
The  segment  midpoint  coordinates  are  not  needed  by  PHASE3  and  are  read  into 
variable  DUMMY.  The  map  distance  conversion  factor  is  not  used  by  the  pro- 
gram either,  since  no  maps  are  drawn. 

Disk  file  TAPE2  contains  refuse  quantity  information  and  node  data. 
All  of  the  data  are  read  by  one  binary  READ  statement.  The  following  is  the 
list  used  in  the  READ  statement: 

NHTOT ,T0TREF , KNODES , (NODNUM ( I ) ,NBRSEG (I ) , XNOD ( I ) , YN0D( I ),I=1 , KNODES ) 
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The  first  three  words  are  the  total  number  of  houses,  the  total  refuse  quantity, 
and  a count  of  the  nodes.  The  node  data  follow,  four  words  per  node. 

TAPE2  contains  the  following  data: 

NHTOT 
TOTREF 
KNODES 
NODNUM(I) 

NBRSEG(I) 

XNOD(I) 

YNOD(I) 

The  index,  I,  corresponds  to  the  Ith  node. 

File  TAPE3  contains  pointers  to  the  first  and  last  segments  in  each 
section  and  the  vehicle  capacity  for  each  section.  All  of  the  data  are  read 
by  one  binary  READ  statement.  The  following  list  is  used  in  the  READ  state- 
ment: 

NSEG,NTRIP, (SINL(I ) ,SFNL( I ) ,TC( I ) ,1=1 ,NTRIP) 

The  first  two  words  are  a count  of  the  segments  and  a count  of  the  sections. 

The  initial  segment  pointer,  the  final  segment  pointer,  and  the  vehicle  capa- 
city follow  and  are  repeated  for  each  section. 

TAPE3  contains  the  following  data: 

NSEG  = count  of  segments 

NTRIP  = count  of  sections 

SINL(I)  = sequence  of  number  of  first  segment  in  section  I 
SFNL(I)  = sequence  of  number  of  last  segment  in  section  I 

TC ( I ) = capacity  of  vehicle  servicing  section  I 

The  pointers  refer  to  the  position  of  the  segment  numbers  on  file  TAPE4. 

File  TAPE4  holds  the  segment  numbers  in  the  order  in  which  they  are 
assigned  to  sections.  The  segment  numbers  are  read  one  at  a time  by  a for- 
matted READ  statement,  using  format  16.  The  numbers  are  stored  in  the  ISTG 
array.  The  segment  numbers  are  grouped  by  section.  All  segments  in  section 
1 are  listed,  followed  by  all  segments  in  section  2,  and  so  forth.  The  seg- 
ment numbers  within  each  section  are  not  necessarily  listed  in  order.  Instead, 
they  occur  in  the  order  in  which  they  were  selected  for  addition  to  the  section 
by  program  PHASE2. 


total  number  of  houses  or  stops 
total  refuse  quantity 
number  of  nodes 
node  number 

packed  bounding  segment  numbers 
x-coordinate  of  node 
y-coordinate  of  node 


These  are  repeated 
for  each  node. 
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2.  SCRATCH  FILE 


W»  l.w 


Disk  file  TAPE7  is  used  as  a scratch  file.  The  data  are  first  written  to 
TAPE7  by  a BUFFER  OUT  statement  with  parity  1 at  statement  86  of  PHASE3.  The 
region  written  extends  from  the  beginning  of  array  NNTF  through  variable  XXX. 
The  saved  data  include  arrays  NNTF,  DTF,  TTF,  and  NSTF.  Each  array  is  a 
double-subscripted  array  with  second  dimension  3.  The  data  in  these  arrays 
describe  paths  from  the  landfill  node  to  each  other  node,  traveled  away  from 
the  landfill;  paths  from  the  garage  node  to  each  other  node;  and  paths  from 
the  landfill  node  to  each  other  node,  traveled  towards  the  landfill.  Array 
NNTF  contains  line  numbers  of  the  next  node  in  each  path.  Array  DTF  contains 
distances  to  or  from  each  node.  Array  TTF  contains  travel  times  to  or  from 
each  node.  Array  NSTF  contains  the  segment  numbers  for  the  next  segments  in 
each  path. 

The  data  are  saved  on  disk  to  free  storage  for  other  purposes  during  the 
execution  of  the  program.  The  data  are  read  by  a BUFFER  IN  statement  in  sub- 
routine EPXP  and  at  statement  390  of  the  main  program.  In  both  cases,  paths 
from  the  landfill  or  garage  to  the  section  are  needed. 


3.  OUTPUT 

Program  PHASF3  produces  disk  and  printed  output.  The  descriptions  of  the 
final  routes  are  written  to  file  TAPE9.  Information  is  printed  to  allow  the 
user  to  verify  user-initiated  changes  in  section  assignments  and  to  permit  de- 
bugging if  errors  should  occur. 

a.  Disk  Output 

Disk  file  TAPE9  contains  information  describing  travel  and  collection 
for  each  vehicle.  The  data  are  written  in  card-image  form  by  eight  formatted 
WRITE  statements  in  main  program  PHASE3.  A description  of  the  path  from  the 
landfill  to  the  garage  is  written  at  the  beginning  of  the  file.  Two  paths 
follow  for  each  section,  one  starting  at  the  landfill  and  one  starting  at  the 
garage. 
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Each  trip  is  broken  into  three  pieces.  The  first  piece  describes 
travel  from  the  garage  or  landfill  to  the  collection  region.  The  second  piece 
describes  travel  and  collection  in  the  collection  region.  The  third  piece  de- 
scribes travel  from  the  last  collection  point  to  the  landfill.  There  is  a 
header  card  for  each  piece.  Each  header  card  contains  a segment  count,  the 
total  distance  in  miles  for  the  piece,  the  section  number,  the  vehicle  capa- 
city, and  the  amount  of  refuse  to  be  collected  from  the  section.  The  header 
card  is  followed  by  one  or  more  cards  that  describe  the  path  for  that  piece 
of  the  trip.  The  path  is  described  by  repeating  three  items:  a node  number, 
a segment  number,  and  the  letter  T or  C to  indicate  travel  or  collection  on 
that  segment.  The  format  used  for  each  card  image  is  8(15, 14 ,A1 ) . Up  to 
eight  segments  can  be  described  on  one  card  image.  The  number  of  the  node 
terminating  the  piece  of  the  trip  follows  the  three  items  describing  the  last 
segment. 

The  header  and  path  data  for  the  path  from  the  landfill  to  the  garage 
are  written  shortly  before  statement  86  of  the  main  program.  The  header  in- 
cludes only  the  number  of  segments  in  the  path  and  the  length  of  the  path.  The 
headers  and  path  descriptions  for  the  three  pieces  of  each  trip  are  written  by 
six  WRITE  statements  from  statement  490  to  statement  500  in  the  main  program. 

The  data  on  TAPE9  are  in  card-image  form  so  they  may  be  punched  in- 
stead of  placed  on  disk.  This  feature  allows  the  user  to  modify  trips  before 
the  data  are  used  by  program  PHASE4.  File  TAPE9  can  also  be  copied  to  the 
output  file  for  visual  examination  prior  to  use. 

b.  Printed  Output 

The  printed  output  consists  of  five  sections:  a listing  of  initial 
section  assignments,  a listing  of  changes  to  section  assignments,  a summary  of 
segment  data,  a summary  of  node  data,  and  a listing  of  suitable  entry  and  exit 
points  for  each  section.  Samples  of  each  type  of  printed  output  are  presented 
in  Appendix  D. 


The  problem  title  is  printed  at  the  beginning  of  the  printed  output. 
The  segment  assignments  produced  by  program  PHASE2  follow.  The  assignments  are 
listed  by  section  in  the  order  in  which  the  segments  occur  on  file  TAPE4.  Up 
to  30  segment  numbers  are  printed  per  line. 


Program  and  user-initiated  section  changes  are  printed  next.  The 
five  columns  of  information  include  the  segment  number,  the  initial  section 
assignment,  the  new  section  assignment,  the  amount  of  refuse  on  that  segment, 
and  the  type  of  adjustment  (program-adjusted  or  user-adjusted).  The  program- 
initiated  changes  are  printed  first.  A summary  of  the  vehicle  loads  and  capa- 
cities for  each  trip  are  printed  following  the  program-initiated  changes.  The 
segment  number  and  new  section  number  in  the  user-initiated  changes  are  read 
from  data  cards.  The  final  vehicle  capacities  and  loads  for  each  trip  are 
printed  following  the  user-initiated  section  reassignments. 


The  first  column  of  segment  data  gives  the  segment  number.  The  sec- 
ond column,  headed  NSTR,  gives  the  street  number  for  the  segment.  Columns  NN1 
and  NN2  give  the  starting  and  ending  node  numbers.  Column  LEN  gives  the  seg- 
ment length  in  miles.  Column  NH  gives  the  number  of  houses  on  the  segment. 
Column  SPD  gives  the  speed  limit  on  the  segment.  Column  NWAY  gives  the  number 
of  directions  of  travel  allowed  on  the  segment.  Column  RQF  gives  the  refuse 
quantity  adjustment  factor.  Column  NSCT  gives  the  final  section  assignment  of 
the  segment.  Segments  containing  no  refuse  have  a zero  as  their  section  as- 
signment. The  last  column,  headed  SF,  gives  the  shape  code  in  octal. 


The  node  numbers  of  the  landfill  and  garage  are  given  at  the  begin- 
ning of  the  node  data  list.  The  first  column,  headed  LINE,  is  the  line  number 
of  the  node.  The  second  column,  headed  NODE,  is  the  node  number.  The  columns 
headed  NNFD,  NNTD,  and  NNFG  give  the  line  numbers  of  the  next  node  on  the 
shortest  path  (1)  to  the  landfill  traveling  away  from  the  landfill,  (2)  to  the 
landfill  traveling  toward  the  landfill,  and  (3)  to  the  garage  traveling  away 
from  the  garage.  Columns  NSFD,  NSTD,  and  NSFG  give  the  segment  numbers  of  the 
next  segment  in  each  of  the  three  types  of  paths.  Columns  DFD,  DTD,  and  DFG 
give  the  distance  of  the  node  from  the  landfill  or  garage  for  the  three  types 
of  paths.  The  travel  times,  in  minutes,  are  given  in  the  next  three  columns, 
headed  TFD,  TTD,  and  TFG.  Data  for  paths  traveled  toward  and  away  from  the 


landfill  are  different  only  when  one-way  streets  cause  the  travel  paths  to  be 
different.  The  column  headed  NUMNBR  gives  a count  of  segments  bounding  the 
node.  The  remaining  six  columns,  headed  NEIGHBORING  SEGMENT  NUMBERS,  give  the 
numbers  of  the  segments  bounding  the  node.  Some  of  the  columns  will  be  blank 
when  a node  has  fewer  than  six  neighboring  segments. 

Entry  and  exit  -nodes  are  listed  for  each  section.  The  line  starting 
FROM  DUMP  gives  the  numbers  of  nodes  that  are  suitable  entry  points  to  the  sec- 
tion for  a vehicle  starting  at  the  landfill.  The  line  starting  FROM  GARAGE 
gives  the  numbers  of  nodes  that  are  suitable  entry  points  to  the  section  for  a 
vehicle  coming  from  the  garage.  The  line  starting  TO  DUMP  gives  the  numbers 
of  nodes  that  are  suitable  points  for  teaving  the  section  to  travel  to  the 
landfill.  From  one  to  ten  node  numbers  may  be  listed  on  each  line.  Variable 
JMAX  in  the  main  program  limits  the  number  of  entry  points,  and  variable  KMAX 
limits  the  number  of  exit-point  selections.  At  present,  both  JMAX  and  KMAX 
are  set  to  4. 

If  file  TAPE9  is  copied  to  the  output  file,  the  printout  will  follow 
the  entry-  and  exit-point  listing. 
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SECTION  V 

PROGRAM  REQUIREMENTS 


1.  SYSTEM 

Program  PHASE3  is  written  entirely  in  FORTRAN  IV.  The  program  runs  on  a 
CDC  6600  computer  using  a SCOPE  3.4.4  operating  system. 

Nine  obvious  types  of  computer-dependent  coding  occur  in  program  PHASE3 
and  its  subroutines.  Main  program  PHASE3  assumes  a 130-character  output  line, 
10  characters  per  word,  and  a 60-bit  word.  System  function  SHIFT  is  used  in 
the  main  program  and  in  subroutines  ADJUST,  TREE,  C0N2ST,  CONNECT,  DISCON,  and 
CL0SE1.  Masking  operations  are  used  in  subroutines  TREE,  C0N2ST,  CONNECT, 
DISCON,  and  CL0SE1.  Asterisk-bounded  text  is  used  in  format  statements  in  the 
main  program  and  in  subroutines  ADJUST,  PRNPCS,  TRACE,  TREE,  FNDPTH,  CONNST, 
DISCON,  EPXP,  SOLV,  and  TRAVEL.  An  Iw.n  format  specification  is  used  in  sub- 
routines PRNPCS  and  EPXP  and  in  the  main  program.  Multiple  replacement  state- 
ments occur  in  the  main  program  and  in  subroutines  M0VE3,  C0N2ST,  FNDPTH, 
CONNECT,  EPXP,  CL0SE1,  GENDM,  TRAVEL,  and  OPTPATH.  A dollar  sign  is  used  to 
separate  FORTRAN  statements  in  the  main  program  and  all  subprograms  except 
subroutine  GENDM,  function  IFIND,  and  function  CUMDIS. 

More  subtle  types  of  machine  dependencies  may  exist,  according  to  the 
machine  used. 

2.  STORAGE 

The  core  requirement  is  slightly  less  than  123,0008  words.  Exactly  3,601 
words  of  data  are  written  to  file  TAPE7.  The  maximum  peripheral  storage  used 
by  file  TAPE9  should  not  exceed  33,800  words. 
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3.  TIME 


The  running  time  for  PHASE3  depends  on  the  number  of  sections  and  the  num- 
ber of  entry  and  exit  points.  Since  the  number  of  entry  and  exit  points  cannot 
be  determined  by  the  user  before  the  program  is  run,  a precise  estimate  of  run- 
ning time  cannot  be  given.  An  upper  limit  on  the  CP  time,  in  minutes,  should 
be  0.36  times  the  number  of  sections.  Very  rough  limits  on  the  10  and  PP  times 
are  0.06  minutes  10  time  per  section  and  0.12  minutes  PP  time  per  section. 
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SECTION  VI 
PROGRAM  LIMITATIONS 


The  only  new  limitation  imposed  upon  the  program  user  by  program  PHASE3 
is  a limit  of  500  on  the  number  of  section  reassignment  cards.  More  cards  are 
permitted,  but  since  there  can  be  no  more  than  500  segments,  there  would  be  no 
need  for  additional  reassignment  cards. 

Limitations  imposed  on  the  map  description  by  the  dimensions  in  PHASE3  are 
500  segments  and  300  nodes.  The  number  of  trips  is  limited  to  50.  These  limi- 
tations should  have  been  observed  when  programs  RCINPT  and  PHASE2  were  run.  The 
segments  connected  to  the  landfill  and  garage  must  not  include  houses. 

Four  limitations  exist  which  may  cause  problems  that  the  user  usually 
cannot  anticipate:  (1)  One-way  streets  may  prevent  the  program  from  finding 
a travel  path  for  a trip,  (2)  storage  has  been  dimensioned  to  allow  pairing 
of  at  most  32  odd  nodes,  (3)  array  dimensions  limit  each  of  the  three  pieces 
of  each  trip  to  100  segments,  and  (4)  array  dimensions  limit  the  total  number 
of  steps  in  the  three  pieces  to  200.  Corrective  procedures  for  dealing  with 
these  limitations  are  discussed  following  the  corresponding  error  messages 
in  the  next  section. 
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SECTION  VII 

ERROR  MESSAGES  AND  CORRECTIVE  ACTION 


1 . STOP  30 

Type:  Fatal. 

Source:  Subroutine  CONNECT. 

Location:  Day  file. 

Meaning:  No  segments  have  been  found  in  section  NTRIP  by  subroutine  CONNECT. 

The  data  on  file  TAPE3  or  on  file  TAPE4  may  be  incorrect,  or  the 
user  may  have  removed  all  segments  from  the  section  by  using  section 
reassignment  cards. 

Action:  If  the  user  has  intentionally  removed  the  section,  the  STOP  30  state- 

ment in  subroutine  CONNECT  should  be  replaced  by  a RETURN  statement. 
If  the  problem  has  been  created  unintentionally,  the  user  should  de- 
termine whether  files  TAPE3  and  TAPE4  correspond  to  the  current  prob- 
lem. If  all  segment-related  data  are  correct,  the  error  may  be  com- 
pletely unrelated  to  the  section  assignments.  Other  indicated  errors 
should  be  corrected. 

2.  STOP  111 

Type:  Fatal. 

Source:  Subroutine  CONNECT. 

Location:  Day  file. 

Meaning:  The  number  of  a starting  or  ending  node  in  the  segment  data  has  not 

been  found  in  array  NODNUM. 

Action:  Verify  that  the  segment  data  and  node  data  permanent  files  corre- 

spond to  the  same  problem.  Correct  any  other  errors  indicated  by 
the  program.  If  the  problem  persists,  a core  dump  of  the  program 
will  be  necessary.  Variable  NNN(J)  in  subroutine  CONNECT  should  be 
found.  If  the  number  is  a valid  node  number,  then  it  should  also 
appear  in  array  NODNUM.  If  it  does  not  appear  in  array  NODNUM,  a 
programmer  should  determine  why  it  is  not  there.  If  NNN(J)  is  not 
a valid  node,  the  programmer  should  determine  how  it  was  obtained 
from  the  ISTG  array. 
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3.  STOP 

Type: 

Source: 

Location: 

Meaning: 

Action: 


4.  STOP 

Type: 

Source: 

Location: 

Meaning: 

Action: 


5.  STOP 

Type: 

Source: 

Location: 

Meaning: 

Action: 


201 

Fatal . 

Subroutine  CL0SE1. 

Day  file. 

A node  number  in  the  segment  data  array  could  not  be  found  in  array 
NODNUM. 

Verify  that  the  node  and  segment  data  permanent  files  belong  to  the 
same  problem.  Correct  any  other  indicated  error.  If  the  problem 
persists.  It  may  be  caused  by  an  unrelated  program  logic  error.  A 
full  core  dump  may  be  necessary  to  find  the  cause. 

230 
Fatal . 

Subroutine  CONNECT. 

Day  file. 

More  than  NPIECE  pieces  were  found  in  section  NTRIP  in  the  loop 
through  statement  230  in  subroutine  CONNECT.  The  piece  count, 
NPIECE,  is  probably  in  error. 

The  computation  of  NPIECE  should  be  examined  for  errors.  NPIECE  is 
set  to  1 at  the  beginning  of  subroutine  CONNECT  and  is  incremented 
following  statement  70.  If  no  error  is  found  in  the  computation, 
the  temporary  section  assignments  in  the  ISTG  array  should  be  ex- 
amined to  find  the  cause  of  an  invalid  piece  number  in  these  data. 

270 
Fatal . 

Subroutine  CONNECT. 

Day  file. 

Piece  number  NTRIP  could  not  be  found  in  the  NNFN  array  in  the  loop 
through  statement  280  in  subroutine  CONNECT.  An  error  in  program 
logic  or  a machine  malfunction  probably  has  occurred. 

Any  other  indicated  errors  should  be  corrected.  Variable  NTRIP  and 
the  first  NPIECE  entries  in  the  NNFN  array  should  be  examined,  either 
from  a dump  or  from  a user-inserted  print  statement.  If  NTRIP  is  not 
present  in  the  NNFN  array,  the  error  has  occurred  between  statements 
240  and  280.  A programmer  should  find  and  correct  the  omission. 
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6.  STOP  501 


Type: 

Source: 

Location: 

Meaning: 


Action: 


7.  STOP 

Type: 

Source: 

Location: 

Meaning: 

Action: 


8.  STOP 

Type: 

Source: 

Location: 

Meaning: 


Fatal . 

Subroutine  CL0SE1. 

Day  file. 

A node  in  section  NTRIP  has  no  neighboring  segments  in  that  section. 
Either  the  node  should  not  be  in  array  NODPTR,  or  one  of  its  bound- 
ing segments  should  be  in  the  section. 

Verify  that  the  node  data  and  segment  data  permanent  files  correspond 
to  the  same  problem.  If  the  error  is  not  in  the  original  node  or 
segment  data,  a program  logic  error  is  probably  the  cause.  The  se- 
lection of  nodes  in  the  loop  through  statement  220  in  the  main  pro- 
gram should  be  reexamined. 

1401 
Fatal . 

Subroutine  CL0SE1. 

Day  file. 

More  than  32  steps  (segments)  are  needed  to  return  to  the  section 
from  an  order-one  node  at  the  end  of  a one-way  street. 

If  the  map  description  is  correct,  the  segment  should  be  reassigned 
to  a different  section.  More  likely,  the  map  description  is  incor- 
rect or  consists  of  two  separate  regions  connected  by  a one-way 
street.  The  description  should  be  corrected,  or  additional  segments 
should  be  added  to  the  map  description  to  connect  the  two  pieces. 

1601 
Fatal . 

Subroutine  CL0SE1 . 

Day  file. 

The  line  number,  LF,  of  the  node  at  which  a path  returns  to  the  sec- 
tion cannot  be  found  in  the  NODPTR  array  in  the  loop  through  state- 
ment 160.  The  program  logic  is  in  error,  probably  in  the  selection 
of  the  return  path. 
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Action: 


The  path  data  in  arrays  IPATHS  and  IPATHN  should  be  examined.  The 
arrays  are  filled  in  the  loop  through  statement  110.  The  type  of 
error  in  the  path  should  indicate  the  probable  source  of  the  error. 


NODE  nnnnn  IS  NOT  IN  THE  NODE  TABLE 
CORRECT  THE  DUMP  OR  GARAGE  NUMBER 
JOB  TERMINATED 


Note:  nnnnn  is  a number  from  1 to  99999. 

Type:  Fatal. 

Source:  Subroutine  TREE. 

Location:  This  message  may  be  printed  in  one  of  three  places:  following  the 

section  ^assignments,  following  the  node  data,  or  following  the 
entry-  and  exit-point  listing. 

Meaning:  If  the  message  is  printed  immediately  following  the  section  reassign- 

ments, then  the  node  number  of  the  landfill  or  garage  specified  on 
the  second  data  card  does  not  match  any  node  number  in  the  node  data. 
Otherwise,  a node  number  from  a starting  or  ending  node  in  the  seg- 
ment data  has  not  been  found  in  the  node  data. 

Action:  If  the  error  message  is  printed  immediately  following  the  section  re- 

assignments, the  node  number  of  the  landfill  or  garage  on  the  second 
data  card  should  be  corrected.  Otherwise,  verify  that  the  segment 
and  node  data  permanent  files  correspond  to  the  same  problem. 

10.  IMPROPER  LINE  POINTER  USED  IN  TRACE 
JOB  TERMINATED 


Action: 


Type:  Fatal. 

Source:  Subroutine  TRACE. 

Location:  This  message  may  occur  in  one  of  two  places:  immediately  following 

the  section  reassignments,  or  following  the  entry-  and  exit-point 
listing. 

Meaning:  Either  the  line  number,  L,  is  outside  the  range  of  permissible  node 

line  numbers  (1  through  KNODES),  or  the  number  of  segments  in  the 
path,  NSG,  exceeds  200. 

Action:  If  the  message  is  immediately  after  the  section  reassignments,  the 

problem  is  in  the  trip  from  the  garage  to  the  landfill.  If  the 
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message  follows  the  entry-  and  exit-point  printout,  the  problem  may 
be  in  the  path  between  any  two  nodes.  The  numbers  of  the  nodes  that 
start  and  end  the  path  are  indicated  in  the  error  message  as  vari- 
ables NSTART  and  NSTOP.  If  the  number  of  segments,  NS6,  exceeds  200, 
the  user  should  compare  the  path  to  the  input  maps.  If  more  than  200 
steps  are  required,  the  dimensions  of  arrays  IPS  and  IPN  in  the  main 
program  should  be  increased.  If  the  path  is  incorrect,  or  if  the 
line  pointer,  L,  is  improper,  verify  that  the  node  and  segment  data 
permanent  files  correspond  to  the  same  problem.  If  the  error  per- 
sists after  all  other  errors  have  been  corrected,  the  user  must  find 
the  cause.  A core  dump  is  provided  after  the  error  message. 


11.  NODE  nnnnn  DOES  NOT  CONNECT  TO  THE  a 

CORRECT  THE  DUMP  OR  GARAGE  NUMBERS  IF  INCORRECT  OR  ELSE 
CORRECT  THE  MAP  DESCRIPTION  CARDS  AND  RERUN  THE  INPUT  AND 
SECTIONING  PROGRAMS 

JOB  TERMINATED 


Note : 

Type: 

Source: 

Location 

Meaning: 

Action: 


nnnnn  is  a node  number  from  1 through  99999. 
a is  either  DUMP,  GARAGE,  or  DUMP  OR  TO  THE  GARAGE. 
Fatal . 

Main  program  PHASE3. 

Following  the  node  listing. 

Self  explanatory. 

Self  explanatory. 


12.  MORE  THAN  mm  STEPS  NEEDED  TO  CONNECT  PIECES  ppppp  AND  ppppp 

IF  THIS  IS  CORRECT,  INCREASE  THE  FIRST  DIMENSION  OF  IPATH  IN 

COMMON  BLOCK  TEMSTG,  SUBROUTINES  FNDPTH  AND  CONNST 

JOB  TERMINATED 

Note:  In  the  current  version  of  the  program,  mm  will  be  50.  ppppp  repre- 

sents a piece  number.  The  tens  and  units  digits  give  the  section 
number,  and  the  thousands  and  ten- thousands  digits  give  the  piece 
number.  Additional  printout  will  give  the  node  and  piece  numbers 
of  the  end  points  and  the  first  15  segments  in  the  path. 

Type:  Fatal. 

Source:  Subroutine  FNDPTH. 
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Location:  Following  the  node  data  listing. 

Meaning:  Self  explanatory. 

Action:  If  a path  shorter  than  50  steps  exists  between  the  two  pieces,  an 

error  In  the  program  logic  is  present,  probably  In  subroutine 
FNDPTH.  A programmer  should  seek  the  error. 

13.  NO  CONVERGENCE  IN  CONNST 

Type:  Warning. 

Source:  Subroutine  CONNST. 

Location:  Following  the  node  data  listing. 

Meaning:  No  paths  could  be  found  between  some  or  all  of  the  pieces  of  the 

current  section.  The  connections  may  be  one-way  streets,  or  the 
map  description  may  be  in  error.  An  error  may  exist  in  the  logic 
in  subroutine  CONNST. 

Action:  Obtain  a core  dump.  The  value  of  variable  ITRIP  in  subroutine 

CONNECT  gives  the  number  of  the  section  containing  the  problem. 

Examine  the  pieces  of  this  section,  using  a map  from  program 
PHASE2  and  the  section  assignments  in  the  segment  list.  If  errors 
are  found  in  the  map  description,  they  should  be  corrected.  If 
one-way  streets  are  the  only  paths  between  pieces  of  the  section, 
the  streets  must  be  made  two-way  streets.  The  paths  generated  by 
PHASE3  on  file  TAPE9  should  be  punched  and  adjusted  by  hand  to  con- 
form to  the  one-way  streets.  If  no  problems  are  found  in  the  map 
description  and  if  no  one-way  problems  exist,  a programmer  should 
examine  subroutine  CONNST  for  logic  errors. 

14.  CONNST  ITER  nn 

Note:  In  the  current  version  of  the  program,  nn  will  be  11. 

Type:  Warning,  but  the  results  from  PHASE3  cannot  be  used  by  PHASE4 

until  the  problem  has  been  corrected. 

Source:  Subroutine  CONNST. 

Location:  Following  the  node  data  listing. 

Meaning:  Subroutine  CONNST  has  been  unable  to  connect  all  of  the  pieces  of  a 

section  in  10  iterations.  Since  at  least  five  pieces  are  connected 
by  each  iteration,  there  are  at  least  50  pieces  to  the  section 
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requiring  more  than  two  step  connections,  or  the  logic  in  subrou- 
tine CONNST  is  in  error. 

Action:  Verify  the  map  description  as  for  error  type  13.  If  the  section 

consists  of  more  than  50  pieces,  section  reassignment  cards  should 
be  used  to  make  the  section  more  compact.  If  the  subroutine  logic 
is  in  error,  a programmer  should  seek  the  error. 


15.  I = ddd  NODE  = nnnn 


Note:  ddd  represents  a subscript  from  1 through  100,  and  nnnn  represents 

a node  number.  A table  of  node  pointers  and  node  orders  follows 
the  message,  and  a core  dump  follows  the  table. 

Type:  Fatal. 

Source:  Subroutine  DISCON. 


Location:  Following  the  node  data  listing. 

Meaning:  An  end-point  node  of  a segment  could  not  be  found  in  the  node  data. 

Variable  I indicates  the  subscript  of  the  node  in  the  node  pointer 
and  node  order  arrays.  Variable  NODE  gives  the  node  number. 

Action:  Verify  that  the  node  and  segment  data  permanent  files  correspond  to 

the  same  problem.  Correct  any  other  indicated  errors.  If  the 
error  persists,  a programmer  should  examine  subroutine  DISCON  for 
logic  errors. 


16.  LINE  NODE  NXND  NXSG  DIST  TIME  SECT 


Note:  A table  of  data  follows  the  indicated  column  headings.  A core  dump 

follows  the  table. 

Type:  Fatal. 

Source:  Subroutine  CL0SE1 . 

Location:  Following  the  entry-  and  exit-point  listing. 

Meaning:  A break  has  occurred  in  a path  generated  by  subroutine  TREE  from  an 

order-one  node  back  to  the  current  section.  A malfunction  in  sub- 
routine TREE  is  the  most  probable  cause. 

Action:  All  other  errors  should  be  corrected.  If  the  problem  persists,  the 

path  should  be  examined  for  possible  map  problems.  The  pointer  to 
the  starting  node  of  the  path  can  be  found  in  the  dump  as  NODPTR(K), 
where  K is  a variable  in  subroutine  CLOSEl . The  number  of  the  final 
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node  of  the  path  is  variable  NORG  in  subroutine  CL0SE1.  To  trace 
the  path,  find  the  line  number  of  the  starting  node  in  the  column 
headed  LINE.  The  number  of  the  node  is  in  the  column  headed  NODE. 
The  next  segment  in  the  path  is  in  the  column  headed  NXSG.  The 
section  assignment  of  the  segment  is  in  the  column  headed  SECT. 

The  line  number  of  the  next  node  in  the  path  is  in  the  column 
headed  NXND.  Columns  DIST  and  TIME  give  the  distance  and  time 
from  the  node  in  column  NODE  to  node  number  NORG.  The  break  in 
the  path  occurs  where  NXND  refers  to  a line  that  has  not  been 
printed,  or  where  NXSG  contains  a zero,  in  the  table  that  com- 
prises part  of  this  error  message. 


17.  THE  NUMBER  OF  ODD  NODES,  nnn,  EXCEEDS  mum 

EITHER  MODIFY  THE  MAP  TO  PAIR  OR  REMOVE  ODD  NODES  OR  ELSE 
SEE  THE  PROGRAM  DESCRIPTION  MANUAL  FOR  INSTRUCTIONS 
NODE/ORDER 


Note: 


Type: 

Source: 

Location: 

Meaning: 

Action: 


Up  to  10  node/order  pairs  per  line  are  printed  following  the  NODE/ 
ORDER  heading.  A list  of  segments  in  each  piece  of  the  section 
follows  this  tabulation. 

Warning,  unless  error  message  20  follows. 

Main  program  PHASE3. 

Following  the  entry-  and  exit-point  listing. 

The  section  contains  more  odd  nodes  than  the  program  is  dimen- 
sioned to  handle. 

No  action  is  necessary  unless  error  message  20  follows  this  mes- 
sage. If  error  message  20  follows,  there  are  two  alternatives: 
either  modify  the  section,  or  increase  the  dimensions  of  arrays 
used  in  odd-node  pairing.  The  segment  list  will  indicate  which 
segments  are  in  the  section.  The  node  and  order  listing  will 
indicate  which  nodes  are  odd  nodes.  At  two  of  the  nodes,  an  even 
number  of  segments  will  meet;  these  are  the  entry-  and  exit-point 
nodes.  The  section  description  should  be  modified  if  section 
reassignment  cards  can  be  used  to  add  or  remove  segments  in  such 
a manner  that  the  number  of  odd  nodes  is  reduced  to  the  maximum 
(mmm) . 
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If  the  dimensions  used  in  odd-node  pairing  are  to  be  increased, 
changes  will  be  necessary  in  the  main  program  and  in  subroutines 
SELORD,  SOLV , NEXTM,  and  PATH.  The  dimension  of  array  IPART  should 
be  increased  to  the  number  of  odd  nodes  in  COMMON  block  PART  in  the 
main  program  and  in  subroutines  SOLV,  NEXTM,  and  PATH.  The  dimen- 
sions of  arrays  ISTP  and  JSTP  should  be  increased  to  the  number  of 
odd  nodes  in  COMMON  block  STP  in  subroutines  SOLV,  NEXTM,  and  PATH. 
The  second  and  third  dimensions  of  arrays  1ST PO  and  ISTPN  should  be 
increased  to  the  number  of  odd  nodes  in  subroutines  NEXTM  and  PATH. 
In  subroutine  NEXTM,  both  dimensions  of  arrays  MTXN  and  MTXO  should 
be  increased  to  the  number  of  odd  nodes.  In  subroutine  SOLV,  the 
last  two  dimensions  in  each  array  in  COMMON  block  TEMSTG  should  be 
increased  to  the  number  of  odd  nodes.  If  the  dimensions  are  in- 
creased to  allow  more  than  36  odd  nodes,  the  dimension  of  array  Z 
in  COMMON  block  TEMSTG  in  the  main  program  should  be  increased  to 
4000  less  than  six  times  the  square  of  the  number  of  odd  nodes. 

In  the  main  program,  the  dimension  of  array  ONDMTX  should  be  in- 
creased to  the  square  of  the  number  of  odd  nodes.  In  a DATA 
statement,  the  value  of  MAXODD  should  be  increased  to  the  number 
of  odd  nodes. 

18.  NCH=  ddd  QUITTING 

Note:  ddd  will  equal  20  in  the  current  version  of  the  program. 

Type:  Warning. 

Source:  Subroutine  SOLV. 

Location:  Following  the  entry-  and  exit-point  listing. 

Meaning:  The  odd-node  pairing  improvement  did  not  converge  in  20  iterations. 

A one-way  street  is  probably  present  on  a path  that  pairs  two  odd 
nodes. 

Action:  No  action  is  necessary  unless  error  message  20  follows  this  mes- 

sage. If  error  message  20  appears,  either  the  one-way  street  must 
be  changed  to  a two-way  street  in  the  data  to  program  RCINPT,  or 
the  odd  nodes  must  be  made  even  nodes  using  section  reassignment 
cards.  If  the  one-way  street  is  made  a two-way  street,  the  TAPE9 
output  from  PHASE3  should  be  modified  if  final  travel  on  that 
street  is  in  the  wrong  direction. 
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19.  KS  ARRAY 


Note: 


Type: 

Source: 

Location: 

Meaning: 


Action: 


A table  of  100  numbers,  20  per  line,  will  follow  this  heading. 

The  arguments  to  subroutine  TRAVEL  and  data  describing  the  seg- 
ments in  the  current  section  will  be  printed  on  the  following 
page.  The  data  are  terminated  by  the  message  TRAVEL  PATH  CANNOT 
BE  FOUND.  RECHECK  NETWORK  FOR  THIS  TRIP. 

Warning. 

Subroutine  TRAVEL. 

Following  the  entry-  and  exit-point  listing. 

The  trial-and-error  search  for  a travel  path  in  the  collection  re- 
gion has  been  unsuccessful  in  subroutine  TRAVEL.  There  are  three 
possible  causes:  (1)  an  error  in  the  map  description,  (2)  one-way 
streets  causing  problems  in  the  odd-node  pairing,  or  (3)  a logic 
error  somewhere  in  the  program. 

No  action  is  necessary  unless  error  message  20  follows  this  mes- 
sage. If  error  message  20  follows,  the  data  printed  with  this 
message  can  be  used  to  determine  the  type  of  problem. 

The  KS  array  indicates  the  number  of  times  the  subroutine  has 
backed  up  to  each  step  in  the  path  during  the  trial-and-error 
search.  Subroutine  TRAVEL  stops  the  search  when  it  has  backed  up 
to  a point  101  times.  All  of  the  path  up  to  this  point  indicates 
the  current  status  of  the  travel  path.  The  nodes  and  segments  in 
the  path  are  on  the  following  page  of  printed  output  under  the 
column  headings  IPN  and  IPS.  All  of  the  segments  in  the  section 
are  listed  following  the  heading  TRAVEL  ARGUMENTS.  Column  head- 
ings indicate  the  starting  and  ending  node  numbers.  Column  NSG 
gives  the  segment  numbers.  Column  NWAY  indicates  the  number  of 
ways  of  travel  on  the  street.  The  column  headed  TIMES  indicates 
the  number  of  times  the  street  is  available  for  traversal. 

The  path  is  to  start  at  the  node  indicated  by  ISTART  and  end 
at  the  node  indicated  by  IST0P.  The  program  user  should  attempt 
to  find  a travel  path  manually,  using  these  nodes  and  segments,  in 
order  to  determine  the  type  of  problem  encountered  by  subroutine 
TRAVEL.  If  one-way  streets  are  causing  the  problem,  they  should 
be  made  two-way  streets  in  the  map  description  to  program  RCINPT. 
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The  output  from  PHASE3  (file  TAPE9)  should  be  adjusted  if  travel 
on  these  streets  is  in  the  wrong  direction.  If  it  is  impossible 
to  find  a travel  path  manually,  an  error  probably  exists  in  the 
selection  of  the  segments  in  the  section. 

20.  NO  PATH  EXISTS  FOR  THE  TRIP  STARTING  AT  THE  a 

Note:  a will  be  either  GARAGE  or  LANDFILL. 

Type:  Fatal. 

Source:  Main  program  PHASE3. 

Location:  Following  error  messages  after  the  entry-  and  exit-point  listing. 

Meaning:  No  travel  path  could  be  found  in  the  current  collection  region. 

Warning  message  17,  18,  or  19  will  precede  this  error  message. 

Action:  The  corrective  action  is  described  under  the  appropriate  warning 

message. 


21.  I=i  LINE=j  LIM=m  NBRS=n 


Note: 

Type: 

Source: 

Location: 

Meaning: 


Action: 


A table  of  line  numbers,  node  pointers,  and  node  orders  follows  the 
message. 

Fatal . 

Subroutine  DISCON. 

Following  the  node  data  listing. 

NODPTR(i)  points  to  a node,  NODNUM(j),  which  should  be  in  the  sec- 
tion, but  none  of  the  m segments  bounding  the  node  are  in  the  sec- 
tion. The  segment  numbers  are  packed  into  each  ten  bits  of  the 
octal  number  n.  The  node  and  segment  numbers  might  be  from  dif- 
ferent problems,  or  the  node  pointers  in  array  N00PTR  may  have 
been  overwritten. 

Verify  that  the  node  and  segment  data  correspond  to  the  same  prob- 
lem. If  the  number  of  nodes  in  the  section  (the  line  number  of 
the  last  line  in  the  table  gives  the  number  of  nodes  in  the  sec- 
tion) exceeds  100,  the  node  pointer  array  has  been  partially  over- 
written by  the  node  order  array.  In  this  case,  the  dimensions  of 
arrays  N0D0RD,  N0D0RDS,  NODPTR,  and  NODPTRS  can  be  increased  to 
200  in  the  main  program.  If  there  are  more  than  200  nodes  in  the 
section,  the  map  description  to  program  RCINPT  must  be  simplified 
or  the  problem  cannot  be  run. 
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If  the  number  of  steps  to  or  from  the  section  exceeds  100,  or 
if  the  total  number  of  steps  from  the  garage  or  landfill  to  the 
section  and  back  to  the  landfill  exceeds  200,  the  route  data  on 
file  TAPE9  will  be  incorrect.  No  error  message  is  given. 

The  number  of  steps  in  each  trip  piece  can  be  increased  to  200 
by  increasing  to  200  the  dimensions  of  arrays  ISFS,  ISTS,  NDFS, 
NDTS,  NODORD,  NODORDS,  NODPTR,  and  NODPTRS  in  the  main  program. 

The  total  number  of  steps  can  be  increased  to  600  by  increasing 
to  600  the  dimensions  of  arrays  ISSV,  NDSV,  and  CORT  in  the  main 
program.  Increases  beyond  these  numbers  of  steps  require  extensive 
program  modifications  and  should  be  considered  beyond  the  capa- 
bilities of  the  program. 


SECTION  VIII 
RECOMMENDED  CHANGE 


In  subroutine  DISCON,  an  error  message  will  not  print  node  numbers 
10000  to  99999  because  of  an  inadequate  field  width  specification.  The 
occurrence  of  14  in  format  45  should  be  changed  to  15. 
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I 


START 


Is 

array  length 
greater  than 
zero? 


Binary  search  array 
for  number. 


START 


Compute  pointer  separation. 


Compute  pointer  end  value. 


Set  initial  and  final  pointers. 


Save  values  at  inspection  point. 


final  location  value 
greater  than  or  equal  to 
initial  location 


Move  initial  values  to  final  location 


Set  final  pointer  equal  to 
initial  pointer. 


Move  initial  pointer  up 
by  spacing  value. 


Subroutine  SHLSRT3 
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Save  values  at  inspection  point. 


Move  initial  values  to  final  location. 


I 

Set  final  pointer  equal  to 
initial  pointer. 

' T 

Move  initial  pointer  up 
by  spacing  value. 


Subroutine  ISHLSRT 


Initialize  loop  on  nodes. 


' 


Subroutine  ADJUST 
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Is  R2  < Rl? 


the  first  segment " — - 
cause  the  second  segment's  vehicle 

capacity  to  be  

" — exceeded?^—- 


Add  the  first  segment's  refuse,  Rl , 
to  the  second  vehicle's  load. 


Remove  the  first  segment's  refuse 
from  the  first  vehicle’s  load. 


Reassign  the  first  segment  to 
the  second  segment's  section. 


Print  a segment 
reassignment  message,. 


Add  the  second  segment's  refuse,  R2, 
to  the  first  vehicle's  load. 


Remove  the  second  segment's  refuse 
from  the  second  vehicle's  load. 

. — — — 


Reassign  the  second  segment  to 
the  first  segment's  section. 


Subroutine  ADJUST 
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nmtnmtun* 


Print  piece  heading. 


Initialize  segment-scanning  loop. 


the  segment  in  the 
desired  piece? 


Increment  the  segment  count. 
Print  the  segment  number. 


Subroutine  PRNPCS 


Subroutine  TRACE 
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' 


START1 


Set  the  front  pointer,  JJ,  to  0. 


Set  the  end  pointer,  KK,  to  one  more  than  the  number 
of  nodes  in  the  node  order  and  pointer  arrays. 


Increment  the  front  pointer. 


Have 

the  pointers 
met? 


Does 

JJ  point  to  an  odd- 
order  node? 


Decrement  the  end  pointer. 


Have 

the  pointers 
met? 


Subroutine  MOVODD 


START 


Define  functions  for  section  number,  starting  and 
ending  node  numbers,  segment  length,  segment  travel 
time,  and  number  of  wavs  of  travel  on  segment. 


Compute  weighting  factor 
for  time. 


Use  I FIND  to  obtain  the  base 
node's  line  number. 


the  1 i ne  number 
valid? 


Print  an  error  message. 


Call  EXIT  to  terminate 
the  program. 


Subroutine  TREE 
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Initialize  the  loop  on  the 
levels  to  be  added  to  the  tree. 


Set  the  new  node  indicator  to  0. 


Move  right  to  the  highest-level 
node  in  the  tree. 


Has 

the  level  for  new 
additions  been 
reached? 


Subroutine  TREE 
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Initialize  the  loop  over  segments 
bounding  the  node. 


Find  a node  that  can  be  reached 
from  the  current  highest-level 
node. 


Compute  the  weighted  time-distance 
value  to  the  new  node. 


Has 

the  node  been 
used  before? 


the  new  time- 

distance  value  greater  than 
or  equal  to  the  old 


Disconnect  the  pointers  to  and 
from  the  node. 


Remove  the  node  from  the 
in-use  category. 


Connect  pointers  from  the 
previous  node  to  the 
subsequent  node,  if  any. 


Subroutine  TREE 


emove  from  the  in-use  category 


the  deleted  node 


Save  the  time  and  distance  to  the  node 


Save  the  weighted  time-distance  value  at  the  node 


Save  the  line  number  of  the  previous  node 


Save  the  number  of  the  previous  segment 


the  previous  node 
have  any  other  bounding 
segments? 


Seek  another  highest- level  node 


Subroutine  TREE 


there  only  one 
bounding  segment? 


Initialize  the  loop  over  the 
bounding  segments. 


Obtain  the  segment  number. 


/loopN 


the  segment  in  the 
current  section? 


Obtain  the  line  number  of  the 
node  at  the  other  end  of  the 
segment. 


Subroutine  C0N2ST 


2 


- — Are  ' 

the  segments  In 
the  same  piece  of  the 
current  section? 


the  two  segments  be 
traveled  In  at  least  one 
direction? 


Compute  the  total  length  of  the 
two-segment  path. 


Have 

any  two-step  paths 
been  saved  yet? 


— Has  ' 

a shorter  path  been 
saved  between  the  same 
two  pieces? 


Save  the  segment  numbers,  piece 
numbers,  and  path  length. 


Subroutine  C0N2ST 
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Have 

all  possible  second 
segments  been 
checked? 


Have 

all  possible  first 
segments  been 
checked? 


Have 

all  nodes  been 
checked? 


Change  the  piece  numbers  of  segments  in  the  higher-numbered 
piece  to  the  lower-numbered  piece  number. 


RETURN 


Subroutine  C0N2ST 
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— " Have ' 
any  nodes  been 
saved  from  the 
same  piece? 


the  time  to  the  new 
node  smaller? 


Select  a place  to  store  the  node. 


Save  the  node,  overwriting 
the  node  with  the  worst  time,  if  necessary. 


Find  the  worst  time  currently 
saved. 


Are 

any  more  segments 
to  be  checked? 


Save  paths  from  each  saved  node 
to  the  starting  node. 


RETURN 


Subroutine  FNUPTH 


Are 

there  any  more 
connections  between 
separate  pieces? 


Are 

there  any  disconnected 
pieces  left? 


RETURN 


Subroutine  CONNST 
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Find  the  next  segment  in  the  section. 


Save  the  starting  node  of  the  segment. 


Call  tRtE  to  build  paths  from 
the  starting  node  using  only 
segments  in  the  section. 


Call  TREE  to  build  paths  to 
the  starting  node  using  only 
segments  in  the  section. 


Increment  by  1000  the  section  number  of 
each  segment  in  the  section  that  is  not 
connected  to  the  starting  node. 


Were 

any  separate  pieces 
found? 


Subroutine  CONNECT 
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2 


Initialize  the  loop  over  segments. 


Does 

the  segment 
connect  two  pieces  of 
the  section? 


Merge  the  connecting  segment 
and  the  two  pieces  into  one 
piece. 


Are 

any  segments  left 
to  be  checked? 


3 


< 


there  only  one 
piece? 


Remove  from  the  section  any 
pieces  that  require  no 
refuse  collection. 


there  only  one 
piece  left? 


Call  CONNST  to  connect 
the  remaining  pieces. 


RETURN 


RETURN 


> 


Set  the  repeat  indicator  to  0 


Initialize  the  loop  over  the 
nodes  in  the  section. 


the  node's  order 
unequal  to  1? 


Find  the  segment  in  the  section 
that  the  node  is  connected  to. 


there  any  houses 
on  the  segment? 


Delete  the  segment  from  the 
section. 


Remove  the  node  from  the  section 


Subroutine  ISCON 


f Buffer  in  \ 
from  TAPE7 
paths  to  and 
from  the  garage 
v and  landfi 11.  . 


< 


Call  SHLSRT3  to  sort  the 
travel  times,  carrying  along  the 
node-pointer  and  node-order  arrays. 


> 


Subroutine  EPXP 
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Have 

any  nodes  been 
saved  yet? 


any  saved  nodes 
lie  on  the  path  from  this 
node  to  the  garage  or 
— landfill?  - 


Save  the  node  as  an  entry  or 
exit  node. 


Subroutine  EPXP 
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START 


thTTe^ionnumber^  - 

differ  from  that  on  the  previous 
cal  1 , or  have  no  paths 
_ been  saved?  - — "" 


Clear  storage  used  by  out-of-section 
path  segments. 


Set  the  count  of  artificial 
segments  to  0. 


Initialize  the  loop  on 
nodes  in  the  section. 


the  node's  order 
unequal  to  1? 


Get  the  neighboring  segment 
count. 


Subroutine  CL0SE1 
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Call  TREE  to  generate  paths 
back  to  the  section. 


Find  and  save  the  shortest  path 
back  to  the  section. 


Initialize  the  loop  over 
starting-point  nodes. 


t*  ' 

Call  TREE  to  generate 
paths  from  the  starting 
node. 

1 

Save  the  distance  from  each 
node  to  the  starting  node. 


NO 


Subroutine  GENDM 


Find  the  nearest  neighbor 
for  each  node. 


i 

Initialize  the  loop  on  nodes. 

i* ::  :■ 

Find  the  largest  (if  IDIR=-1 ) or  smallest  (if  ID  I R=+ 1 ) 
near-neighbor  distance  among  unpaired  nodes. 


Pair  the  nodes. 


Initialize  the  loop  on  nodes. 


Subroutine  SELORD 


the  node's  near 
neighbor  paired? 


Select  a new  near  neighbor 
for  the  node  from  the 
remaining  unpaired  nodes. 


any  nodes  left  to 
be  checked? 


Are 

any  nodes  left  to 
be  checked? 


RETURN 


Subroutine  SELORD 


Subroutine  ISCON 


START 


Set  the  profitable  cycle  indicator 
to  0. 


Initialize  the  loop  over  nodes. 


Move  the  diagonal  entries  of  the 
old  path  and  profit  matrices 
to  the  new  path  and  profit  matrices, 


Examine  all  cycles  starting  and 
stopping  at  the  node  for  a 
more  profitable  cycle. 


Was 

a more  profitable 
cycle  found? 


pYES 

Save  the  new 

profit  and  path. 

Subroutine  NEXTM 
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Compute  the  number  of  nodes 
that  can  be  packed  into  48  bits 


Compute  the  number  of  path-doubling 
iterations  necessary  to  produce  the 
longest  possible  path. 


Form  the  path  and  profit  matrices 
for  two-step  paths  between  nodes. 


Initialize  the  path-doubling  loop 


Call  NEXTM  to  find  the  profits  and  most 
profitable  paths  between  nodes  when  the 
maximum  number  of  steps  is  doubled. 


a profitable  cycle 
found? 


Subroutine  SOLV 


Has 

the  longest  possible 
path  been  examined? 


RETURN 


Increment  the  iteration  count 


fewer  than  20  iterations 
been  performed? 


Print  an  error  message 


Select  the  alternate  steps  in  the 
profitable  cycle  as  a new  odd- 
node  pairing. 


Subroutine  SOLV 


Subroutine  TRAVEL 


15? 


the  segment's 

starting  node  the  previous 
node  in  the  path? 


the  segment  a 
one-way  segment  or  is  its 
ending  node  not  the  previous 
— ^node  in  the  path? 


Append  the  segment  and  its 
starting  node  to  the  path. 


Decrement  the  segment's 
availabil ity  count. 


any  segments  left 
to  be  examined? 


Subroutine  TRAVEL 


Append  the  segment  and  its 
ending  node  to  the  path. 

Decrement  the  segment's 
availability  count. 

LH  ■ ~ 

2 


Get  the  line  number  of  the 
previous  segment  in  the  path. 


Can 

the  segment  be 
used  again? 


Append  the  segment  to  the  path 
segment  array. 


Append  the  next- to- last  node  in  the 
path  to  the  path  node  array. 


Decrement  the  segment's 
availability  count. 


Does 

the  current  path  end 
on  the  correct  node? 


Count  the  segments  remaining  to 
be  used. 


$ 


Subroutine  TRAVEL 
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RETURN 


Was 

a travel -only  segment 
found? 


Seek  the  next  collection  segment. 


Was 

a collection  segment 
found? 


Subroutine  OPTPATH 
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1 


there  more  than  one 
segment  in  the  traveled 
stretch? 


Call  TREE  to  find  the  minimum 
time  path  from  the  beginning 
to  the  end  of  the  traveled  segments. 


Are 

any  U-turns  present? 


the  U-turn  cause  

an  appreciable  saving  of 
distance? 


FEWER 


new  path  fewer, _ 

the  same,  or  more  segments 
than  the  original 


Subroutine  OPTPATH 
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/ Read  \ 
file  TAPE3 
for  pointers  to 
the  segments  in 
each  section. 


Print  the  title  and  the  heading  for 
the  initial  section  assignments. 


Initialize  the  loop  over  sections 


f Read  \ 
file  TAPE4  > 
for  the  numbers 
of  the  segments 
in  the  section . , 


Print  the  section  and  segment 
numbers  


Program  PHASE3 
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1 


Are 

any  sections  left? 


<Call  I 
order. 


ISHLSRT  to  put  the  segments  into  numeri 
. Carry  along  the  section  assignments. 


HCa,\ 


Compute  the  vehicle  load  for 
each  section. 


Count  the  segments  bounding 
each  node. 


Read  the  node  numbers  of  the 
landfill  and  garage. 


Print  the  program-initiated 

section-change  heading. 

< 


Call  ADJUST  to  reassign 
segments  bounding  order-two 
nodes  into  the  same  section 


Print  the  new  vehicle 
s. 


Read  tne  user-specified 
section  changes. 


Print  the  user-specified 
section  changes.  


Program  PHASE3 


Print  the  segment  data 


Call  TREE  to  generate 

minimum  time  paths  to  all 

nodes  from  the  landfill  and  garage 


< 


Call  TRACE  to  obtain  the 
path  from  the  garage  to  the 
landfill. 


> 


/ Buffer  out  \ 
/ to  TAPE7  the  \ 
/data  for  paths  to\ 
» and  from  the  I 
V landfill  and  J 
\ garage. 


Print  the  node  numbers  of  the  landfill  and  garage. 


Print  the  data  describing  paths  to  and  from  the 
landfill  and  garage.  


Initialize  the  loop  over  nodes. 


Does 

the  node  connect  to 
the  landfill? 


Print  an  error  message. 


Program  PHASE3 


^nod 


'use  function  I FIND  and  subroutine  M0VE3  to  save  the 
node  numbers  and  pointers  for  the  nodes  in  the  sectic 


< 


Call  DISCON  to  delete  dangling 
segments  that  require  no 
collection 


Call  EPXP  to  obtain  choices 
for  entry  and  exit  nodes. 


Save  the  node  orders  and  pointers. 


he^v. 

tion^ 


> 


> 


Does 

the  entry  node  number 
equal  zero? 


Initialize  the  loop  on  the  exit  nodes. 


Does 

the  exit  node  number 
equal  zero? 


Restore  the  node  orders  and  pointers. 


Increment  the  orders  of  the  entry 
and  exit  nodes. 


Set  to  1 the  number  of  traversals 

for  each  segment  in  the  section 

or  to  0 for  each  segment  not  in  the  section. 


/loop\ 


Program  PHASE3 
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170 


Jl 


Call  SELORD  to  generate  an 
initial  pairing  of  odd  nodes. 


Are 


Are 

any  odd-node  pairs 
left  to  be  examined? 


Remove  pairs  of  traversals  from 
any  two-way  segment  traversed 
more  than  twice. 


^^Call  TRAVEL  to  obtain  a travel  path  over  all  \ 
segments  in  the  section,  starting  at  the  entry  node  and 
ending  at  the  exit  node. 


Mas 

a path  found? 


Print  a warning 
message.  - 


Mas 

an  error  found? 


Call  OPTPATH  to 
improve  the  path. 


Call  TRACE  to  obtain  the  path 
from  the  section  to  the  landfill 


Use  CUMDIS  to  obtain  the  distances 
for  the  three  parts  of  the  trip. 


the  current  trip 
distance  greater  than  or  equal 
to  the  saved  distance? 


Save  the  new 

trip  distance. 

Save  the 

new  path. 

Are 

there  any  more 
exit  node  choices? 


Program  PHASE3 


APPENDIX  B 
PROGRAM  LISTINGS 


Subroutine  FLIP 
Subroutine  M0VE3 
Function  IFIND 
Subroutine  SHLSRT3 
Subroutine  ISHLSRT 
Subroutine  ADJUST 
Function  CUMDIS 
Subroutine  PRNPCS 
Subroutine  TRACE 
Subroutine  MOVODD 
Subroutine  TREE 
Subroutine  C0N2ST 
Subroutine  FNDPTH 
Subroutine  CONNST 
Subroutine  CONNECT 
Subroutine  DISCON 
Subroutine  EPXP 
Subroutine  CL0SE1 
Subroutine  GENDM 
Subroutine  SELORD 
Subroutine  PATH 
Subroutine  NEXTM 
Subroutine  SOLV 
Subroutine  TRAVEL 
Subroutine  OPTPATH 
Program  PHASE3 
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APPENDIX  C 

DEFINITIONS  OF  IMPORTANT  VARIABLES 
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Page 

Subroutine  FLIP 

240 

Subroutine  MOVE3 

240 

Function  IFIND 

240 

Subroutine  SHLSRT3 

240 

Subroutine  ISHLSRT 

240 

Subroutine  ADJUST 

241 

Function  CUMDIS 

241 

Subroutine  PRNPCS 

241 

Subroutine  TRACE 

242 

Subroutine  MOVODD 

242 

Subroutine  TREE 

242 

Subroutine  C0N2ST 

243 

Subroutine  FNDPTH 

243 

Subroutine  CONNST 

244 

Subroutine  CONNECT 

244 

Subroutine  D I SCON 

245 

Subroutine  EPXP 

245 

Subroutine  CLOSE  1 

245 

Subroutine  GENDM 

246 

Subroutine  SELORD 

247 

Subroutine  PATH 

247 

Subroutine  NEXTM 

247 

Subroutine  SOLV 

248 

Subroutine  TRAVEL 

248 

Subroutine  OPTPATH 

249 

Program  PHASE3 

249 

239 


Note:  A single  variable  symbol  may  have  different  meanings  in  relation  to  the 
various  subroutines.  For  this  reason,  variables  are  defined  below  for 
each  subroutine  and  for  program  PHASE3. 


SUBROUTINE  FLIP 


Al,  A2,  A3 
IF 


A count  of  the  items  to  be  reversed. 
The  array  to  be  reversed. 


SUBROUTINE  MOVE 3 

Arrays  to  be  moved. 

Subscript  data  are  moved  to. 
Subscript  data  come  from. 


IARRAY 


FUNCTION  I FIND 


Array  being  searched. 
Length  of  IARRAY. 
Number  being  sought. 


SUBROUTINE  SHLSRT3 

Arrays  reordered  as  array  X is  sorted. 

Count  of  words  to  be  sorted. 

If  SGN  = 1.0,  X is  sorted  into  increasing  order. 
= -1.0,  X is  sorted  into  decreasing  order. 

Array  to  be  sorted. 


SUBROUTINE  ISHLSRT 

Array  reordered  as  array  IX  is  sorted. 

The  spacing  between  words  to  be  sorted. 

If  I SGN  =1,  IX  is  sorted  into  increasing  order. 
= -1,  IX  is  sorted  into  decreasing  order. 


If  SGN 


If  I SGN 
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IX 

NW 


ISTG 

KNODES 

NBRSEG 

NH 

NRQF 

NSC1 , NSC 2 
NSECT 

NSG1,  NSG2 

Rl,  R2 

STG 

TC 

TL 


ISEG 

ISTG 

NSG 

STG 


ISTG 

ITRIP 

MSECT 

NPIECE 

NSEG 

NTRIP 


Array  to  be  sorted. 

Count  of  words  to  be  sorted. 

SUBROUTINE  ADJUST 

Array  of  integer- valued  segment  data. 

Count  of  nodes. 

Array  of  packed  neighboring-segment  numbers. 

Symbolic  subscript  of  ISTG  array  for  house  count. 

Symbolic  subscript  of  STG  array  for  refuse-quantity 
adjustment  factor. 

Section  assignment  numbers. 

Symbolic  subscript  of  ISTG  array  for  section  assignment. 
Segment  numbers. 

Refuse  quantities. 

Arrays  of  floating  point-valued  segment  data. 

Array  of  vehicle  capacities. 

Array  of  vehicle  loads. 

FUNCTION  CUMDIS 

Array  of  segment  numbers. 

Array  of  integer-valued  segment  data. 

Count  of  segments. 

Array  of  floating  point-valued  segment  data. 

SUBROUTINE  PRNPCS 

Array  of  integer-valued  segment  data. 

Piece  number. 

Symbolic  subscript  of  ISTG  array  for  temporary  section 
assignment. 

Count  of  pieces  whose  segments  will  be  printed. 

Count  of  segments  In  map  description. 

Section  number. 


T 


KNODES 

LEND 

LPREVN 

NODNUM 

NPREVS 


NSTART 

NSTOP 


SUBROUTINE  TRACE 

Array  to  be  filled  with  node  numbers  in  path. 
Array  to  be  filled  with  segment  numbers  in  path. 
Count  of  nodes  in  map  description. 

Line  number  of  terminal  node  of  path. 

Array  of  line  numbers  of  next  node  in  path. 
Array  of  node  numbers. 

Array  of  numbers  of  next  segment  in  path. 

Count  of  segments  in  path. 

Node  number  of  start  of  path. 

Node  number  of  end  of  path. 


NODORD 

NODPTR 


SUBROUTINE  MOVODD 

Count  of  nodes  in  section. 

Count  of  odd  nodes  in  section. 

Array  of  orders  of  nodes  in  section. 
Array  of  pointers  to  nodes  in  section. 


INSECT 


KNODES 

LA 

MAXLVL 


SUBROUTINE  TREE 

The  weighted  time  and  distance  value  from  a node  to  the 
starting  node. 

An  jrray  of  distances  from  the  nodes  to  the  starting  node. 

If  zero,  all  segments  may  be  used;  otherwise,  only  segments 
in  section  INSECT  are  used. 

An  array  of  four  pointers  for  each  element  in  the  tree  (node 
pointer,  left-branch  pointer,  right-branch  pointer,  and  back 
pointer) . 

Array  of  Integer-valued  segment  data. 

Count  of  nodes  in  map  description. 

Back-pointer  value. 

Maximum  number  of  steps  allowed  away  from  the  starting  node. 

Variable  that  indicates  whether  any  new  nodes  have  been  added 
to  the  tree. 

Line  number  of  next  available  free  storage  in  the  1ST  array. 
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NIU  Array  of  weighted  time  and  distance  values  from  nodes  to  the 

starting  node. 

NIULOC  NIULOC(I)  gives  the  line  number  in  the  1ST  array  of  node 

NODNUM(I). 

NODNUM  Array  of  node  numbers. 

NPREVN  Array  of  line  numbers  of  next  node  in  path  to  starting  node. 

NPREVS  Array  of  segment  numbers  of  next  segment  in  path  to  starting 

node. 

ST6  Array  of  floating  point-valued  segment  data. 

TIME  Array  of  travel  times  from  starting  node  to  each  node. 

W Coefficient  of  distance  in  weighted  time  and  distance  value. 


DIST 

ISEG 

ISTG 

JSEG 

KNOOES 

MSECT 

NBRSEG 

NCON 

NP1 , NP2 
NS 

NSEG 

NS1 , NS2 
NTRIP 
NUMNBR 
STG 


SUBROUTINE  C0N2ST 

Length  of  two-segment  path. 

Number  of  first  segment  in  two-segment  path. 

Array  of  integer- valued  segment  data. 

Number  of  second  segment  in  two-segment  path. 

Count  of  nodes  in  map  description. 

Symbolic  subscript  of  ISTG  array  for  temporary  section 
assignment. 

Array  of  packed  neighboring-segment  numbers. 

Count  of  pieces  connected  by  two-step  paths. 

Arrays  of  piece  numbers. 

Count  of  entries  in  NP1  and  NP2  arrays. 

Count  of  segments  in  map  description. 

Arrays  of  segment  numbers  of  segments  in  two-step  paths. 
Section  number. 

Array  of  counts  of  neighboring  segments. 

Array  of  floating  point-valued  segment  data. 


SUBROUTINE  FNDPTH 

IP  Number  of  piece  containing  node  NORG. 

IPATH  Double-subscripted  array  to  be  filled  with  segment  numbers  of 

up  to  five  paths  connecting  pieces  of  section  NTRIP  to  node 


ISTG 

KNODES 

LNF 

LPREVN 

NDSV 

NODNUM 

NPF 

NPREVS 

NTRIP 

TIM 


IPATH 

ISTG 

LNF 

MSECT 

NDSV 

NORG 

NPF 

NPIECE 

NSEG 

NTRIP 

TIM 


ISTG 

ITRIP 

KNODES 

MSECT 

NTRIP 


Array  of  integer-valued  segment  data. 

Count  of  nodes  in  map  description. 

Array  of  line  numbers  of  nodes  starting  paths  to  node  NORG. 
Array  of  line  numbers  of  next  node  in  path. 

Count  of  paths  saved  in  array  IPATH. 

Array  of  node  numbers. 

Array  of  piece  numbers  from  which  paths  start. 

Array  of  numbers  of  next  segment  in  path. 

Section  number. 

Array  of  travel  times  for  paths. 


SUBROUTINE  CONNST 

Two-dimensional  array  of  segment  numbers  in  paths  to  node  NORG. 

Array  of  integer-valued  segment  data. 

Array  of  line  numbers  of  nodes  from  which  paths  in  array 
IPATH  start. 

Symbolic  subscript  of  ISTG  array  for  temporary  section 
assignment. 

Count  of  paths  ending  at  node  NORG. 

The  number  of  a node  in  the  current  section. 

An  array  of  piece  numbers  from  which  paths  to  node  NORG  start. 
Count  of  pieces  of  section  NTRIP. 

Count  of  segments  in  map  description. 

Section  number. 

Array  of  travel  times  for  paths  in  array  IPATH. 


SUBROUTINE  CONNECT 

Array  of  integer-valued  segment  data. 

Piece  number. 

Count  of  nodes  in  map  description. 

Symbolic  subscript  of  ISTG  array  for  piece  number. 
Section  number. 
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SUBROUTINE  DISCON 

ISTG 

Array  of  integer-valued  segment  data. 

KN 

Count  of  nodes  in  section. 

NHRSLG 

Array  of  packed  neighboring-segment  numbers. 

NODNUM 

Array  of  node  numbers. 

NODORD 

Array  of  orders  of  nodes  in  the  current  section 

NOOPTR 

Array  of  pointers  to  nodes  in  current  section. 

NTRIP 

Section  number. 

SUBROUTINE  EPXP 

DISTS 

Array  of  travel  times  to  nodes  from  the  landfill  or  garage. 

I 

Section  number. 

KLIM 

Maximum  number  of  entry-  or  exit-point  nodes 

to  be  saved. 

KN 

Count  of  nodes  in  section  I. 

NNTF 

Two-dimensional  array  of  line  numbers  of  the 
path  to  or  from  the  landfill  or  garage. 

next  node  in  the 

NODORD 

Array  of  orders  of  nodes  in  the  section. 

NODPTR 

Array  of  pointers  to  nodes  in  the  section. 

SUBROUTINE  CL0SE1 

DIST 

Array  of  distances  from  nodes  to  the  order-one  node. 

IDIR 

IDIR  = 1 for  travel  away  from  the  starting  node  or  -1  for 
travel  toward  the  starting  node. 

IPATHN 

Array  of  line  numbers  of  nodes  in  the  path  from  the  section 
to  the  order-one  node. 

I PATHS 

Array  of  segment  numbers  in  the  path  from  the  section  to  the 
order-one  node. 

ISTG 

Array  of  integer-valued  segment  data. 

KN 

Count  of  nodes  in  the  section  before  CL0SE1  is  called. 

KNODES 

Count  of  nodes  in  the  map  description. 

KNX 

Count  of  nodes  in  the  section  after  CL0SE1  is  called. 

LPREVN 

Array  of  line  numbers  of  the  next  or  previous  node  in  the 
path  to  or  from  the  order-one  node. 

MAXLVL 

The  maximum  number  of  steps  allowed  in  the  path  from  the 
order-one  node  to  the  section. 

MSECT 


NODNUM 


NODORD 


NODPTR 

NORG 


NOTHER 


NPREVS 


NTRIP 


KNODES 


NODNUM 


NODPTR 


NTRIP 


ONDMTX 


Symbolic  subscript  of  the  ISTG  array  for  temporary  section 
assignment. 


Word  of  packed  neighboring-segment  numbers. 
Array  of  node  numbers. 

Array  of  orders  of  nodes  in  the  section. 
Array  of  pointers  to  nodes  in  the  section. 
The  node  number  of  an  order-one  node. 


Number  of  the  node  at  the  other  end  of  the  segment  from  the 
order-one  node. 


Array  of  numbers  of  the  next  or  previous  segment  in  the  path 
to  or  from  the  order-one  node. 


Count  of  segments  in  the  path  from  an  order-one  node  to  the 
section. 


Symbolic  subscript  of  ISTG  array  for  number  of  times  the 
segment  is  traversed. 

Count  of  segments  saved  in  IPSSV  array. 

Section  number. 


Array  of  travel  times  between  nodes  and  the  order-one 
node. 


SUBROUTINE  GENDM 


Array  of  distances  from  node  NORG. 

Count  of  nodes  in  map  description. 

Array  of  line  numbers  of  next  node  in  path  to  node  NORG. 
Count  of  odd  nodes. 

Array  of  node  numbers. 

Array  of  pointers  to  node  numbers. 

Number  of  node  used  as  the  start  of  a minimum-distance  tree. 


Array  of  the  number  of  the  next  segment  in  the  path  to  node 
NORG. 


Section  number. 

Matrix  of  distances  between  odd  nodes. 


Array  of  travel  times  from  nodes  to  node  NORG. 
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SUBROUTINE  SELORD 

DIST  Matrix  of  distances  between  odd  nodes. 

IDIR  If  IDIR  = 1,  the  nodes  with  the  smallest  near-neighbor  distance 

are  paired  first;  if  IDIR  = -1,  the  nodes  with  the  largest  near- 
neighbor distance  are  paired  first. 

IPART  Array  of  sequence  numbers  of  other  nodes  in  odd-node  pairs. 

NNN  Dimension  of  distance  matrix. 

NOOD  Count  of  odd  nodes. 

SUBROUTINE  PATH 

II  Sequence  number  of  starting  node  of  first  path. 

IPART  Array  of  sequence  numbers  of  other  nodes  in  odd-node  pairs. 

ISTP  Array  of  sequence  numbers  of  nodes  in  path. 

ISTPN  Matrix  of  packed  sequence  numbers  of  nodes  in  new  path. 

ISTPO  Matrix  of  packed  sequence  numbers  of  nodes  in  old  path. 

JJ  Sequence  number  of  node  ending  second  path. 

KK  Sequence  number  of  node  ending  first  path  and  starting  second 

path. 

MID  Count  of  nodes  through  the  end  of  the  first  path. 

NEDGE  Count  of  odd  nodes. 

NG  Error  indicator  returned  by  subroutine  PATH. 

NNPW  Maximum  number  of  sequence  numbers  that  can  be  packed  into 

one  computer  word. 

SUBROUTINE  NEXTM 

ISTP  Array  of  sequence  numbers  of  nodes  in  the  current  path. 

ISTPN  Matrix  of  sequence  numbers  of  nodes  in  the  new  paths  between 

odd  nodes. 

ISTPO  Matrix  of  sequence  numbers  of  nodes  in  old  paths  between 

odd  nodes. 

JSTP  Array  of  sequence  numbers  of  nodes  in  the  most  profitable  cycle. 

MTXN  Matrix  of  profits  in  new  paths  between  odd  nodes. 

MTXO  Matrix  of  profits  in  old  paths  between  odd  nodes. 

NEDGE  Count  of  odd  nodes. 

NG  Error  indicator. 

NRES  Profitable-cycle  indicator. 


SUBROUTINE  SOLV 


DIST 

Matrix  of  distances  between  odd  nodes. 

IPART 

Array  of  sequence  numbers  of  other  odd  nodes  in  odd-node  pairs 

ISTPI 

, ISTP2 

Matrices  of  packed  sequence  numbers  of  nodes  in  paths  between 
odd  nodes. 

JSTP 

Array  of  sequence  numbers  of  odd  nodes  in  profitable  cycle. 

MTX1 , 

MTX2 

Matrices  of  profits  in  paths  between  odd  nodes. 

NCH 

Iteration  count. 

NNPW 

Maximum  number  of  node  sequence  numbers  that  can  be  packed 
into  a word. 

NODU 

Count  of  odd  nodes. 

NRES 

Profitable-cycle  indicator. 

SUBROUTINE  TRAVEL 

IERR 

Error  indicator. 

IPN 

Array  of  numbers  of  nodes  in  path. 

IPS 

Array  of  numbers  of  segments  in  path. 

ISP 

Number  of  node  at  which  path  ends. 

1ST 

Number  of  node  at  which  path  starts. 

KLINE 

Array  of  line  numbers  of  segments  in  path. 

KS 

Array  of  backup  counts. 

KTIMES 

Array  of  counts  of  times  segments  are  available  for 

traversal 

MAXKS 

Maximum  backup  count  allowed  before  path  search  is 

terminated 

NEXT  1 1 

Array  of  line  number  to  be  used  next  in  adding  a segment  to  a 
particular  step  of  the  path. 

NN 

Count  of  nodes  in  path. 

NN1 , NN2 

Arrays  of  numbers  of  nodes  starting  and  ending  the 

segments . 

NSEG 

Count  of  segments  in  section. 

NSG 

Array  of  segment  numbers. 

NTIMES 

Array  of  number  of  times  each  segment  must  be  used. 

NWAY 

Array  giving  the  number  of  ways  a segment  may  be  traversed. 
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ON 

DO 

IDEL 

IS 

ISC 

ISTG 

LEN 

NH 

NNP 

NS 

NSC 

NSP 

STG 

TRTY 


CORT 

DFD 

DFG 

DIST 

DTD 

DTF 

IERR 

IPAIR 

IPN 

IPS 

ISFS 

ISSV 

ISTART 

ISTG 


SUBROUTINE  OPTPATH 

Length  of  new  travel  stretch. 

Length  of  old  travel  stretch. 

Difference  between  new  and  old  counts  of  steps  in  travel 
stretches. 

Count  of  steps  in  old  travel  stretch. 

Section  number. 

Array  of  integer-valued  segment  data. 

Symbolic  subscript  of  STG  array  for  segment  length. 

Symbolic  subscript  of  ISTG  array  for  house  count. 

Array  of  node  numbers  in  path. 

Count  of  segments  in  new  travel  stretch. 

Symbolic  subscript  of  ISTG  array  for  section  assignment. 

Array  of  numbers  of  segments  in  path. 

Array  of  floating  point-valued  segment  data. 

Array  of  travel -or-col lection  indicators  for  each  segment  in 
path. 

PROGRAM  PHASE 3 

Array  of  travel -or-col lection  indicators  for  path  saved  for 
output. 

Array  of  distances  from  landfill. 

Array  of  distances  from  garage. 

Total  length  of  a trip. 

Array  of  distances  to  landfill. 

Matrix  of  distances  to  and  from  garage  and  landfill  for  entry- 
and  exit-point  nodes. 

Error  indicator. 

Array  of  sequence  numbers  of  other  nodes  in  odd-node  pairs. 
Array  of  numbers  of  nodes  in  travel  path. 

Array  of  numbers  of  segments  in  travel  path. 

Array  of  numbers  of  segments  in  path  from  section  to  landfill. 

Array  of  numbers  of  segments  saved  for  output  in  final  travel 
path. 

Number  of  entry-point  node. 

Array  of  integer-valued  segment  data. 
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ISTOP 

ISTS 

JMAX 

KMAX 

KN 

KNODES 

MAXODD 

MSECT 

NBRSEG 

NDFS 

NDMP 

NDSV 

NDTS 

NGAR 

NH 

NN1 , NN2 

NODD 

NODNUM 

NODORD 

NODPTR 

NSECT 

NSEG 

NTRIP 

NUMNBR 

ONDMTX 

STG 

TC 

TLOAD 

TRTY 


Number  of  exit-point  node. 

Array  of  numbers  of  segments  in  path  from  garage  to  section. 

Maximum  number  of  entry-point  nodes  to  be  used. 

Maximum  number  of  exit-point  nodes  to  be  used. 

Count  of  nodes  in  section. 

Count  of  nodes  in  map  description. 

Maximum  number  of  odd  nodes  that  can  be  paired. 

Symbolic  subscript  of  ISTG  array  for  temporary  section 
assignment. 

Array  of  packed  bounding-segment  numbers. 

Array  of  numbers  of  nodes  in  path  from  section  to  landfill. 
Node  number  of  the  landfill. 

Array  of  numbers  of  nodes  in  travel  path  saved  for  output. 
Array  of  numbers  of  nodes  in  path  from  garage  to  section. 

Node  number  of  the  garage. 

Symbolic  subscript  of  array  ISTG  for  house  count. 

Symbolic  subscripts  of  ISTG  array  for  starting  and  ending  node 
numbers. 

Count  of  odd  nodes. 

Array  of  node  numbers. 

Array  of  orders  of  nodes  in  section. 

Array  of  pointers  to  nodes  in  section. 

Symbolic  subscript  of  ISTG  array  for  permanent  section 
assignment. 

Count  of  segments  in  map  description. 

Count  of  sections. 

Array  of  counts  of  neighboring  segments. 

Matrix  of  distances  between  odd  nodes. 

Array  of  floating  point-valued  segment  data. 

Array  of  vehicle  capacities. 

Array  of  vehicle  loads. 

Array  of  travel -or-col lection  indicators  for  current  travel 
path. 
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GLOSSARY 


Air  Force  Refuse-Collection  Scheduling  Program:  a set  of  four  computer  pro- 
grams that  perform  residential  refuse-collection  scheduling  and 
produce  printed  schedules  and  maps  of  the  routes. 

binary  tree:  a data  structure  in  which  each  element  can  be  linked  to  at  most 
two  other  elements,  referred  to  as  left  and  right  branches. 

entry  node:  a node  in  a collection  region  at  which  a path  from  the  garage  or 
landfill  reaches  the  section. 

exit  node:  a node  in  a collection  region  at  which  a path  to  the  landfill 
starts. 

node:  a numbered  point  on  a street  at  which  some  characteristic  of  the  street 
changes. 

odd  node:  a node  at  which  an  odd  number  of  segments  meet. 

profitable  cycle:  a cyclical  path  of  paired  odd  nodes  that  reduces  the  total 
pairing  length  if  new  pairs  are  formed  from  the  consecutive  unpaired 
nodes . 

section:  a group  of  segments  serviced  by  the  same  collection-vehicle  trip. 

segment:  a portion  of  a street  between  two  nodes. 
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